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ABSTRACT 
The purpose of this study was to investigate the kinetics of the reverse Claus 
reaction, 


2H Ones) Xeoxe——-> 211.5 +.SO;, 


This was done by obtaining space-time/conversion data at various partial pressures of 
sulphur and water, from which initial rates of the reaction are calculated, based on 
octatonic sulphur. 

Kaiser S-201 catalyst, (primarily gamma-alumina) was used in an Autoclave Berty 
(internal recycle) type reactor. Liquid sulphur and water were vapourised and mixed with 
nitrogen as a carrier gas to form the feed stream. The reaction products, water and 
nitrogen were analysed on a sulphur—free basis by gas chromatography and an electronic 
integrator was used to measure the peak areas. 

A total of 173 steady state runs were performed, yielding 163 finite rate 
measurements and 18 initial rates. Ten initial rates were calculated at a constant water 
partial pressure of 3720 Pa and at partial pressures of sulphur varying from 266 Pa to 
2955 Pa, which was the full range that was possible with the experimental apparatus. 
The remaining eight initial rates were obtained at a constant sulphur partial pressure of 
756 Pa, and a variable water partial pressure, from 1537 Pato 21,443 Pa The balance 
of the feed streams was always nitrogen. For each set of experimental runs, the space 
velocities ranged from 50 to 220 hr-! and the pressure and temperaturewere held 
constant at 1.46 atm (1.479 * 10° Pa) and 592 K respectively. Two more initial rates 
were obtained at temperatures of 624 K and 653 K, employing an additional 10 finite 
rate measurements. 

A rate function of the form, 


a b #) 
-Rw = A [exp(-E/R T)] (Pw) (Ps,) /[1 + K, Pwl 


where, 
—Rw = rate of reaction with respect to water, mol/hr g cat, 


Pw = partial pressure of water, Pa, 
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Ps, = partial pressure of S,, Pa, 


was found to fit the data) The following values of the parameters were obtained, 


A =0.5004 
-E/R = -7994 
aati 089 
by i=20:39 


K, 2.742 * 10-5 


The temperature is in Kelvin, and the partial pressures are in Pa. This rate equation 
was found to be consistent with earlier work by McGregor (197 1) and Liu (1978), who 
studied the forward reaction rate at the University of Alberta. The form of the rate 
equation is the same as that suggested by them, and the values of the parameter K, 
obtained in this work are similar. The initial rate of the reverse reaction is less by a 
factor of about 10 than that of the forward reaction, which agrees with the results 
obtained by them for the forward reaction. However, a reaction rate equation 
incorporating both a forward and reverse rate term using all the finite rate data obtained 
was not developed. A direct comparison of the two expressions was not therefore 
made. 

A total of 52 finite rate data at high catalyst loadings (10.0066 g) were used to 
estimate the maximum conversions obtainable in a Berty type reactor (CSTR). The 
resulting experimental conversions were found to be greater than the predicted 
thermodynamic conversions at the same temperature. Some possible explanations for 


this discrepancy were examined. 
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1. INTRODUCTION 


1.1 Background 

Sulphuric acid is unchallenged as the world's leading industrial chemical (26,57). 
Indeed, its total annual production may be used as a measure of the industrial prosperity 
of a nation (26). It is not then surprising that, of the fifty million tons of sulphur 
produced annually (1976), the great majority is used in the production of sulphuric acid 
(57). It follows that sulphur production is a particularly important facet of the chemical 
industry. 

During the eighteenth and nineteenth centuries, the elemental sulphur deposits 
found in Sicily were the main source of sulphur. Towards the end of the nineteenth 
century, however, the Frasch process (after its inventor, Hermann Frasch) became almost 
the sole source of sulphur. This process produced cheap, high purity sulphur from the 
unique dome deposits of Louisiana and Texas. This sulphur was of such high quality and 
low cost that, even when shipped across the Atlantic to Europe, it was cheaper than the 
conventionally mined Sicilian product (57). 

At about the same time that Frasch was granted a patent for his process, a British 
inventor, C.F. Claus, obtained one for a process that entailed the burning of hydrogen 
sulphide in air over bauxite to produce elemental sulphur and water (17). This process 
only became industrially significant as a major source of elemental sulphur from 1960 
onwards, with the development of the natural gas industry at Lacg in France and in 
Alberta, Canada (57). Since 1970, however, half of all the sulphur produced worldwide 
annually has been recovered sulphur (from hydrocarbon processing), obtained mostly by 
means of the Claus process (57). 

The natural gas reserves of Lacq and Alberta are designated as being sour’. That 
is to say, they contain significant quantities of hydrogen sulphide —— in some cases as 
much as forty percent. This is a highly toxic and corrosive gas that readily oxidises to 
sulphur dioxide, which is only slightly less poisonous (57). Yet sulphur dioxide is more 
hazardous because it is soluble in water and forms sulphurous acid. The acid rain of 
environmental concern is essentially sulphurous acid (61). Clearly, sour gas is unsuitable 


either as a fuel or as a feedstock to other chemical processes. The Claus process has 
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proved to be a successful method of “sweetening” natural gas, particularly at the 
wellhead, so that the gas may be safely transported by pipeline to a processing plant. As 
a result, the process has become a major means of producing raw sulphur even though, 
in essence, this element is merely a byproduct of the necessary gas sweetening process. 
There have been several modifications to the original process as patented by 
Claus (57). These have been reviewed in the next chapter. The most important 
modification was the splitting of the reaction into two independent steps: in the first, the 
hydrogen sulphide was partially burned in a waste heat boiler and in the second, the 
sulphur dioxide produced was reacted with the remainder of the hydrogen sulphide in a 
Catalytic reactor. By this means, the large heat of reaction was dissipated without the 
damage to the catalyst which had been a problem in previous versions of the process. 


The process can be summarised by the following chemical reactions : 


GiGi 2) ete) SS See Hee (1.1) 


2.H;Sigit SOMole => 22:H,Olg)+ 18 /81S;(o] muna ee eee (1.2) 
AH°.9, = 108.2 kJ 


Sulphur vapour has been known for many years (4), to exist as several molecular 
species containing different numbers of sulphur atoms The basis of equation 1.2 above 
is octatonic sulphur. In actuality, the sulphur species may range from monotonic to 
octatonic and even higher, depending on the temperature and partial pressure of sulphur. 
As a consequence, sulphur vapour is made up of a distribution of these species. At low 
pressure, the vapour is almost entirely S, while at high temperatures, S, predominates. 
At intermediate temperatures, the distribution of species may be predicted by the Gibbs 
free energy of the system. McGregor (56) developed a computer program whose 
function was to minimise the Gibbs free energy of the sulphur system. This was revised 
by Truong (81) and also by Razzaghi (72) to include the "best available” thermodynamic 
data and also other species encountered in the Claus reaction. Many workers have 


reported that the equilibria predicted using thermodynamic data do not correspond to 
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those achieved experimentally (28,53,56). This has been considered to be due to 
imprecise thermodynamic data (56). The experimental difficulties in obtaining accurate 
data on the various sulphur species are great. The data on the odd numbered sulphur 
species in particular, are subject to large discrepancies depending on the source. In 
order to minimise the effect of such inaccuracies, the present work will consider the 
sulphur species S,, S,, S,, and S, to be present in the reacting system. There is a much 
higher degree of unanimity between sources on thermodynamic data on these species 


(81). The Claus reaction then takes the following form: 


PS PUSe Ne aCe y eu aH@n Wee ge bare e (1.3) 
2A OE GOW Re ey (1.4) 
Ca ee ee eae a (1.5) 
De SO Se ee ee ee (1.6) 


Thus, if the kinetics of the reaction are determined via equation (1.3), the equations (1.4) 
to (1.6) are assumed to be in equilibrium. The odd number sulphur species have been 
ignored, since the uncertainty in the free energy data of these species is much higher 
than that associated with the even number species. Using the data that are available, it is 
found that the odd numbered species are a small percentage of the total sulphur (72). An 
illustration of the difference between the two distributions is given in Appendix 2. The 
sulphur equilibrium reactions are assumed to equilibrate rapidly and reaction (1.3) is 
assumed to be rate-determining. It should be noted that, while historically the overall 


reaction : 


ASee) 2005 => Tt ONG) tins (ine le kere meena eres re (122) 
AH?®,,, = 628.3 kJ/mol 


has been called the Claus reaction, it is reaction (1.2), the catalytic reaction itself, which 
has come to be known as the Claus reaction, although it should more properly be 
referred to as the modified Claus reaction. Since reaction (1.2) is sulphur yield limiting, it 


shall be referred to as the Claus reaction in this work. 
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Further improvements to the Claus process have had the aim of increasing the 
sulphur yield.) For example, increased yield has been achieved by means of improved 
boiler design (20,42), H,S burner design (25), and control of the H,S/air ratio (9). The 
development of improved catalysts (20,66) and the use of more optimal operating 
conditions in the catalytic reactor have also improved the yield. 

Stricter regulations governing the permitted levels of sulphur emission have 
forced further improvements to sulphur recovery plants (57,65,66). Since impurities 
such as ammonia (from amine treatment piants) or hydrocarbons (for example, methane 
from natural gas) can significantly reduce the conversion of hydrogen sulphide, the 
pre-treatment of the feed stream for their removal is required. This can also greatly 
increase the concentration of hydrogen sulphide (16,30,31). More important than these 
improvements, though, has been the development of Claus tail gas clean-up units to 
process the Claus effluent stream even further (33,65,76,88). Such units are now used 
extensively and can increase the sulphur recovery from approximately 97% to over 99%. 

In an exothermic reaction, the equilibrium yield is favoured by low temperatures. 
Decreasing the temperature also increases the sulphur molecularity towards eight. This 
will reduce the partial pressure of sulphur in the product stream and as a result, increase 
the equilibrium yield. Obviously, therefore, it is advantageous to operate the Claus 
catalytic reactor at as low a temperature as possible. Modern developments have been 
directed at taking advantage of this (32,33), even to the extent of operating the Claus 
converters at a temperature below the dew point of sulphur, as in the Cold Bed 


Absorption (CBA) process. 


1.2 Purpose of this Study 

The equilibrium of the catalytic reaction (1.2) has been found to be the limiting 
factor governing the yield in the Claus reaction. There is a significant discrepancy 
between the experimental equilibrium conversions obtained by several workers (4,23,28) 
and the conversions predicted by the "best available” thermodynamic data It would be 
useful to obtain empirical data on the forward and reverse rates of the catalytic reaction 
in order that the equilibrium may be more accurately evaluated since at equilibrium, the 


forward and reverse rates are equal. With this additional knowledge, the design and 
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control of Claus plants may be made more precise (20). 
The rate of an equilibrium reaction such as the Claus reaction will be regarded as 


the addition of two separate rates representing the forward and backward reactions, viz., 


In the particular case of the Claus reaction, the forward reaction has been the subject of 
much research and this has been surveyed in Chapter (2). However, there has been no 
study of the reverse Claus reaction published, although it has been considered briefly in 
some forward reaction studies. 

The purpose of this study is to determine the kinetics of the reverse Claus 


reaction : 
Zin sOr F316eS5 St ASSO eee eee ae (1.9) 


The kinetics of the forward reaction have been satisfactorily determined and the rate 


equation is of the form (22) : 


1 Y 
(Ph, s) 


—-Rso, =k (Pso,) 
 Cpiaaeiame oe oe ee ee (1.10) 


It is the intention of this study to seek a comparable expression for the reverse reaction. 
If the rate expression is to exhibit thermodynamic consistency with equation (1.10) 


above, a form such as equation (1.11) is suggested : 


a b 
—Rw = k, (Pw) _(Ps) 
1RKNRWH 0 NRE yee oe (1st) 


The above equations apply individually, only to the special case of the initial rate 
of reaction, where no products are present to take part in the appropriate reverse 


reaction. The forward and reverse rate expressions may be combined in a general 
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equation of the following form, which represents the overall rate of the Claus reaction : 


fe d a b 
RW =k (Phos). (PS0,) = Kediewin les.) 


eae | 1 Ky Cw meee Tee ae 112) 


The parameters K, and K, will be equal since they represent the adsorption equilibrium 
constant of water and hence will be a function of the temperature only. 


This empirical expression may be useful in the design and control of Claus plants. 


1.3 Experimental Approach 

An equilibrium reaction is perhaps most easily investigated by using the method of 
initial rates. This approach was employed by McGregor (56) to investigate the forward 
Claus reaction. The initial rate is that rate of reaction which occurs in the absence of any 
products, that is: 


a b 
Se RWiSsKeP WES | ON esis See eee ne eee (ets) 


By obtaining finite rate data at various space velocities, it is possible to derive the 
initial rate. The data is extrapolated to obtain the rate at zero conversion. The theoretical 
analysis is outlined in Chapter four. Once the initial rates have been determined, the 
parameters of the rate equation may be derived by applying a combination of linear and 
non-linear regression analyses to an assumed form of the rate expression. There are 
three effects studied in this work : The variation of the rate of the reverse reaction with 
(1) the partial pressure of sulphur; (2) the variation of the rate with partial pressure of 


water; (3) the effect of temperature on the rate. 
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2. LITERATURE SURVEY 


2.1 Introduction 

Although the literature contains many references to the Claus process, these are 
in the main of a descriptive rather than a theoretical nature. This survey will concentrate 
on that section of the literature which is pertinent to the kinetics of the Claus reaction 
itself. Other aspects of the Claus process have been extensively reviewed elsewhere, 
notably by Truong (8 1), Liu (53), Cho (15), Liu (52), and McGregor (56). Other reviews are 
available for the period up to 1970 by Pilgrim and Ingraham, (69) and for 1970-75 by 
Chandler (13). 


2.2 The Claus Reaction 
The reaction, 


ripe) se CyP2 (OF, emer e Oleic se WWI or cs ee (2.1) 


had been well known as a laboratory demonstration reaction from the early nineteenth 
century (18), but it was C.F. Claus who, in the 1880's, first designed a commercial 
process using it (17). His patent of 1887 describes a process for the recovery of 
elemental sulphur from calcium sulphide. The sulphide, when treated with carbonic acid, 
produces hydrogen sulphide which subsequently reacts with air over one of several 
Catalysts, including bauxite, to produce sulphur. 

In 1887, the Chance brothers (12) were granted a patent for a commercial 
process for the recovery of elemental sulphur from the waste gases of the Le Blanc 
process (57). This process for producing sodium carbonate was responsible for a great 
deal of the pollution in nineteenth century England during the height of the Industrial 
Revolution (26). The Chance application of the Claus process represents its first use in 
the field of pollution abatement. 

Reaction (2.1) is highly exothermic, however, and the Chance—Claus process 
proved very difficult to control (69). As a consequence, in 1937, |.G. Farbenindustrie 


modified the process and effectively eliminated the problem of temperature control (2). 
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This was achieved by burning one-third of the hydrogen sulphide with air in a waste heat 
boiler, and producing sulphur dioxide. This sulphur dioxide was then reacted with the 
remainder of the hydrogen sulphide to produce the final products, elemental sulphur and 


water. Hence the original reaction takes place in two stages: 


RS) Cae NO meer NO) Peal OP UG RCLOIO NC) ES tye (222) 


Padme SO err Oa Cyp ese PASO IE) ess eon (2.3) 
A =O Ce a OU KoS 


The major fraction of the total heat of reaction is thus dissipated in the furnace section. 
This process became known as the modified Claus process, and is the basis of all 


present Claus plants. 


2.3 Claus Reaction Kinetics 

The adoption of the Claus process as a primary method of desulphurisation of 
natural gas and, as a result, of sulphur production, has prompted research on the kinetics 
of the Claus reaction. 

Pilgrim and Ingraham (69) reviewed the process and its kinetics up to 1970, while 
Chandler (13) reviewed subsequent developments between 1970 and 1975. The 
emphasis, however, of these reviews is on the practical aspects of the process itself. 

The first reference to the reaction in the literature is by Cluzel (18) in 1812. He 
noted that no sulphur would be produced if the gases were first dried by passing them 
over calcium chloride. Lewis et al., (49,51,70) studied the equilibrium of the H,S/SO, 
system in a series of papers. Their work, published in 1918, demonstrated that the 
reaction is reversible and also that it is heterogeneous. The rate of the reaction was 
reported as being high in both directions. An equilibrium constant was calculated, but this 
was subject to large deviations. 


1 The heats of reaction are given as ranges since they represent the actual process ~ 
reactions. The temperatures and pressures are not therefore precise. The distribution 
of sulphur species will also be subject to variation. 
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The first quantitative work was performed by Taylor and Wesley in 1927 (79). 
This work confirmed the heterogeneity of the reaction. Despite reported difficulties in 
analysing the reaction products, and in particular, determining the sulphur partial pressure, 
they ventured an estimate of the rate equation for the reaction on glass, which they 
reported to be of the form: 


Seite Sea (tis S)2! (9 SO.) eee ee (2.4) 


The authors discussed what the nature of the reaction mechanism was likely to be but 
this was essentially speculation since, at that time, there was a dearth of experimental 
data on the adsorption of the various species on glass. | 

In 1937, Undintseva and Chufarov reported (83,84) that hydrogen sulphide and 
sulphur dioxide would not react homogeneously at temperatures between 250 and 
350 °C. They also found glass, aluminium, and aluminium oxide to be good catalysts. In 
addition, iron and iron oxide were reported to exert only a weak catalytic effect. 

Murthy and Rao (62), in 1952, studied the reaction over metallic sulphide catalysts 
using a batch recycle reactor at 25 °C. The condensation of sulphur on the catalyst 
prevented steady state from being attained since the surface area of the catalyst was 
changing. They reported that, at this temperature, the reaction would not proceed with 
any of the catalysts tested unless water was present. 

The work of Gamson and Elkins, reported in 1953, represents the first 
comprehensive study of the Claus reaction (28). The reaction was studied in an integral 
reactor over Porocel catalyst (an activated bauxite). Space time versus conversion data 
_ were obtained but no attempt was made to derive a rate equation. A hypothetical sour 
gas plant feed stream, comprising 6.75% H,S, 3.39% SO,, 26.83% H,0, and the balance, 
N,, by volume, was used. Four space velocities, and three temperatures ranging from 
230 to 300 °C, were employed. The authors were careful to caution that the kinetic data 
of this study were inconsistent with their thermodynamic analysis of the Claus reaction in 
that the experimental conversions were higher than those calculated from free energy 
data. 

In 1957, Hammar (38) studied the reactions involved in a sulphur plant technology, 


using a glass differential reactor. He failed, as had his predecessors, to derive a rate 
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equation for the Claus reaction. 

Cormode, in 1965, studied the Claus reaction using a recycle flow reactor 
containing Porocel catalysts (19). He suggested that H,S itself was decomposing in his 
apparatus but the analysis was not precise enough to allow him to come to a firm 
conclusion. 

Munro and Masdin (60) reported, in 1967, that water significantly retarded the 
rate of reaction. This effect was attributed to Le Chatelier's principle. They employed 
activated alumina and molecular sieve as catalyst, over the temperature range of 
100-300 °C at a space velocity of 70 vol H,S/cat vol hr. The water content in the feed 
was varied from 5 to 20 %. 

Landau et al., (47) concluded in their results, published in 1968, that diffusional 
effects were rate limiting. Unfortunately, they did not state whether pore or film 
diffusion was occurring and, in addition, they did not describe their experimental 
equipment nor the catalyst charge sizes. It is impossible, therefore, to determine what 
type of diffusion, if any, was occurring. 

A report by the Allied Chemical Corporation in 1970 (1) on Claus process kinetics 
contains a comprehensive literature survey. A dynamic flow system and Porocel LPD 
catalyst were used. The experimental conditions were varied over the following ranges : 
the temperatures were between 200 and 370 °C, residence times between 0.125 and 
2.50 seconds, and feed concentration of H,S between 1 and 6 %. It was reported that 
the rate of reaction was so fast that equilibrium was attained within 0.5 seconds. The 
observed equilibrium was found to be higher than that which had been predicted using a 
mathematical model that considered the sulphur vapour to be composed solely of S,. 


The following rate constant for the reaction was derived ; 


EG 1GG.expis0b ia Gene ae eae (2.5) 


assuming the reaction to be first order in all reactants. The temperature was measured in 
Kelvin, but pressure units were not given and a full rate equation was not developed. 
The work of Struck et al., (74) was carried out at low temperatures, below the 


sulphur dew point. The reaction was found to go almost to completion, in accordance 
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with the predicted thermodynamic equilibrium at these temperatures. At 182 °C, the 
achieved conversion was higher than that predicted, and at 100 °C, the conversion was 
close to 100%. 

The first serious attempt to derive a rate equation for the Claus reaction was 
made by McGregor in 1971 (56), as part of an extensive programme of research being 
carried out at the University of Alberta. In a series of 80 experiments involving the use 
of an external recycle reactor containing Porocel catalyst, the reaction was studied by 
means of initial rate experiments over the range 208-287 °C. The following parameters 


were derived : 


Ea = 7589 + 451 cal 
0.828 + 0.0951 


order wrt H,S 


0.417+0.111 


order wrt SO, 


Various mechanistic equations were fitted to the data but no attempt was made to 
discriminate between the models. Water was observed to exhibit a barely perceptible 
autocatalytic effect at low partial pressures, and a marked retardation at high partial 
pressures. The catalytic activity was found to vary with the external area of the catalyst 
pellets, and it was concluded that the reaction took place primarily on the exterior 
surface of the catalyst. Finally, no evidence of a homogeneous component to the rate of 
reaction was found. 

By applying non-linear data reduction techniques to all the available data, an 


empirical rate expression incorporating a reverse reaction term was obtained : 


—Rso, = 1.292 exp[-—7 100/RT] [(Ph,s)%9°* (Pso,)475 — 0.504 (Pw) 98] en) 


However, as the variances were large, there was a need for more data. 
The empirical expressions derived by McGregor were found to correlate the data 
as well as any of the mechanistic models he proposed. The results reported by Hammar 


and by Cormode were compared with it. The expressions of Hammar (38) correlated the 


= 


12890 = BESO = BH rw ssbie oon 
Pht 2S FlhO = Oe tweed. 


re . a 
ch stint aa Agata -ccinad ttt ty aN BAMA re encetaaRe aPunannee MOEN, 


sichtaso ws vie ee midi gram see MAW, Bigcaie ortt eaevevtedt ‘raid 
isiaG Mig somncawestaity dita wol re taste scan 
raylatcn ort tos sae avy e oto | file ine Give Stems AAT peruano | ; 
concn sept Test gar isesd caw 1 8 eaea 
iq-ahes. set ot manognoy iC TOTO tage on viet Seylatep are to emiwe | 
Srka* Baily TRY 

maaia’ atal> en 7 


[hog SEED a Sa) ae sas.t — 


- 7 ; j 7 eal raed —_ mee = 
> ues mew ieee) 
— oy ae na ~~ a <— : _ 


ae 
ar - eh ee Por * Pw 


Zz 


a 


data well, although the actual data could not be compared. The work of Cormode (19) 
was in reasonable agreement, though the actual rates were higher. 

In 1972, Karren (39) was unable to obtain steady state reaction data because he 
found that the liquid sulphur collected by the sulphur condenser catalysed the reaction. 
The retarding effect of water was noted, as was the extent of sulphur loading on the 
Catalyst. 

Dalla Lana et al., in 1972, took the data analyses of McGregor further (22), and 
concluded that, on statistical evidence, the reverse reaction term of equation (2.5) could 
not be justified. They postulated that the retarding influence of water was more likely to 
be due to adsorption of water on the catalyst. They derived a different equation using 


MceGregor's data : 


-Rh,s = Ae 1.121 exp[-7440/RTI (Ph,s) (Pso,)°5 
Ap [is+0:0042 sible (2.7) 


The term Ae/Ap was introduced by McGregor to allow for variations in specific surface 
area between different catalysts or even between different forms of the same catalyst. 
The orders of the reaction with respect to H,S and SO, were rounded off to agree with 
the stoichiometry of the reaction. 

In 1978, George (30) performed an initial rate study of the Claus reaction using 
cobalt molybdate catalyst. A fixed bed integral reactor was used. A temperature range 
of 240-450 °C and pressures of up to 1.3 atmospheres were employed. The following 


rate equation in the form of equation (2.6) above was developed : 


ERsor= & (Phisi/ 1120.10 PWi... wee. (2.8) 


The actual experimental orders were : 
ness) eee MSIE) SOS Is) 
SOM O2e OOS 
The activation energy was found to be: 


Ea = 5.5 kcal/mol 
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lt was concluded that sulphur did not catalyse the reaction as had been suggested 
by Hammar (38). Water again retarded the reaction and this effect was allowed for by 
the adsorption term, (the denominator) in equation (2.7). In addition, film diffusion effects 
were found to be negligible and it was suggested that pore diffusion was significant. 

From his experiments in 1975 to evaluate a copper—doped catalyst Cho (15) 
concluded that conversions greater than the theoretical equilibrium were obtained. The 
effect of external concentration gradients was negligible. The reverse reaction was 
thought to be significant although its visible rate was negligible. 

Dalla Lana et al. (21) derived, in 1976, arate equation for the Claus reaction using 


data obtained by Liu in his work on the pure gamma alumina catalyst. 


—Rso, = 0.460 exp[-7350/RT] (Ph,s) (Pso,) ° 
1 2 


hes 006 Pwiiae © yee eee ere (2.9) 


This equation is nearly identical to that derived from McGregor's study (22). The similarity 
between equations (2.6) and (2.8) was remarkable since each was derived from data on 
two different catalysts. 

The likely mechanism of the reaction was investigated and a 
Langmuir—Hinshelwood rate expression was derived, based on the following surface 


reaction step which was considered irreversible and rate controlling : 
H,S.AI* + O.AI* —>H,OAIy + SAI =... (2.10) 


where Al” is a lewis acid site on the catalyst surface. 


The postulated mechanism led to the following rate equation : 


Y2 


0.5/n Z 
-Rso, = k (Ph,s) (Pso,) / [Psn/K,] [1 + Pw/K,] 


where n, the average size of the sulphur molecule, is approximately 7. The sulphur term 


is then negligible and the expression is similar to equation (2.8). 
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Kerr et al. in 1976, reported that they had found a second order equation to apply 
in their experiments involving bauxite and activated alumina catalysts (42). This was 


reported as being, 


d(H,S)/dt = [FESHISON + ps fe! toni ee (2.12) 
where, 

= 24149 exp[ -2526/T ] 
and, 


Ea = 5.02 kcal/mol 


The experiments were carried out in a macroscopic vessel under isothermal conditions 
with flow rates, temperatures, compositions, and pressures corresponding to those 
attained in the second, third, or fourth catalytic converters of an industrial plant. 

In 1977, Blanc et al. (8) developed an essentially empirical reaction model which 


included terms to account for the reverse reaction. The equation was given without 


units : 
=soy — k expi-E/RT] ~ (Ph,s) (Pso,)*= — (1/K%) (Pw) (Ps)? 
[Te Pwil te Veer ae © ee. (2.13) 
where, 
K = (Pw)? (Ps) 
(Ph,s)? (Pso,) 


This equation was applicable at temperatures greater than 270 °C. The power of the 
denominator term could be varied over a large range without the sum of the squared 
deviations between the experimental and the computed values changing greatly. This was 
explained by the fact that the factor multiplying the water adsorption term, and the 
power of the denominator are very highly correlated. 

The catalyst used by Blanc et al. was activated alumina, type 5-10, and the data 
were collected over the temperature range 270-370 °C. Feed concentrations of H;S 
ranged from 2.5 to 10 %, while the ratio of SO, to H,S ranged from 0.4 to 0.6. The 
amount of water in the feed was varied from 5-35%. General results were as follows : 


the overall order is greater than 1 ; the order with respect to sulphur is not zero ; the 
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overall rate increases with temperature ; the activation energy is small ; water vapour 
favours the reverse reaction and decreases the forward reaction ; the reaction is very 
fast ; there are no film diffusion effects ; and pore diffusion effects are apparent. 

In 1978, Liu (53) studied the Claus reaction as catalysed by gamma-alumina (Alon), 
using a recycle reactor and employing infra-red analysis to measure the reaction rate. 


The following rate equation was derived : 


-Rso, = by, exp{-E,./RT ] (Ph,s) (Pso,) °5 
su EB. cexpE,/RT) Pail PS no EXPI-E,/RT) (PWIP (PS, = we ee (2.14) 


where, 
Digs 7a Oe 
Deceit #21 
E,/R = 3509 + 254 
eryisi = Shileys) 


This equation is very similar to that developed by McGregor, who used the Porocel 
(bauxite) catalyst. 

Evidence of the reverse reaction was sought by introducing Claus reaction 
products. Since hydrogen sulphide and sulphur dioxide were not subsequently detected, 
it was concluded that the reverse reaction was negligible. This test was performed at 
210 and 320 °C. The nitrogen flow used was 40 ml/min. with saturated sulphur at 
300 °C. Water was introduced at 1.95 * 10-* gmol/min. The model developed by 
McGregor was found to correlate this data well. The term (Ph,s ee ) was found to 
be vital to a good fit. Water was observed to have a significant retarding effect on the 
reaction, as did sulphur. Film diffusion was negligible and it was concluded that pore 
diffusion could be disregarded providing the catalyst weight was below 250 mg. 

In 1978, Grancher (34,35) reported on the work performed by the Aquitaine 
research group at Lacq in France. A rate expression was determined, but it was not 
reported in these papers. The author reported that the expression had the following 


features: the order with respect to H,S was 1, and with respect to SO,, 0.5; and there 
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was inhibition due to water. The reaction was reported to be limited by internal diffusion 
resistances. The validity of the model was checked by comparison with the Lacq plant 
itself. 

Finally, the paper by Liu et al., published in 1980 (50), summarises the state of the 
art with respect to kinetics studies of the Claus reaction. The authors conclude that use 
of arate equation is an aid in making process calculations and that these are "sensitive” to 
the various rate expressions. In particular, the reliability of these rate expressions Is 
increased when they are based on surface chemistry studies. There is still, however, no 
reliable method for predicting either the thermodynamic properties or the dynamics of 
the sulphur vapour equilibrium, and this is a liability. Liu'’s final conclusion ts that the 
particular catalyst used is merely a matter of preference, unless special conditions 
otherwise dictate. For example, the presence of hydrocarbons in the feed stream may 
introduce significant quantities of carbon sulphides, CS, and COS. If this is the case, it 
would be advantageous to employ a catalyst which promoted the decomposition of these 


species in addition to catalysing the Claus reaction. 
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3. EQUIPMENT DESCRIPTION AND OPERATION 


3.1 Introduction 

This experimental study is essentially similar to that performed by McGregor (56) 
in his study of the forward reaction in that an initial rate study was made of the Claus 
reaction system. There are, however, several important differences. In particular, the 
reactants, water and sulphur, are not vapours at room temperature and pressure. As a 
result, the feed and reactor are different. Several different versions of the experimental 
equipment were employed in an effort to determine which, with respect to smoothness 
of flow, ideal mixing, et cetera, was the most satisfactory. The most successful is 
shown in Figure 3.1 below. The salient features of some of the alternative systems 
employed are discussed in the final section of this chapter. This is the system used to 
obtain the rate data of this study. Another significant difference between the McGregor 
study and this one is the use of a Berty—type internal recycle reactor. The external 
recycle reactor of McGregor’s work has several disadvantages for the water/sulphur 
reaction system. The merits of the Berty reactor vis-a-vis the external recycle reactor 


will be discussed in detail. 


3.2 Experimental System 

A diagrammatic representation of the experimental system is shown in Figure 3.1. 
It shows a stream of dry nitrogen which is pre-heated to 375 °C in an oven[11]*. Liquid 
sulphur from a syringe pump [10] is added in the same preheat oven[11]. The 
temperature is above the sulphur dew point and the sulphur vapourises. The water 
component is added by a syringe pump [12] and is heated to the temperature of the feed 
line and vapourised. It is then added to the feed stream immediately before it enters the 
reactor [13]. The reactor effluent is passed to a sulphur condenser [14] which eliminates 
most of the unreacted sulphur. The product stream is then split, with a small portion of it 
being led to the gas chromatograph by way of a final sulphur condenser [15], and then is 


expelled through a vent to the atmosphere. 


?The numbers in square brackets following an item of equipment refer to Figure 3.1. 
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Figure 3.1 Experimental System 


Equipment List 


—_ 


Nitrogen cylinder 


7h Drier 

Sf Rotameter 

4. Needle vaive 

5 Differential flow controller 

6. Feed pressure and temperature sensors 
The Differential pressure cell 

8. Feed control valve 

Sh Reactor pressure sensor 

10. Sulphur syringe 

11. Preheat Oven 

12. Water syringe feeder 

13, Berty reactor 

14. Sulphur condenser 

15: G.C. secondary condenser 

16. Secondary sulphur condenser 
17, Water condenser 

18. Reactor pressure control valve 
19. Sulphur Syringe Reservoir 

20. Sulphur Syringe Oven 
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The main stream continues through to a final sulphur condenser [16] and a water 
condenser [17] before being vented via a control valve [18] to the atmosphere. 
The system had several drawbacks and some advantages which will be elaborated 


upon in the individual sections below. 


3.3 Nitrogen Feed System 

Linde Gas Products high purity nitrogen was regulated to 60 Ib/in? [1], a pressure 
which was adequate for the entire range of flow rates employed. Each new nitrogen 
cylinder was analysed for oxygen using the gas chromatograph. A molecular sieve 
column was used at room temperature. Every nitrogen cylinder was found to be free of 
oxygen contamination using this method. A single peak was the only response, indicating 
that nitrogen only was present. The nitrogen was passed through a 12 inch bed of 
calcium sulphate [2] in order to remove any water that was present. A small amount of 
"drierite’ [calcium sulphate coloured with cobalt] was used as an indicator of the 
efficiency of this operation. It was found unnecessary to replace the calcium sulphate. 
The stream was then passed through a Matheson rotameter [3] of tube size 602. The 
valve of the rotameter was kept fully open and the rotameter served only as a visual 
indicator of the flow. Immediately after the rotameter, a Nupro needle valve [4MG] was 
employed [4]. This valve provided the fine metering of the nitrogen flow. A Moore 
differential flow controller [63SD-—-L] was placed downstream of the needle valve to 
smooth any small fluctuation in flow [5]. This type of controller provides constant mass 
flow independent of pressure changes. 

The temperature and pressure of the nitrogen feed stream were measured [6] 
immediately before the differential pressure cell [7] which measured the volumetric flow. 
Temperature was measured using an iron-constantan thermocouple, and pressure by 
means of a USG pressure gauge, which was accurate to 0.1 psig. The volumetric flow 
rate was measured by a Foxboro differential pressure cell [613 DL] and controlled by a 
Foxboro controller and monitored by a Foxboro chart recorder. Volumetric flow rates 
of between 0-2 litres per minute at a pressure of 8.5 psig were obtainable. 

The nitrogen stream was then fed into the preheat oven [11] where it was passed 


through a 12 foot coil of 1/8 inch stainless steel tubing and was preheated to 375 °C 


1 7 ; 
_ ein 25) @nras 
U > <4 beds iio. Nig 46 ane oe the) Gots 
> \iz es * 2 ° ~— a ee - = > : 
a 


amet enisy te pe? C'S ug «| & een Obey 


7 oe ~~ oe Seta een GM Tee UoYES we 


bv emrt «df of.Ghas" et ect non a uraenta ks trys Paks be 
on eer oP an ar 2 wo wim n bo- te)! 207 Gee ; Mgt be 
ot : +] ro ) hitbeabuy Sraotue! (Ae Ree ins retowes GF 
ma 8 - or _ s p¥eenat OFr 19 ev Gin ier-gue wee Bs 
oe 4 et t bacw / taiea 40 Seer 0S sng oe 
eta = #i) €34~w1--F <2 pecbtird Seas t mts pect ra ar 
of: =f re sa to | Sloe? seated § rege 2: heres 777 ans rap’? 


ee ee eS a ee oot ere he evinn 


Ss _— a we “ 
of “4 fas We stetetover? Vorte i\eruoer “ni “aul @? ta 3 
See vt | , «< 2 4 


noe My = peppers si hi Petvos peuvext” {=i peyog 


emus & : 
’ "vt o's 


) @etey etree * meen eeighi sé: wate <= | Seta orale On teore 7 


bil (eepcrne i> silt eet Is a virlt 2 «cuatro Lone ee nage! 
eae vol Ss Pures 10 re apie al 


} 
“we adh x80 me pine Coe swear — : 


eaq trarst= 
in bawesre™ 4, 5 
avegit onnarAP OY of Se e-ciee! repetace FY hg 


oom nena 2a oma 
heat guides 


peep igre nt ated at 
rqites’ tierce tad 


before it was led to the sulphur mixing vessel. 


3.4 Sulphur Feed System 

The sulphur component of the feed was added by means of a liquid syringe 
feeder driven [10] by a constant speed motor . The liquid sulphur was led into the 
preheat oven [11] which was maintained at a temperature of 375 °C, where it entered a 
mixing vessel via a syringe-type nozzle and contacted the preheated nitrogen. This 
mixture then entered a surge tank to allow it to equilibrate. 

The cylinder of the sulphur syringe had an internal diameter of 1 inch and a 
working length of 5 inches. The entire syringe [10], including the reservoir [19], was 
enclosed in a temperature-—controlled oven [20] that was maintained at 140 °C. This 
temperature was chosen to ensure that no plugging occurred at the oven exit and also to 
maintain the liquid sulphur in the low-viscosity state. Sulphur undergoes molecular 
restructuring at about 159 °C and there is a dramatic increase in viscosity (80). 

The sulphur syringe was filled from a liquid sulphur reservoir [19]. The reservoir 
was sealed from the syringe by a valve. A second valve was inserted between the 
syringe and the mixing oven. This arrangement allowed an on-stream refilling of the 
syringe which disturbed the steady-state conditions only momentarily. This arrangement 
also ensured that no air would enter the syringe while it was being filled with sulphur. 
The seal of the syringe piston was made using two teflon "O-rings. The material of 
construction was 316 stainless steel. The reservoir valve was a Whitey ORM2, rated to 
232 °C. The syringe valve was a Whitey SS-NBS4-HTP, rated to 316 °C. 

The syringe pump itself was driven by a constant speed motor. This drove a gear 
train which could be altered to supply any one of six different speeds. The rotary 
motion of the gears was translated into horizontal displacement of a push rod by means 


of a worm gear. For the arrangement of the gears, see Figure A—1.1 in Appendix 1. 
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3.4.1 Operation 

Over a period of time, it was found necessary to replace the teflon "O-rings 
because they tended to wear. Great care was taken to ensure that no air remained in the 
syringe after the replacement of the "O-rings was complete. This was achieved by 
repeatedly filling the syringe with sulphur and then expelling it into the reservoir. This 
was done a minimum of three times. An added precaution was taken by filling the 
reservoir only as required, shortly before refilling the syringe itself. Additional problems 
were encountered due to the reservoir valve leaking and this required replacement on 


two occasions. 


3.5 Water Feed System 

A Sage Instruments model 355 syringe pump [12] was used to add tne water 
component of the feed stream. A 20 mi syringe was used for all runs. This size was 
sufficient to supply water for at least one run, and yet was small enough that the back 
pressure was not so great that it slowed the pump. In contrast with the sulphur feed, the 
water flow could be adjusted fractionally over a wide range of flow rates (v/de 
Appendix 1)>. The water was vapourised by heating the final 10 cm of the water line to 
the temperature of the reactor feed line, about 325 °C. The water was then added to the 
nitrogen and sulphur stream at its point of entry into the reactor [13]. 

The mixing of the two streams was achieved using a Swagelok "T" fitting. A 
section of 1/16th inch tubing was used to introduce the vapourised water directly into 
the centre of the main feed stream via a syringe-type nozzle. Perfect mixing of the feed 
stream was ensured by operating the reactor at high speed, at or above 1200 rpm. This 
was checked by conversion versus speed tests, (v/de Preliminary Experiments, 


Section 4.2). 


3.5.1 Operation 

The syringe pump itself was found to perform satisfactorily in all ranges. There 
was a need, however, to keep the water flowing long after the sulphur flow had been 
stopped to prevent the sulphur back diffusing up and plugging the 1/16th inch tubing 


which served to deliver the vapourised water. 


tetew of tee Gers are bec i R8e igheee epee ange & ” 
“a gad giegeebien OF A tremite bet? otto manaaenog/ "al 
ed ett Int Qs Varre cov. tay tonal amnieataeal 2 ond sate (Mique of mnensitiue | a 
ott isen usta aot ty 3 it Aarlipbe avons Aida dee ah eee oeageng Te 
st iNhegin? Qh" G agian é ed VighoTah tt bipe-vGa a7 eae, WR ; 

Rice ken stow act ool oe 
en <img ca eudewoqnetadlt 7 
tot jp ares Ai Sv NY 
irae aw a le omen ot -_ - 

A as Sr ' anaes oo wary is mares : 


PO a ee A ew 


ot well (ies of (oe edh ele 
ent os fehé Gert aqvy whey aT 2° Se 
nett qrapevient ole 9 

& grat doigeke «gna Om 

ers wird “aber inst oncrtcs rt te 
bed ame ra.griin Metts auaan ag 
ev wip ONS! ave ye SET Np 


a Nes ote + 
per net 


_ . = 
ie 


_ 


3.6 Reactor 

There are a wide choice of reactor types available to the researcher in the field 
of chemical reaction kinetics (10,73). The reasons why a particular reactor is chosen 
depend on a large number of factors including the nature of the reaction, the type of 
data required, and the reaction regime that is of interest. Several reviews have been 
published to aid in the discrimination between the various reactors (3,43,75). Carberry 
(11) suggests that the “ideal” laboratory reactor for the study of catalytic kinetics is one 
which operates isothermally over a wide range of conversions at steady state and with a 
well-established residence time, while facilitating the direct determination of the rate law. 
Accepting these criteria, it appears that the continuous stirred tank reactor has the most 
suitable characteristics. As a consequence, much effort has been devoted to designing 
laboratory reactors which approach continuous stirred tank (CSTR) operation, and which 
will also accommodate a heterogeneous catalytic reaction. Kiperman (44), credits Temkin 
with first proposing the term "gradientless reactors” for these reactors, since the quality 
of being well mixed ensures the elimination of both temperature and intraphase 


concentration gradients. 


3.6.1 Gradientless Reactors 

The use of gradientiess reactors in particular, has been reviewed by Sunderland 
(75) and Kiperman (43). The first reactor of this type was the external recycle reactor in 
which a pump recirculates a fraction of the reactants through a packed bed. The studies 
of McGregor (56) and Liu (52) at this university were performed using a reactor of this 
type. Butt et al. (9) and Perkins and Rase (68) discuss the use of this particular type of 
reactor. 

The external recycle reactor has several disadvantages. The external tubing and 
pump represent considerable dead volume compared to the volume of the packed bed. 
This results in long residence times which may be undesirable, especially in the kinetic 
study of a fast reaction. Also, depending on the temperature of the reaction, it may be 
necessary to cool the recycle stream before it reaches the pump and reheat it before it 
enters the reactor. This would be awkward, especially if there were a danger of one of 


the components condensing at a cold spot in the recycle loop, as in the case of the 
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partial oxidation of methanol to formaldehyde (27). In this reaction, the 
paraformaldehyde may condense (27). Obviously, this may be a problem with reverse 
Claus reaction studies, since not only are high temperatures involved ( > 600 kK), but 
significant sulphur partial pressures are also employed. Sulphur is notorious for its 
tendency to condense at inconvenient places. 

In order to avoid the limitations of the external recycle reactor, internal recycle 
reactors were developed. The first such was the so-called Carberry, or spinning basket 
reactor (Tajbl et al, 78). This reactor is similar to the simple CSTR except that the agitator 
is constructed to form a basket in which the catalyst is contained. Such a reactor was 
designed and built in this department by Turner (82) who also described its operation. 
The main disadvantage of this reactor is that it is impossible to arrange the catalyst and 
baffles, et cetera, to ensure that each catalyst pellet is equally accessible to the reactant 
phase since there are inevitably velocity gradients in such a mixed vessel. 

The inherent problems of both the external recycle and the spinning basket 
reactors may be eliminated by using an internal recycle reactor. There are many reactors 
of this type, and most have the common feature that the vapour has a definite path for 
circulation around the reactor vessel, during which it passes through a packed bed. 
These reactors have been reviewed by Bennet et al. (3) and by Sunderland (75). Each 
variation on the internal recycle theme has its own particular merits and demerits. 

The reactor developed by Berty et al. (6) was adopted for this study since this 
reactor is well documented (5,6,55) and, most importantly, the spinning basket reactor 
(SBR) developed by Turner (82) could be adapted to the Berty reactor configuration 
relatively easily. An important advantage of the Berty reactor, although not used in this 
study, is that it is possible to measure, or calculate, the actual mass velocities through the 
catalyst bed and thus reproduce, if desired, industrial reactor conditions . The mass 
velocities may be calculated by directly measuring the pressure drop across the catalyst 
bed and using Leva's or some similar correlation. This was the method used by Berty 


et al. (6) to evaluate the performance of this type of reactor. 
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3.7 Berty Reactor 

The spinning basket reactor of Turner was converted to a Berty reactor by 
extending the reactor casing and fitting internal baffles and an impeller. Extensive 
experimentation using a full size Lucite model determined that a forward curved bladed 
impeller performed best, though high speeds of over 2500 rpm were required to 
achieve visible flow through the packed bed even at a low catalyst loadings (1-2 g). To 
achieve the required flow pattern, it was found that the clearance of the impeller had to 
be as low as possible so that bypassing of the catalyst bed was reduced to a minimum. 
This testing was performed using a hot wire anenometer with smoke from hot solder 
flux, as an indicator of flow. The Turner SBR system could be operated at speeds of up 
to 4000 rpm and it was anticipated that this would be well in excess of that required to 
ensure idea! Berty reactor operation. The residence time distribution of the reactor was 
assessed by using a He pulse and checking the response with a gas chromatograph. This 
test indicated that ideal mixing was achieved at speeds as low as 500 rpm. A subsequent 
test, however, in which conversion was plotted against impeller speed, indicated that this 
was due to simple mixing with bypassing, and not to circulation. The results of the 
conversion vs. rpm test determined that ideal mixing and circulation could be ensured 
only at impeller speeds of greater than 3000 rpm. 

Attempts to obtain steady state kinetic data were frustrated because the catalyst 
activity appeared to increase over the course of even one run. This was found to be due 
to carbon dust originating from the lower high-temperature graphite bearing, entering 
the reactor and providing additional catalytic material. The stainless steel shield whose 
use was recommended by Turner, had been installed but the wear on the bearing was too 
great and the dust overloaded the shield. It would have been desirable to reduce the 
impeller speed —— by increasing the impeller diameter, for example —— but the reactor 
oven was too small and this was not possible. Ultimately, the frequency of maintenance 
that was required to clean the reactor vessel and also replace the carbon bearing made it 
desirable to use another reactor. 

An Autoclave Berty type reactor was obtained to replace this reactor. This 
reactor was similar to that developed by Berty et al. (6), and used by Mahoney (55) in 


studies of n-heptane isomerisation. It was further evaluated by Berty (5). This is the 
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reactor shown in Figure 3.2 and its specifications are : 


Volume = 664cm? 

Catalyst Volume = 25.7 cm: 

Max RPM = 1500 

Max Pressure = 2000 psi 

Max Temperature = 2500 K 
Catalyst Bed Diameter = 2.54 cm 


Impeller Diameter = 10.16 cm 


Kuchcinski and Squires (46) have developed an empirical set of equations for use 
in determining whether or not the Berty reactor will function properly (i.e. without 
significant temperature and concentration gradients), during low pressure operation. 
Berty (5) had suggested that, because at low pressures and corresponding low gas 
densities it is difficult to maintain sufficiently high mass velocities in the reactor, in order 
to eliminate temperature and concentration gradients, a lower pressure limit of 45 psig 
(3.1 * 105 Pa) should be observed. The mass velocity of the gas through the catalyst bed 
governs the thickness of the stagnant boundary layer at the catalyst surface and in order 
to maintain gradientless conditions , high mass velocities are necessary. Kuchcinski (46) 
had established that the Berty reactor could operate satisfactorily at 1 atmosphere by 
examining the reactor performance in the reaction of CO with O,. 

The equations that Kuchcinski and Squires derived were based on the following 
reasoning : 

If an adiabatic catalyst bed is assumed and the conditions are at steady state, the 


heat generated by the reaction is equal to the heat removed by the gas phase, 


nAQH=mCpAT 


Using the Colburn / factors (90), the mass and heat transfer coefficients can be 


estimated and the heat and mass balances yield. 
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Figure 3.2 Autoclave Berty Reactor 
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From these, it is possible to calculate AX and AT directly. 

Using the data obtained from the oxidation of CO reactor study, empirical 
equations were derived to estimate the change in temperature and conversion across the 
catalyst bed for any set of reaction conditions. 


The Reynolds number is defined as being, 


V3 
Re = (79.9 * 10-8) (M P/T) (rpom/1500\4 Dp/ L) /Nayu) 


ore (3.1) 
for Re < 50, 
AT=112rp AH Pros (T/M P) ° (1500/rpm)*5 (4/(L/Dp)?5 / 
(CD) ay iS 612, |= |e ee eee Sez) 
AX = 47000 rp » (T/P)}8/ 
[aw )+5 (Dimm)9-67 (M/ & )946 (rpm /1500)5 (4 Dp/L)-25] 
eee (3.3) 


The full calculation and the definitions of the variables are given in the Appendix 2 and 
the nomenclature. For this work, the Reynolds number in the reactor is about 20. This 
results in a conversion of 0.05% across the bed and a temperature difference of 0.02 K. 
Thus these calculations support the assumption in this work of gradientless operation on 


the Berty reactor. 


3.7.1 Autoclave Berty Reactor Description 

The reactor is shown in Figure 3.2. The entire reactor portion was constructed 
of 316 stainless steel including the gasket which sealed the reactor vessel. The impeller 
was driven by a Magnedrive drive shaft which operates by a magnetic coupling, thus 
sealing the reactor from the main drive. 

Heat to the reactor was supplied by three separate heating elements. These were 


fitted inside a heating jacket which slipped over the reactor housing. The two main 


\ wget as) “rep OOEy sabarionn dear Str «7S 
sit) (tii spiel} eran? | a 


i. 7 7 
a | 


sont «Gq nQGOTE = RL 7 


patie aig0et nay Sail anil 


2in7 Oo weds «| so9se” ah 290 


WROD Faso wih exilegmat se 
ee roliINgD cenrsivery 10 HQW atl AmB 


Ah Reon. Aigw salt we exabinamnn ent 
23 8 AO% tc nats sinavnae es in, 


N 
CO 


heaters, each rated at 2000 watts, covered the outer wall while the third, rated at 1000 
watts, covered the top. The two main heaters were linked in series and supplied by a 
208 volt supply. The top heater was connected to a 110 volt supply. The main heaters 
were controlled using a Foxboro controller, whereas the top heater was controlled using 
a Variac. 

The design had several safety features. A high temperature cut-out was used to 
protect the heaters from overheating. This was set to 500 °C. The Magnedrive shaft 
was protected from the hot reactor by a water jacket. The water flow rate was 
monitored by a conductivity switch which would cut the heater power should the flow 
cease. The reactor temperature was monitored with Chrome!—Alumel thermocouples. 
Lastly, iron—constantan thermocouples were inserted between the oven casing and the 


reactor wall to act as the sensors for the controller. 


3.7.1.1 Reactor Operation 

‘Ie The reactor screws were tightened to 175 ft-lb, using a torque wrench applied 
in increments of 25 ft-lb cross-—rotationally. 

2 The reaction temperature was reached, starting from 4 cold reactor, after about 
six hours. The reactor was kept hot over a period of several runs. 

By All runs were carried out at 1200 rpm and the tachometer was found to be 
accurate. This instrument was checked using a stroboscope. 

4. Cooling to room temperature took about three to four hours. 

C100 anti-galling lubricant was applied to the screws at each dismantling. This was 
carefully cleaned off each time and no problems were encountered with seized screws. 
Throughout this work, the reactor performed without any trouble, except 

towards the end when it was discovered that sulphur had back diffused past the top 
bearing and condensed on the shaft. The added friction caused the magnetic coupling to 
slip. The amount of sulphur was very small, however, in comparison to the length of time 
the reactor had been in operation and it was extremely unlikely that it had an effect on 


the experimental data. 
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3.8 Sulphur Condenser 

The sulphur condenser used in this work was that designed by Liu (52), with some 
modification to adapt it to high sulphur flows. The main addition was a sulphur reservoir. 
This collected all the sulphur that condensed and had sufficient capacity to require only 
infrequent emptying. This vessel was added because it was found that liquid sulphur 
would plug the line leading to the top of the apparatus. This would cause pressure 
fluctuations eventually leading to excessive entrainment and plugging of the downstream 


system. 


3.8.1 Operation 

The sulphur condenser’s performance was satisfactory considering that it was 
found impossible to remove all the sulphur in one condenser even at 115 °C (just above 
the melting point of 113 °C). It is likely that this was due to the formation of sulphur 
mist, a problem in industrial plants (9 1). 

Approximately every two to three sulphur syringe fillings, the reservoir was 
removed, and emptied with the aid of a propane torch. It was found that vibration from 
the circulation fan motor tended to loosen some of the Swagelok fittings over a period 
of time and it was important to ensure that these were tight at the beginning of a series 


of runs in order to eliminate leaks. 


3.9 Secondary Condensers 

A small amount of sulphur was found to entrain from the main sulphur condenser. 
Therefore, it was necessary to make provision for this as it would rapidly plug the line at 
the first cold spot encountered. A system of three secondary condensers (15,16, 17) 
was used to cope with this problem. This was operated in the following manner (the 


vessel numbers refer to Figure 3.1) : 


(f- Vessel 16 was unheated. It was packed with glass wool in the top section and 
with 5mm glass beads in the bottom. This effectively condensed the remaining 
sulphur and some of the water vapour, too. 


2. Vessel 17 was required to condense the remaining water, otherwise it would 
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condense in the line to the control valve. Both vessels had valves so that the 
condensed water could be drawn off rapidly. 

3 Vessel 15, in the GC sample line, was filled with 316 stainless steel turnings and 
served to remove the entrained sulphur in this line. This vessel! and the line 
downstream were heated to approximately 80 °C to ensure that no water 


condensed before it reached the GC. 


The line from the sulphur condenser to vessels 15 to 17 was maintained at a 


temperature similar to that of the sulphur condenser. 


3.9.1 Operation 

The main condensing system worked satisfactorily and the load on the secondary 
condensers was rarely too high. When it was, however, the entrance to vessel 15 
became plugged. Usually, this was found to be due to the fact that the temperature of 
the main condenser had risen. The vessel flange could be removed and the plug 


eliminated by melting It, using a propane torch. 


3.10 Gas Chromatograph 

A Hewlett Packard 5700A gas chromatograph was used in conjunction with a 
Hewlett Packard 3700A integrator to perform the product analyses. 

The product stream with the sulphur component removed was sampled manually 
for 12 seconds. The sample was passed through a five foot Chromosorb 108 
1/8th inch column, using Helium at a flow rate of 35 ml/min. The column oven 
temperature was 90 °C. A thermal conductivity detector (TCD) analysed the sample 


composition. The temperature of the TCD was 115 °C and the current 160 mA. 
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Figure 3.3 Sample Chromatogram. 


Column Chromosorb 108 1/8 in. by 5 ft. 
Oven Temperature : 90°C, 

T.C.D. Temperature: 115 °C, 

Detector current: 160 mA, 

Carrier flow rate : 35 ml/min (He). 


Chart Speed : 0.5 cm/min 


Chart Attenuation : Logarithmic, 4 decades 


TOP 
AREA & 
TYPE AREA 
M 1556539 92.68 
T 4725 . 261 3 
3239 5 Abie ky, 


75 0109 6.848 
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The GC. was calibrated several times throughout the course of the experimental 
programme in order to ensure that the response factors remained constant. The 
overheating of the column oven required the replacement of the columns and the TCD 
filaments. The new columns and operating conditions altered the response factors 
slightly. Otherwise, the response factors remained constant. The calibration procedure 
and results are reported in Appendix 1. 

A sample chromatogram is shown below in Figure 3.3. Sample analysis, at a 
sensitivity of 0.03 mV/min, and allowing for the slight tailing of the water peak, took 
between five and eight minutes depending on the water concentration. On the average, 
ten samplings per hour could be made. 

Previous studies performed by Liu (53) and Cho (15) had determined that a column 
consisting of Chromosorb 104 (4ft) and Poropak Q-S would satisfactorily separate the 
Claus reaction components. A column of this type was used initially, however, 
Chromosorb 108 gave a better separation and was adopted instead. The tailing of the 
water peak was less and the elution time shorter. This advantage may have been due to 


the low concentrations of H,S and SO, involved. 


3.11 Process Measurements 


3.11.1 Pressure 
jg Atmospheric : measured by a Fischer mercury barometer in the neighbouring 


laboratory, reading in mmHg. 


oe Nitrogen : measured using a USG gauge, calibrated at 8.5 psig, reading accurate to 
0.1 psig. 
So Reactor : measured just upstream of the sulphur mixing, using a Dynisco 


pressure transducer PT25—30 and a Foxboro chart 
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3.11.2 Temperature 
1 Feed : sensed by an Iron—Constantan thermocouple. 


oe Reactor : sensed by a Chromel—Alumel thermocouple. 
All thermocouples were connected to a Honeywell 24 point chart recorder 
calibrated to read O-30 mV. The reference temperature was 0 °C, using an ice—water 


bath. 


3.11.3 Flow Rates 


iF Nitrogen : Volumetric flow rate measured by a differential pressure cell. 
Zi Water : Mass flow rate obtained from the syringe pump calibration. 
3: Sulphur : Mass flow rate obtained from the syringe pump calibration. 


3.11.4 Process Control 


iF Nitrogen flow : Moore differential flow controller. 
2. Reactor pressure : Foxboro controller. 
3 Reactor temperature: Foxboro controller. 


3.12 Calibrations 
tt Pressure Transducer 
em ro =10189557 28h G101379 


where R ts the chart reading. 


Pe D/P Cell 
Q? (I/min)? = 0.089984 * R% — 0.053501 
3: Thermocouples 
Feed : Accurate to within recorder reading, +0.15 mV. 
Reactor : Accurate to within recorder reading, +0.15 mV. 
4. Water Syringe Pump 


Range 0.01:— G (g/min) = 4.66522 * 10% * R% — 9.3262 * 10° 
Range 0.001:- G (g/min) = 4.56200 * 10° * R% — 1.0665 * 10° 
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ee Sulphur Syringe Pump 
Gear g-atoms S,/min 
L1 = 4.686 * 10? 
H 1 = POLGU Ox 
aes = 9S79Ga0" 
H2 = 1.4817 * 107 
LS = 1.9420 * 107 
H3 - 3.0729 * 9102 

6. G.C. Response Factors 


2 
H,O = 0.8747 + 0.01918 * mol% — 0.00091 * (mol%) 
HS = 1.841969 
SO7=91:2077 


3.13 Materials 
Specifications supplied by the manufacturers. 

Nitrogen : Linde Gas products High Purity 99.99% .......... O, main impurity 
Sulphtr Fisherm ni orade sublimatec |. fs. aaa eee H,O0 main impurity 
VValtenmenniplecisuiGdi24):-99, 99700. in see ee eee N, main impurity 
Catalyst : Kaiser S201 alumina 

area 380 m?’/g 
chemical composition, 
SiO, 0.02% 
Fe,Q, 0.02% 
TiO, 0.002% 
Na,O, 0.3% 
Al,O, 93.6% 
loss on ignition 6.0% 


size, 8-39 mesh, 2.00-2.38 mm 
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3.14 Alternative Equipment 

Over the course of this study, several forms of equipment were either modified 
or discarded. The reactor used in the preliminary stages of this work and the reasons 
for discarding it have been discussed. This was replaced mainly because it required 
frequent maintenance. The original gas chromatograph, a Microtek 2003, was replaced 
with a Hewlett Packard 5700 for the same reason. The recorder baseline often became 
noisy. This problem was attributed to an indetectable air leak and was not encountered 
with the Hewlett Packard gas chromatograph. 

The main disadvantage of the experimental system was the irregularity that was 
observed in the water feed when both the sulphur and water pumps were operating. 


Two approaches could be taken in the attempt to solve this problem. 


Ue Sulphur Feeder 

Since the syringe pump itself appeared to operate smoothly, the root of 
the problem was likely to be the successive drops of sulphur evaporating 
irregularly in the sulphur mixing vessel. The surge vessel was inserted to smooth 
the flow. This vessel's capacity was large enough to accommodate the flow 
between one sulphur drop and the next. No difference in the flow characteristics 
was observed between high and low flows of sulphur and so it was unlikely that 
an even larger vessel would have eased the problem. 

The replacement of the pump itself with a sulphur saturation system was 
considered. Many operational problems were anticipated, however, and the use 
of such a system would have required the entire rebuilding of the reactor feed 
system. It was decided to continue with the syringe pump in order to determine 


whether sufficiently precise data could be obtained with the original equipment. 


2 Water Feed 
It was possible that the feed fluctuation was due to the actual variation in 
the water flow itself. It would have been possible to eliminate this using a water 
saturator. Since this equipment would have been considerably easier to build and 


test than the sulphur saturator, preliminary tests were made employing a water 
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saturator without the equipment required for the fine metering of the flow. This 
apparatus consisted of two water spargers in series, immersed in a 
thermostatically controlled water bath. 

Unfortunately, several problems were encountered. The system had a high 
pressure drop which restricted the flow rates that could be obtained with the 
normal pressure available from a standard nitrogen cylinder (60 to 70 psig). The 
fixed sulphur flow rates meant that the water flow had to be easily adjustable. 
This was especially desirable in view of the amount of data required in an 
experimental initial rate study. The main reason, however, that the saturation 
system was discarded was the fact that the scale deposited in the sulphur mixing 
vessels proved to be catalytic, and it was necessary to add the water feed just 
before the feed stream entered the reactor. In order to allow a sulphur saturation 
to operate, two streams, one of water/nitrogen, and the other of 
sulphur/nitrogen, would have to be blended. By careful control of the 
compositions of each stream and the ratio of their flow rates, the required feed 
stream composition could be achieved. A much more complex feed system 
would therefore have been required in order to improve on the one used in this 


work. 


To summarise, it is likely that the above problems could be solved, but that a very 
much more complicated and expensive feed system would be necessary. This would 
probably involve the use of water and sulphur saturators, and include the use of two 
separate feed streams, S/N, and H,O/N,. Each of these streams would require on-line 
analysis and precise flow measurement and control. A sulphur saturation system would 
require precise control of high temperatures, 300 —- 400 °C, and of nitrogen flow rates, 


in order that the particular feed flows required could be quickly and reliably obtained. 
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4. EXPERIMENTAL PROGRAMME 


4.1 Introduction 

The Claus reaction is a reversible reaction with the equilibrium favouring the 
reverse reaction, that is, the production of sulphur and water. With such a reaction, it is 
necessary to consider both the forward and the reverse components in any kinetic 
analysis. This requirement necessarily doubles the number of parameters to be 
determined by the data reduction. It is advantageous, therefore, to simplify the rate 
equation by considering the special case of the initial rates. This will set the product term 
in the overall rate equation at zero, in which case the reaction may be regarded as a 


single, irreversible one. For example, for the general gas phase reaction , 


A+B = C+D 


(4. 1) 
the general rate equation is of the form 
—ha—'K far kare rd 
(4.2) 
In the initial rate case, the last term —~> O, 
—Riinitial rate) = k, Pa Pb 
(4.3) 


The number of experiments required to collect initial rate data is far greater than 
that required for irreversible reaction studies, since a series of finite rate data points are 
required to obtain one initial rate datum. The initial rate method, however, reduces 
considerably the total number of parameters that must be estimated from the data This 
allows the data to be used to obtain more specific information on the reaction rate. Also, 
the data is not wasted since, once the forward and reverse terms have been determined, 
all the data Botlectad may be used to estimate the form of the overall rate equation. The 
present study does not include such an analysis of all the data. 

McGregor (56) used the initial rate method to investigate the forward Claus 
reaction. He followed the data analysis suggested by Mezaki and Kittrell (59). This work 


will adopt that same approach. Essentially, the following argument is applied (48) : 
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The mass balance for a CSTR yields ; 


-R = Fa Xa/Wt 


Rearranging this we obtain ; 


Xa = —R Wt/Fa 


A plot of Xa versus Wt/Fa will yield a smooth curve, the slope of which, at a 
given conversion and reciprocal mole flow rate, is equal to the rate of reaction. 

However, this rate of reaction is not the overall rate of reaction since the molar feed rate 
to the catalyst bed is not the same as the molar feed rate to the reactor. In the limit, as 
Wt/Fa -—> 0, Fa —-> infinity , the mole feed to the bed and the reactor are almost 
identical. As aresult, the slope of this curve at the origin, that is, zero conversion, will be 
equal to the initial rate of reaction. 

The slope at the origin may be derived from the experimental data by fitting the 
points with some smooth curve that passes through the origin by using a least squares 
parameter estimating technique. Mezaki and Kittrell (569) suggested using the nonlinear 
curve, 


Y = A tanh(B X) 


where X is the independent variable and Y the dependent variable. A and B are 
parameters which must be determined by a least squares analysis. This curve has the 
advantage that its derivative at the origin, which is equal to ( A*xB ) , may be calculated 
analytically (Appendix 3). Thus, the initial rates may be obtained quickly from a non-linear 
least squares fit of this function to the data) A series of conversion vs. 
(mole flow rate) -! curves will therefore yield the initial rate vs. partial pressure data 
required to perform an empirical study of the reverse reaction. 

Once the orders of the reaction have been ascertained, the rate constant of the 
reaction can be calculated. By conducting initial rate runs at different temperatures, it is 
possible to determine the apparent activation energy and the temperature dependence of 


the rate constant by a plot of log(k) vs. 1/T (48). 
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4.2 Preliminary Experiments 


Several preliminary experiments were necessary to ensure that the rate data 


obtained were reliable. 


4.2.1 System Inertness 

Initially, and at several times throughout the experimental study, the reactor 
equipment was tested for catalytic activity. This was done simply by feeding the 
reactants into an empty reactor and testing for reaction products. The spinning basket 
reactor system originally used (82) was found to be catalysing the reaction. By careful 
cleaning of the equipment, the reaction products, attributable to residual carbon dust in 
the reactor, were reduced to a negligible amount. Though the Autoclave reactor was 
found to be unreactive at the beginning of its use, after a period of time reaction was 
discovered to be taking place in the sulphur mixing section. This was attributed to a 
carbonaceous deposit which was found to have coated the inside surface of the mixing 
vessel. It is believed that this materia! originated in the sulphur syringe itself, possibly 
from the teflon "O’-rings, and possibly from the vaive packings. The valves themselves 
had to be replaced on two occasions because they leaked. Complete cleaning of the 
interior of the sulphur mixing vessels proved to be impossible, so they were replaced. 
The problem was effectively eliminated by placing the water feed just before the 
reactor. The deposit was effectively trapped by the mixing vessels and there was no 
further problem with extraneous reaction. Unfortunately, however, this arrangement 
precluded the use of water saturators for the water feed. 

Particular attention was paid to the possibility of reaction within the sulphur 
condenser. Karren (39) had reported that liquid sulphur catalysed the reaction. The tests 
for system inertness were carried out with the full experimental system which included 
the sulphur condenser. Since liquid sulphur was present in the condenser and no 
hydrogen sulphide or sulphur dioxide were detected, it was concluded that liquid sulphur 


did not catalyse the reaction. 
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4.2.2 |ldeal Mixing and Bulk Diffusion 

All reactors of the CSTR configuration must be evaluated to ensure that the basic 
assumption of ideal mixing is valid. The Berty reactor can easily be tested by running the 
reactor at constant feed composition and varying the impeller speed. Above a certain 
impeller speed, the conversion will remain constant, indicating that the mass velocity of 
the reactants across the catalyst bed is no longer influencing the reaction rate. Figure 4.1 
shows a plot of conversion versus impeller speed. This shows a regime in which the 
impeller speed has no effect on the conversion. As the figure indicates, ideal mixing is 
achieved at impeller speeds of greater than 1000 rpm. 

The feed conditions used in this test were more severe than those used in the 
kinetic runs. A higher space velocity and catalyst loading ensured that the actual kinetic 
experiments would be carried out under ideal mixing conditions. An impeller speed of 
1200 rpm was therefore chosen for all runs. 

The impeller speed experiment also determined that neither interphase, nor bulk 
diffusion effects, were rate limiting. 

(A similar experiment with the Turner spinning basket reactor determined that 
much higher speeds were required to ensure ideal conditions pertained inside the reactor 


vessel because of the small diameter of the impeller.) 


4.2.3 Catalyst Activity 

It is essential in any heterogeneous catalytic study that the catalyst activity remain 
constant over the period of the experiments, otherwise it becomes an additional variable, 
one which is difficult to quantify. The analysis of the deactivation of catalysts Is a 
specialised field of study. In this study, Kaiser S201 alumina catalyst was used. This ts 
reported to be a very robust catalyst (67) which maintains its activity over long periods 
of time. 

lt was found that the apparent activity of the catalyst in the Turner spinning basket 
reactor system increased with time. This was attributed to carbon dust rather than to any 
real increase in activity. 

With the Autoclave reactor, the activity of the catalyst was found to remain 


constant over the period of the experimentation. 
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Xw, Fractional Conversion of Water 
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Figure 4.1 Conversion versus Impeller Speed 
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4.2.4 Volume of the Reactor 

The volume of the reactor was measured in the following way. For an ideal gas at 
constant temperature, P V = constant. From this relation, the volume of a vessel may be 
calculated. The theoretical analysis is given in Appendix 1. 

Figure A- 1.2 in the appendix illustrates the apparatus used to determine the 
reactor volume. A cylinder of known volume (500 in3) was filled to about 70 psig and 
connected via a valve to the reactor. Atmospheric pressure was recorded and the valve 
was opened siowly to allow the gas to expand into the reactor ideally. This procedure 
was repeated six times and the volume of the reactor was calculated each time. The 


actual reactor volume was taken as the average. 


4.2.5 Repeatability of G.C. Peak Areas 

The accuracy of any reaction study is limited by the reliability of the product 
sampling techniques. It was found that the water syringe feed gave a constant water 
flow rate as shown by the GC. water peak area. For example, for one particular flow 


rate, 


Average H,O Peak Area % = 1.920 + 0.0590 (3.07%) 


The syringe was calibrated for the entire range of flows to be encountered, and the 
actual water composition in the stream was known. It was this method that was used to 
calibrate the G.C. water peak. The introduction of the sulphur feed, however, caused the 
water peak to fluctuate by as much as 10%. As a consequence, it was decided to 
calculate the water feed concentration directly from the products using the hydrogen 
sulphide peak in order to establish the actual water feed. 

From the product analysis, the molar flows of H,S and H,O were calculated. From 
the reaction mass balance, the molar flow rate of water in the feed stream must be equal 
to the sum of the H,S and H,O product flows. 

This approach was necessary since the introduction of the sulphur feed affected 
the water feed rate and as a result the actual water feed flow was unknown. The 


material balance across the reactor can now be checked by the g-atom balances on 
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oxygen and sulphur only. 

An important consequence of using this method is that the observable error 
associated with the product analysis is distributed between the calculation of the water 
feed and the product, rather than being concentrated in the product analysis. This will be 


discussed more fully in the discussion section 5.8. 


4.2.6 Ideal Gas Law 

In the material balance calculations used in deriving the rate data in this study, the 
ideal gas law has been assumed to be valid. This assumption is shown to be valid by 
calculation of the component fugacities, (Lewis (50)). Pure component fugacity 
coefficients may be estimated using generalised compressibility charts and the principle 
of corresponding states. The critical constants and the reaction system are shown 
below. The fugacity coefficients are estimated from the generalised property charts of 
Newton (64). It is apparent that, at the temperature and pressure of the reactor, the 
fugacity coefficients are very close to unity. Though the charts are not precise at such 
low reduced pressures, the error is but a few percent. At room temperatures, 
McGregor (56) indicates the gases may all be treated as ideal. Sulphur is not part of the 


material balance at room temperature. 


Table 4.1 Critical Constants and Fugacity Coefficients. 

Species TeX) Pc (atm) Ze Pr ite (GO OK mae 
(1.46atm) 

N, 126.0 335 0.291 0.0436 4.76 1.0 

me S756 88.9 0.284 0.0164 1.606 1.0 

SO, 430.7 Uke: 0.269 0.0188 Shee 1a) 

H,O 647.4 2.5.3 0.230 0.0067 0.927 1.0 

S) USE SAI ASL O02 0.0081 0.457 0.98 


The critical property data were obtained from JANAF (36) thermochemical tables 


for all species except sulphur whose properties were obtained from Rau et al.(7 1). 
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4.2.7 Sulphur Vapour Equilibrium 

It has been assumed in this study that the sulphur vapour has sufficient time to 
reach an equilibrium distribution of species. Although data is available on the equilibria of 
the various sulphur species, none is available on the rates of these reactions. Berkowitz 
(4) considers these rates briefly. He reports that there is evidence that the rate of 
equilibration of the sulphur species is slow. He suggests that the activation energy of the 
equilibration is probably related to the strength of the sulphur—sulphur bond in the ring. 

A large surface area is reported as being an aid in achieving equilibrium sulphur vapour. 
The nature of that surface was also a factor. Qualitative experiments to determine the 
rate of equilibration were made over the range of 300-350 K. 

The temperature and design of the sulphur mixing system were such that they 
would tend to improve the rate of equilibration. The sulphur mixing vessel was tilted at 
an angle of approximately 25° from the vertical. This would ensure that the liquid sulphur 
droplets would fall directly onto the vessel wall providing a large surface area for 
evaporation. The sulphur inlet line was 1/16 inch stainless steel tubing which projected 1 
inch into the vessel and which had been filed to resemble a syringe needle. This 
arrangement ensured that the incoming sulphur contacted the nitrogen flow directly and 
the needle itself provided as large an area as possible for evaporation. A surge tank was 
also added to increase the time available for the mixture to equilibrate. There is no 
simple test, however, to determine whether the assumption that the sulphur vapour 


equilibrium is attained is a reasonable one. 


4.2.8 Intra-phase Pore Diffusion 

Diffusional resistances within the catalyst pellet may be rate-limiting (10, 48, 73). 
These effects depend on the relative rates of diffusion and on the surface reaction. The 
single pore model, as analysed by Levenspiel(48), shows that the rate of reaction will be 
diffusion-limited if (L k/D) < 0.5. Thus, for small L, pore length (that is, small catalyst 
pellet), slow reaction rate, and high diffusion, it is less likely that pore diffusion will be 
rate-limiting. 

It is possible, then, to test whether such an effect is present by a series of 


experiments measuring the conversion versus catalyst pellet size. McGregor (56) 
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performed this test with two catalyst sizes, 0.6 mm and 0.25 mm. The rate of reaction 
was much higher with the smaller catalyst pellet. McGregor concluded that the rate was 
proportional to the external surface area of the catalyst. George (30) and Landau et al. 
(47) determined that, with pellet diameters of over 2 mm, all the rate function 
parameters, except the rate constant, remained constant. Liu (53) analysed his ir. studies 
of the Claus reaction using criteria suggested by Hudgin (89) and also by estimating the 
effectiveness factor. He used catalyst wafers 2 mm thick. He concluded that pore 
diffusion was insignificant in his studies. 

No test for pore diffusion effects was made in this study. McGregor’s 
conclusion that the reaction rate was proportional to the external area of the catalyst is 
assumed to be valid in this study. The catalyst size used in McGregor'’s work was 0.6 mm 
as compared to 2 mm in this study. The reaction rate, however, is much less —- by a 
factor of ten , (see Appendix 2) — so that the possibility of rate-limiting pore diffusion is 
much less. The possible effect of significant pore diffusion resistance in the light of the 


data obtained in this work will be examined in the discussion section, 5.8. 


4.3 Kinetic Experiments 

The reverse Claus reaction is a function of the partial pressures of water and 
sulphur, and of temperature. Since it was expected from the results of previous studies 
on the forward reaction (30,53,56 ) that water adsorption would inhibit the reaction, the 
easiest parameter to isolate and evaluate was that of the order of the reaction with 
respect to sulphur, since in that case the water partial pressure is kept constant as is the 
reaction temperature. Having established the sulphur order of reaction, the dependence 
of the reaction on the partial pressure of water can be investigated with a view to 
establishing whether or not water does indeed inhibit the reaction, and whether an 
adsorption term in the rate expression can satisfactorily explain this inhibition. Additional 
runs at higher temperatures may then be made in order to investigate the temperature 


dependence of the reaction. 
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4.4 Experimental Procedure 


4.4.1 Start-up 

Initially, the entire system was flushed with nitrogen. The temperatures of the 
various sections were adjusted when necessary. The entire system including the preheat 
oven, sulphur condenser, et cetera, took about ten hours, when starting from room 
temperature, to come up to operating temperature. About one hour before running, the 
flow rate of nitrogen was set to that required for the experimental runs in order to allow 
the temperatures to equilibrate. 

The flow rates of N, and H,O were determined by the particular sulphur syringe 
gear that was engaged. The settings that were required to obtain these flow rates were 
calculated from the calibration equations for the D/P cell and the water syringe pump 
respectively. 

Once the temperatures and pressures were as required, the water feed and 
reactor fan were started. It was important to check at this stage that the system was 
functioning as expected, since it was possible that the water feed line had become 
blocked. The gas chromatograph was used to check that the required partial pressure of 
water was being achieved. 

Once the system had come to reasonable steady state, as determined by the 
temperatures, flow rates and GC. analysis, the sulphur flow was started. The appropriate 
gear having been engaged, the sulphur feed valve was opened. Unfortunately, the only 
indication of sulphur flow was the appearance of reaction products in the effluent 
stream. When starting from a sulphur-free system, it would be one to two hours before 
the sulphur feed was fully established and steady state attained. 

The required feed rates having been attained, quantitative analysis of the product 
stream was possible. The reaction products were monitored over a period of about an 
hour in order to establish the true steady state reaction. The general procedure was to 


obtain six analyses at the steady state which might then be averaged. 
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4.4.2 Running 
Once the system was operational, the steady state could be altered quite rapidly 
by adjusting the feed rates of the reactants. The reactor returned to steady state within 
15 minutes of its being changed. In this manner, several runs could be made in one day. 
It was important that the temperatures and pressures of the system be monitored 
regularly to ensure that any problems, especially plugging by sulphur, were detected 


quickly. 


4.4.3 Shutdown 

After a series of runs were completed, the system was shut down. It was 
important that this be done carefully to avoid sulphur condensing anywhere other than 
inside the sulphur condensers. This was achieved by immediately shutting off the sulphur 
feed and closing the sulphur feed valve, thus isolating the syringe from the main system. 
The relative amount of sulphur passing through the reactor could be estimated by a 
continuation of the monitoring of the product stream with the GC. On average, it took 
about 60 minutes for the amount of reaction products to become negliglible. At this 
point, the water feed was stopped. Full nitrogen flow was maintained for another 30 to 
60 minutes and then was reduced to a few cc's per minute. The drain valves of the final 
condensers were opened at this point to drain off the collected water and also to ensure 
that the low flow nitrogen had free passage. 

The entire system was left at the operating temperature overnight so that the next 
set of runs could be begun quickly. This was also found to be extremely useful as an 


extra precaution against the condensation of any residual sulphur in the lines. 


4.4.4 Maintenance 

The principal maintenance of the reactor system involved the periodic cleaning of 
the sulphur condenser. The condenser was not arranged so that this could be performed 
on-line. Though on-line cleaning would have been an advantage, the modification 
described in the equipment section considerably increased the period between drainings, 


and therefore decreased the need for cleaning. 
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The build up of scale inside the sulphur mixing vessels was the only other serious 


problem. This was bypassed by placing the water feed just before the reactor. 


4.5 Data Reduction 

Most of the experimental measurements of temperatures, pressures, and flow 
rates, were obtained in terms of a percentage of full scale readings on an electronic 
recorder. The rest were obtained by a variety of other means. Atmospheric pressure, 
for example, was read from a barometer in mm Hg. Feed water flow rates were 
obtained from the calibration equation of the syringe pump (see Appendix 1) and also 
from the reaction mass balance. Nitrogen feed pressure was read in psig directly. The 
reactor and feed temperatures were converted from the mV readings obtained from the 
temperature recorder, to degrees Celsius using the thermocouple tables from the CRC 
handbook (14) and the calibration factor. The pressure data were converted to absolute 
pressures in atmospheres according to the calibration equations as derived in Appendix 1 
The nitrogen flow rate was converted to litres per minute by means of the calibration 
equation as shown in Appendix 1. 

The product compositions were calculated from the gas chromatograph peak 
areas for each component as shown in Appendix 2. These gave a mole per cent on a 
sulphur free basis. The actual mole percentages were obtained using the mass balances. 
The ideal gas law was applied to all species and the nitrogen gas stream served as a basis 
for calculation since this component was inert and the same amount comprised the feed 
and the product streams. From the calculated nitrogen flow and the product 
compositions, the molar flow rate of each of the products was obtained. 

The equilibrium distribution of the sulphur species was estimated by using the 
Gibbs free energy minimisation method and by assuming that only the sulphur species 
S,, Sy, S,, and S, were present in the vapour. It is known that this is only an approximate 
method of estimating the sulphur species distribution (72) but these species were chosen 
with the view that the free energy data on them are the most reliable available . 

The fractional conversion was obtained from the difference in the molar flow 


rates of each of the reactants. 
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The partial pressures of all species were calculated using the calculated sulphur 
species distribution, and assuming these species to be in equilibrium with each other, the 
mole fractions of the other species, and the total pressure. 

A computer program for carrying out these calculations, including the processing 
of the raw data and the reactor material balance, is shown in Appendix 3. A listing of this 
program, a flow chart describing its operation, and a detailed calculation of one set of 


raw data is given. 


4.5.1 Initial Rate Calculation 

In order to obtain initial rate data, it was necessary to plot conversion versus the 
inverse molar flow rate and obtain the slope of the resulting curve at the origin. This was 
done by fitting the curve Y =A tanh (B X) by means of a non-linear curve fitting 
technique. The method used was that of Meyer and Roth (58) and it proved extremely 
effective in the calculation of a least squares estimate of the parameters of the Tanh 
function, provided that the experimental data had been scaled to ensure that the 
exponentials were small. The initial rate is given by the derivative of the function at the 
origin. This function may be differentiated analytically as shown in the Appendix. Thus, 
the initial rate is obtained by multiplying the non-linear least squares estimates of the 


parameters A and B. Appendix 3 lists the program used to fit this data. 
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5. RESULTS 


5.1 Introduction 

A total of 173 finite rate experiments were performed. Of these, 113 were 
performed at a constant temperature and water partial pressure. From this data, the 
initial rate Of reaction is derived for 10 sulphur partial pressures. In addition, the 
equilibrium conversion is estimated for 6 different sulphur partial pressures. Of the 
remaining runs, 50 were carried out at constant sulphur partial pressure and a further 8 
initial rates were obtained. An additional 10 runs were performed at two higher 
temperatures, yielding another two initia! rates. 


Each individual run is identified by a five character name, viz. 


SIO4A, - Sulphur Initial rate data set, 
ie. sulphur partial pressure varying. 
04, 10 * nominal total sulphur percentage in feed (0.4%). 


A, individual run identifier 


hence, 
WI01A —- Water Initial rate data set, 
ie. water partial pressure varying. 
01, nominal percentage of water in the feed (1%). 
A, individual run identifier. 
and, 


WT6HA — Water Temperature run, 
ie. initial rate with high water partial pressure and higher 
temperatures for the series of runs. 
6, nominal percentage of water in the feed, (6%) 


A, individual run identifier. 


A total of 53 runs were made at low space velocities during the development of 
the reactor system. These determined the experimental conditions required to obtain 


repeatable data and also provided a series of data sets from which it is possible to 
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estimate the equilibrium conversions. 

Prior to reducing the quantity of catalyst employed in order to allow the taking of 
initial rate data, the GC. oven ran "wild”. Inexplicably, the oven heater was suddenly given 
full current, and the temperature rapidly increased, baking the column and fouling the 
thermal conductivity detector filaments. The column and the filaments were replaced and 
the GC. recalibrated. The problem recurred on two occasions, however a lead thermal 
fuse was used to protect the column from overheating and there was only a short 
interruption in obtaining the experimental data. As a result of this problem, all the 
experimental runs used to obtain the initial rate data were performed after the G.C. had 
been re-calibrated, from October 1980 on. Runs conducted prior to this date were not 
used in the initial rate analysis but were used to determine the equilibrium conversion. 
The GC. calibration for these preliminary runs is therefore slightly different. 

Appendix 3 contains both the computer programs and sample outputs that were 
used to calculate the values of the kinetic parameters of the reaction. The computer 
program EXPTCAL was used to analyse the raw data obtained from each run. A sample 
of the main output from this program is given after the program listing. A summary of 
each and every run that was performed is contained in Appendix 4. The program 
REGRSS was used to analyse the conversion/space time experimental data. The 
subroutine NLSYST established the non-linear least squares fit to the data. The output 
from this program illustrates the experimental data points and the least squares estimated 
points. This output is contained in Appendix 4. The program IRANAL was used to 
calculate the least squares estimate of the parameters of the function which it had been 
suggested would describe the water dependency of the reaction (equation 1.1). The 


output from each of the various estimates of these parameters is given in Appendix 4. 


5.1.1 Sulphur Dependence 

Initial rate experiments were performed to determine the order of the reverse 
Claus reaction with respect to sulphur. These were performed at a constant pressure, 
temperature, and partial pressure of water. The partial pressure of sulphur was varied 
over the entire range that the experimental system would permit. The mole fraction of 


sulphur(8) ranged from 0.14 - 2.0%, corresponding to a partial pressure of 
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266 — 3000 Pa The partial pressure of water was held constant at about 

2.5%, ie. 3700 Pa. The temperature was 592 K and the pressure, 1.46 atmospheres. A 
total of 55 runs were made at 10 different partial pressures of sulphur yielding 10 initial 
rates. At least 3 runs were made for each set of conditions at molar feed flow rates as 


high as possible in order to obtain data near to zero conversion. 


5.1.2 Water Dependence 

The water dependence of the reaction was analysed by performing initial rate 
experiments at a constant sulphur partial pressure of 753 Pa The temperature and 
pressure were equal! to that used in the sulphur experiments, that is 592K and 1.46 atm 
respectively. The mole fraction of water used varied from 1 to 14.5% and eight 
different partial pressures were employed. The partial pressure of S, was about 0.5% of 
the total pressure. A total of 50 runs were made over a range of inverse molar feed 
rates of 1 —- 40 g hr/mol. 

A minimum of 5 runs were made at each set of conditions. Five runs were 
possible since the water feed rate could be smoothly varied over a wide range. As with 


the Si experimental series, the highest molar feed rates possible were employed. 


5.1.3 Temperature Dependence 

Two additional initial rate data sets were carried out at temperatures above 
592 K, namely 624 K and 653 K. Higher temperatures were chosen to avoid any 
problems that might have arisen by operating close to the sulphur dew point. At the 
partial pressures of sulphur employed, the sulphur dew point was around 550 K. The 
partial pressures of the reactants were kept close to that of the 592 K set of runs, 
WIO6, although exact correspondence could not be achieved since the higher 
temperatures are expected to alter the sulphur species distribution. An additional 10 


runs were performed which, combined with the WI06 data set, give the initial rate of 


reaction at three temperatures. 
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5.2 Analysis of SI data 

A total of 62 runs were performed at a constant water partial pressure in order 
to establish the dependence of the reaction on the partial pressure of sulphur. These 
runs are summarised in Table 5.1 below. The full raw data are given in Appendix 4. The 
conversion and space time data calculated in the initial data analysis, are used to obtain 
the initial rate of reaction. Figure 5.1 illustrates this data for one data set SIO4. The 
results of the non-linear curve fitting program REGRSS are contained in Appendix 4 and 
these show both the experimental data and the fitted data) The run SI25B was not used 
in the data analysis since the calculated conversion was much lower than that suggested 
by a smooth curve through the remaining points. An inspection of the raw data reveals 
that the pressure in the experimental system was much lower than usual for an 
experimental run. This would suggest for example, that a drain valve had been left open. 
This would reduce the molar feed rate of nitrogen and water increasing the relative 
amount of sulphur and as a result this run cannot be used for initial rate analysis. 

Having obtained the initial rate data, the partial pressure of sulphur(8) may be 
plotted as a function of the rate of reaction to illustrate the relationship between the two. 
A simple plot of rate against partial pressure yields a curve. The form of the rate 
equation that is thought to pertain (equation 1.11) suggests that a log-log plot will show a 
straight line, the slope of which is the order of the reaction with respect to sulphur. The 
log-log plot (Figure 5.2) shows the data collected at various sulphur concentrations and 
indicates a straight line correspondence. This plot is used to derive the order of reaction 


with respect to sulphur by linear regression, 


order wrt S, = 0.397 + 12% 


The limits of accuracy have been determined from the slope variance as shown in 


Section 5.6. 1. 
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TABLE 5.1 DATA SETS SI (summary) 
RUN PRESSURE TEIVIP WATER 

Atm K percent 
SIO4P 1.4783 Soule to 1.8063 
S!040 1.4736 591.15 2.1248 
SIO4N 1.4783 Doi tS 2505 14 
S!04Q 1.4754 eye Mt eas Pale 
SIO04R 1.4757 SoS 2.0736 
SI04S 1.4780 So lade 2.0490 
Avgs 1.4767 Ooiee 2.0062 
SIO5H 1.4757 ays iii [8S 2.8147 
S!05| 1.4809 592715 Sul2S 
Si05u 1.4792 59265 2.4115 
SIO5K 1.47932 eke itis) 2.4417 
Avgs 1.4789 O92 2.4235 
SI06! 1.4839 So 215 217162 
SIO6K 1.4792 ope) eats) 25026 
SIO6J 1.4816 59715 2.8954 
SIO6L 1.4796 5S i515 1.9773 


Avgs 1.4813 591.4 2.5452 


SULPHUR 
percent 
Ogles2 
ON7 23 
0.1781 
0.1862 
0.1786 
0.1858 
0.1784 


0.2594 
0.2435 
O250)7 
0.2194 
0.2433 


0.2896 
0.2936 
0.3057 
0.3030 
0.2980 
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Wt/MOLES X 

g hr/mol Conversion 
16y035 0.0903 
213.74 0.1089 
630.05 0.1343 
126.98 0.0876 
80.914 0.0586 
62.844 0.0463 
ShAG9 0.1341 
204.43 0.1128 
60.70 0.05416 
Ache, 0.06101 
119.84 093s 
5ea6 0.05022 
74.81 0.07074 
38.07 0.04053 


Table 5.1 continued. 


RUN 


SIO8L 
S!I08M 
SIO8N 
S!080 
Si08Q 
SI08R 
SI08S 
SIO8T 
SI08U 
SIO8V 
SIOSW 
Avgs 


S!110K 
$!110J 
S!I10L 
SI10ON 
S!110M 
SOR 
S$!10Q 
$1100 
S!10S 
S!10R 
SI10U 
Sus 


Avgs 


PRESSURE TEMP 
Atm K 
1.4608 Se he Ns) 
1.4728 594.15 
1.4683 siskee) 
1.4681 Sooo 
1.4759 BU 1s15 
1.477 1 Do ale 
1.4749 OStaS 
1.4820 Boae1S 
1.4763 5Sd15 
1.4792 59115 
1.4821 5USals 
1.4755 5926 


1.4669 ae) S5 1h) 
1.4680 Soauo 
1.4669 594.15 
1.4751 Soaeho 
1.4681 B9B.1'5 
1.4749 ete mode: 
1.4820 See Ws 
1.4820 sie ihe 
1.4793 SRS) aS) 
1.4397 Sees. 
1.4596 594.65 
1.4584 abe) ie 
1.4680 ase) 


WATER 
percent 
2.0413 
22316 
2.0014 
2.1379 
2.1380 
Zt oS 
LBvva 
2.0055 
2.5693 
2.4744 
21059 
2.1584 


2.1451 
Cay S| 
2.1586 
2,032 
2.2994 
BEA As: 
1.9474 
230 
20900 
2.6130 
2.2145 
PLETE 9G) 
23oie 


SULPHUR 
percent 
OS939 
0.3742 
0.4072 
0.4066 
0.4127 
0.4014 
03792 
0.3809 
0.4062 
0.4145 
0.4027 
0.4013 


0.6120 
0.6048 
0.6037 
0.6141 
0.6114 
0.6287 
0.6056 
Dios a2 
0.5741 
0.6036 
05737 
0:5857 
0.6008 
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Wt/MOLES X 

g hr/mol Conversion 
469.54 0.0887 
309.96 00299 
235.43 0.0864 
235.47 0.0876 
188.48 0.0874 
147.80 0.0801 
74.96 0.0544 
96.94 0.0616 
71.88 0057s 
56-Gn 0.0494 
37.04 0.0401 
268.22 0.0708 
443.66 0.0728 
224.63 0.0696 
144.25 0.0667 
223.00 0.0762 
TOL Ve 0.05090 
45.60 0.04327 
20:29 0.06182 
Sie lele) 0.05608 
178.70 O10697, 
21.0) 0.03038 
34.76 0.04138 
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Table 5.1 continued. 


RUN 


S115G 
S!I15H 
S!151 

Sin 5d 
S!15K 
SI15M 
Shor 
SI15N 
Si50 
S!115Q 
SITS 
Avgs 


RheooURE JEME 
Aim K 
1.4669 ee Shriis) 
1.4633 Sonia 
1.4751 sie Sie 
1.477 1 SoD 
1.4783 Bgie5 
1.4820 Ss ig) 1's) 
1.4796 esi 
1.4691 he Ps, MS: 
1.4742 Bg4515 
1.4617 ash) NS 
1.4596 Soins 
1.4729 Somes 


WATER 
percent 
23038 
ZAG 
page 21S 1S 
2.3089 
2.3948 
2.4319 
2.2564 
2.2650 
3.1418 
2.4200 
2.7280 
2.387 1 


SULPHUR 
percent 
0.8776 
0.8684 
0.8456 
CEO 
0.6323 
0.8580 
0.8772 
0.9070 
0.8688 
0.9011 
0.9276 
0.9015 


Wt/MOLES X 

g hr/mol Conversion 
275.42 0.0613 
ATS26 0.0574 
T3853 0.0586 
1336 0.8576 
86.80 0.0554 
43.65 0.0440 
67.30 0.0472 
50.94 0.0379 
33.68 0.0455 
16.48 0.0247 
25:39 C0312 
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Table 5.1 continued. 


RUN 


SI2Z0A 
SIZ0e 
S!I20B 
Avgs 


SI25A 
SIZoB 
SI25C 
Avgs 


SIB0A 
SI30B 
SI30C 
Avgs 


SI35A 
SI35B 
See 
Avgs 


PRESSURE TEMP 


Atm 
1.4648 
1.4523 
1.4596 
1.45909 


1.4832 
1.3477 
1.4756 
14795 


1.4832 
1.4814 
1.4826 
1.4824 


1.4767 
1.4791 
1.4767 
1.4776 


BS Taio 
abe Fs) 
SOGaia 
ISP 7 


she ie) 
592Za9 
ays eh) 
Soa 


DISS 
Sees tS 
SS)5} 1S) 
B9G2 


WATER 
percent 
2.6142 
2.6342 
2.8231 

2.6905 


DTT! 
Zoi 
Z.0008 
2.5184 


2.5436 
2.1694 
1.9308 
2.2146 


2.4213 
2.5302 
2.7816 
2. Td. 


SULPHUR 
percent 
eee) 
2535 
ie27o¢ 
12503 


1.4018 
1.6016 
E3233 
1.3626 


1.7430 
e772 
1.7534 
eZ 


2.0464 
20207 
1.9433 
ZU 059 
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Wt/MOLES X 

g hr/mol Conversion 
31:95 C0367 
15.84 0.0227 
24.47 0.0327 
Gio 0.0336 
23.64 0.0166 
Bo.93 0.0232 
30.60 0.0306 
2376 G:0295 
15.45 0.0223 
1518 0.0218 
23.06 0.0267 
S052 0.0280 
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tars 
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Table 5.2 


Data Set 


S104 
S105 
S106 
S108 
S110 
S15 
$120 
SI25 
S130 
SI35 


Bo 


Sulphur Initial Rate Data 
Pressure = 1.49352 * 105 Pa 
Temperature = 592 K 
Water = 3489 Pa 


Model: Y = A tanh (B X) 


A B AxB pps; sum of sqs 

mol/g hr * Pa 

100 
Oy 385 0.5300 7.3405 266.9 3Sea10n- 
01370 0.6201 8.4954 364.6 4.806* 10> 5 
0.1247 0.8152 10.165 447.3 Paes TG) = 
O09 1s 0.9947 SHO, 600.0 232 el Ome 
0.07081 1.658 11.740 8937/7 TOOS=10a: 
0.05499 2567 14.116 1345 2129%1 05" 
0.04586 S10 16.124 1848 S62 24 0s 
0.03802 4.459 16.953 2043 3.602% 11054 
O10S356 sx Sio}6) isHeles) 2639 1.147*10- ° 
@02937 6:353 teo5g 3077 8.954*10- * 


gf ets5 
C7 @a44.; 
£ ee Me ed 
MeL uee 
mt etn) 

Pas 


a6] 
ecu 


Ht. 


i 


Ai 


il 
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0.125 


(o) 
a 
(o) 
(o) 


Data Set SI04 
Fitted Curve, Y = A tanh(B X) 
O fitted points 


@ data points 


Fractional Conversion of S, (Xs,) 


Pressure = 1.4767 atm 
Temperature = 591.2 K 
0.025 Partial Pressure S, = 266.9 N/m? 


Partial Pressure of water = 3490 N/m? 


Catalyst weight = 0.6079 g 


0 50 100 150 200 250 300 


Wt. Cat./mol feed S, (g hr/mol) 


Figure 5.1 Sulphur Initial Rate Plot 


*100 


Rate (S,) (mol/g hr) 


20 


= experimental data 
Pressure = 1.47 atm. 
Temperature = 592 K 


Partial Pressure of Water = 3490 N/m? 


Weight Catalyst = 0.6079 g 


100 200 500 1000 


Partial Pressure of S, N/m? 


Figure 5.2 Order with Respect to Sulphur 


2000 


61 


62 


5.3 Analysis of WI data 

A total of 51 finite rate data points were obtained in order to derive the initial 
rates used to establish the water dependence of the reaction. Some 8 of these were 
discarded and were not used in the determination of the initial rates. The conversions of 
runs WIO6C, WIO08C, WI10C, WI12C, WI15C and WI15D were all suspected as being 
too low and when they were repeated, this suspicion was confirmed. These runs were 
all performed consecutively and the probable cause of the error was problems with the 
sulphur feeder. It is likely that the “low” range of sulphur syringe gears had been 
erroneously engaged. This would result in less reaction than would be expected at the 
assumed sulphur feed rate and as a result the conversions were correspondingly low. 

The conversion calculated for the run WIO6F is 20% less than would be predicted 
by a smooth curve through the remaining data. An analysis of the raw data indicates that 
the molar feed rate of water was more than 10% less than the average for the run. This 
compares with a variation of +5% for the other water feeds. As a consequence, this 
datum was not included in the initial rate analysis. 

The W! data sets are analysed in the same way as the S| data’) Due to the more 
flexible control of the water flow rate, it was possible to obtain at least five and usually 
six separate data for each data set. Unfortunately however, the proximity to the origin 
was limited by the sulphur flow rates which were fixed. 

A summary of the space time/conversion data is given in Table 5.3 below. Figure 
5.3 illustrates the experimental data for the data set WI01. The initial rates that are 
derived from the experimental data are given in Tabie 5.4. The full fitted data are given in 
the Appendix. 

Figure 5.4 illustrates the logarithmic relationship between rate and water partial 
pressure. The data have been fitted by eye and a curvilinear relationship has been 


assumed. 
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TABLE 5.3 


RUN 


WIO1A 
W101B 
W101C 
WI10 1D 
WIO01E 
WI0 1F 


Avgs 


WIO2A 
WI02B 
WI02C 
WI02D 
WIO02E 
WI02F 
WI02G 
Avas 


WIO04A 
WI04B 
WIO04C 
WI04D 
WIO4E 


Avgs 


DATA SETS WI (summary) 


PRESSURE TEMP 
Atm K 


1.4737 BO Sabe 
1.4710 ais ae 
1.4648 OSS 
1.4753 BIGn to 
1.4761 DoS 
1.4708 he Tie lil 
1.4720 So22 


1.4714 SIosS 
1.4612 Shs ell 8) 
1.4730 Bo2a ne 
1.4722 eps) 1183) 
1.4761 ase 
1.4720 Srolih, he) 
1.4696 spe) hs 
1.4708 oN Pe 


1.4737 Stel ihe) 
1.4652 SoZ aha 
1.4730 Domi o 
1.4722 BI3.13 
1.4761 BOS. lia 
1.4721 592.4 


WATER 


percent 


1.1436 
O.98i15 
She iais 
1.038 1 
1.0070 
0.9819 
1.0266 


2.2368 
ESii35 
1.9366 
1.8204 
2.0106 
1.8234 
£9074 
PISoS 


3.9481 
S3993 
SASS 
3.3023 
3.6638 
3.6820 


SULPHUR 


percent 


0.4846 
0.4838 
0.5040 
0.5048 
Gb0S5 
GS ts 
0.5006 


0.4906 
0.49898 
0.5031 
C505. 
0.4986 
0.5141 
0.3945 
0.4865 


0.4905 
0.4953 
G50 ve 
0.5024 
0.5063 
0.4999 


63 


Wt/MOLES X 

g hr/mol Conversion 
47.52 Cals S7 
SOeG7 Os 
2/69 0.1114 
17.46 0.0957 
3707 0.0843 
SXSL ess, 0.1208 
24.55 0.1242 
19.46 0.1044 
14.25 0.096 1 
9.84 0.0756 
6.89 0.0587 
13:60 0.0987 
14.92 9.0995 
RSH 0.1103 
10.24 0.0964 
7.438 0.0838 
5.46 0.0636 
3.83 0.0501 


- 1 7 hoe. 
7 rel JAE D ord thy oT OF 7 reset 
aida 8 


Table 5.3 continued, 


RUN 


WIO6A 
WIO06B 
WIO6C 
WIO6D 
WIO6E 
WIO6F 


Avgs 


WIO8A 
W108B 
WIO8C 
WI08D 
WIO8E 
WIO8F 


Avgs 


WI10A 
W110B 
WI10C 
W1I10D 
WI110E 
W110F 
W110G 


Avgs 


PRESSURE TEMP 


Atm K 

1.4648 59S 
1.4531 SS ake) 
1.4742 sper ets) 
1.4664 ire), (la) 
1.4702 Bg sr5 
1.4696 5S inl5 
1.4637 pg22 
1.4652 Se LS 
1.4519 SSaa> 
1.4730 59 leo 
1.4687 Sh) 
1.4772 59Sso 
1.4696 Domes 
1.4666 592.0 
1.4628 ssp tS) 
Pace bo iat5 
1.4648 SSK IG) 
1.4746 Oza 
1.4761 SS sao 
1.4720 fo}° ahs) 
1.4708 aise ie) 
1.4695 5915 


WATER 
percent 
5.7800 
5.4334 
Sp eho FA J 
5.2939 
Shenley. 
4.9378 
bocae 


7.0672 
S:s0sA 
7.8088 
7.7498 
7.5824 
T2660 
7.4338 


ST FIAT 
9.0704 
9.4982 
9.6183 
816578 
9.2684 
93607 
9395" 


SULPHUR 
percent 
0.4968 
Oanliom 
0.5039 
0.51038 
0.5096 
0.5189 
0.5076 


0.4981 
0.5788 
0.5039 
0.5671 
0.5082 
O24.63 
0.587 


0.5011 
OSIRIS) 
0.5041 
OS Gil 
0.5100 
OS183 
0.5170 
ORioe 


Wt/MOLES X 
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g hr/mol Conversion 
9.48 0.11416 
6.57 0.09410 
5.09 0.04214 
gas OOSSA7 
2:55 0.04911 
54.0 0.06596 
TAT 0.09028 
5.29 0.07833 
Son 0.04519 
2.58 O0Sass 
1.86 0.04358 
aan0 0.06404 
6.330 G11130 
SReAS, 0.07023 
2920 0.03694 
SHne 0.04455 
1.465 0.03167 
seed 0.06786 
3.002 OOS7o7 


et 7 i 7 4 7 
i . : 
aS J 
7 


; oes ar (ee 


Table 5.3 continued, 


RUN 


WI12A 
W112B 
WI12C 
W112D 
W112E 
W112F 
W112G 


Avgs 


WI15A 
W115B 
W115C 
W115D 
WI15E 
WI15G 


avgs 


PRESSURE TEMP 


Atm K 
1.4616 Bao 1s 
1.4554 Shs Ps. fe) 
1.4730 eke 3. 1h 
1.4746 592e15 
1.4749 Sih 1S 
1.4708 She) aie 
1.4696 DO Tals 
1.4692 boZa0 
1.4616 SOSus 
1.4554 5S lao 
1.4730 Se smio 
1.4722 SES Ks, PS 
L425 BOSE 
1.4696 So lao 
1.4677 Bol 


WATER 
percent 
11.6000 
11.4576 
11.4909 
TARA AS 
iiroooG 
12.3749 
12.1966 
162s 


14.2480 
14.3046 
14.0629 
14.1088 
14.8363 
1.37828 
14.318 


SULPHUR 
percent 
0.4872 
O.51'85 
0.4967 
OS ii2 
0.5134 
0.5096 
05152 
05126 


0.4875 
0.5180 
0.4950 
0.4975 
Ora sS | 
0.4930 
0.5029 
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Wt/MOLES X 

g hr/mol Conversion 
4.702 0.09142 
3.134 0.0665 1 
Za 1e\ 0.03074 
fea 5 0.0383 1 
12200 0.03285 
2.879 0.06362 
2296 0.05043 
Shere) 0.07816 
eae 0.06148 
1.94 0.02909 
1226 0.01844 
Ose 0.02700 
te 5) 0.04322 


=1% 3 bi 
ace*t.! 


GBI Suni 


Table 5.4 Water Initial Rate Data 
Pressure = 1.48810 * 105 Pa 
Temperature = 592 K 


Sulphur = 756 Pa 


Model: Y= A tanh (B X) 


Data Set A B A*B PP water sum of sas 
mol/g hr * Pa 
1000 
W101 Oni 62 0.05544 63965 We Glers PEGS VIA) 
W102 Oz 18 O073827 8.9243 2880.2 2.098% 10 - 
W104 Ova S 7, OF ee 7; 13.4902 5479.1 122308 
WIO06 O1272 0.1496 To.0250 S230 1S37*10% 
W108 0:08271 0.2442 22.6400 11062 2603710 
W110 0.1883 0.1046 haie962 13983 1.186*10 - 
WI12 O W198 O205i 2a. 57141 17594 LED /* 10" 
WI15 0.1063 0.2447 26 1582 2107, 3.402*10 - 


Fractional Conversion of Water-(Xw) 


0.100 


0.075 


0.050 


0.025 


0 10 20 


Data Set WI01 

Fitted Curve, Y = A tanh (B X) 
O = fitted data 

EB experimental data 
Pressure = 1.47 atm 


Temperature = 592 K 


Partial Pressure S, = 756 N/m? 


Partial Pressure of Water = 1536 N/m? 


Catalyst Weight = 0.6079 g 


30 40 50 


Wt. Cat./mol feed (water) (g hr/mol) 


Figure 5.3 Water Initial Rate Plot 
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AA ne . oe Cot e °. 
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Figure 5.4 


Rate (water) (mol/gm.hr) * 1000 


4 = experimental data 
Pressure = 1.47 atm 
Temperature = 592 K 


Partial Pressure of S, = 756 N/m? 


Catalyst Weight 0.6079 g 


1000 2000 5000 10,000 20,000 


Partial Pressure of Water (N/m?) 


Logarithmic Plot of Rate versus Water Partial Pressure 
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Rate (water) (mol/g hr) * 1000 
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b 2 
Y=a(Pw) /(1+¢ Pw) 


— = experimental points 
Pressure = 1.47 atm 

Temperature = 592 K 

Partial Pressure of S, = 756 N/m? 


Catalyst Weight = 0.6079 g 


0 eG 10 “15 20 


Partial Pressure of Water (N/m2)* 0.001 


Figure 5.5 Rate (water) versus Water Partial Pressure) 
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Figure 5.5 illustrates the relationship between the partial pressure of water and the 
reaction rate. This curve is the fit given by the modified damped least squares method 
(Meyer and Roth (58), see Appendix 3) of the function : 
b =p 
TS eera (EW) ae tds C(Py))) gimmie eee ee oat) 
where n = 0,%,1,2,3 
In order for the computer program to run successfully, it was necessary to scale 

the initial rate data so that the dependent and independent variables were of the same 
order of magnitude. This was done by multiplying the rate of reaction by 10 and dividing 


the partial pressure by 10000. The actual function being fitted is therefore ; 


b n 
Rea RW On a(Pwy 104) eit tet (i104) 9) sen Ce eee (5:2) 
Table 5.5 shows the values of the parameters that were obtained after descaling the 


results of the program. 


Table 5.5 Comparison of Rate Equation Parameters 

n a k,*107 b c(K,* 104) Sum of Sas. 
0 0F1925 117.84 0.5228 0.0000 DroZIer Oe 
We 0.2839 28367, O77 96 1040 O32 10: 
1 013557, Ms238 0.8790 0.7461 GOUie 10-8 
2 053358 6.9140 0.8911 0.2742 5.449 *10- 
3 OS 215 6.1696 OS8754 Oule 7.6 B17 Oras 


The choice of scaling factors affects the solution obtained. However, a thorough 
investigation of the effect of different sets of scaling factors determined that those 
chosen gave the lowest sum of squares. An alternative solution reached by scaling the 
partial pressure by a factor of 100 and the rate by 1000, results in an order for water 
of about 2. The log-log plot of this data (Figure 5.4) does not give a straight line, 


indicating the retarding influence of the denominator term on the rate. The slope of this 


9 it Paes 
ai Pore =P 
| iq 


Zo2 ta ree ti 


MObe tet G0) 
“Oi« “9as sav u 
-aiecaas EATS O* 


-“}1 = OFT-€ 


ri 


line as Pw --> 0 will be the order of the reaction with respect to water. Taking the first 
three points of the plot, the slope is about 0.7, and this permits both orders of 1 and 2 
to be valid. The value of 1 is the more probable. 

The parameters estimated by the program are the "a", "b" and '"c" in equation 5.3 
above. The parameter "a" is equal to the rate constant times the sulphur partial pressure, 
(k, Ps , where "x" has been determined to be 0.39, from the SI data analysis). Thus the 
"k," in Table 5.5 is the rate constant of the reaction as calculated from the partial 
pressure of sulphur. Parameter "b” is the order of the reaction with respect to water. 
Parameter "c” corresponds to the water adsorption equilibrium constant. It should be 
borne in mind that this parameter is temperature dependent. No attempt was made to 
determine the nature of this temperature dependency in this work. 

A listing of the main program IRANAL, and the subroutine NLREG is given in 


Appendix 3. 


5.4 Analysis of Temperature Data 

Two additional data sets were obtained at higher temperatures. A summary of the 
experimental data is shown in Table 5.8. Note that the data set WT6L is the same as the 
set WIO6. Table 5.7 lists the initial rate data) The data was analysed in light of both of 
the previous data and the rate constant derived from that data. 

The rate constants were calculated from the initial rate data by substituting the 


values of the kinetic parameters obtained in the rate equation, viz. 


ko = Rw is 0.274 24Pw) 7 Psa)? eae (553) 


The constant in the rate equation denominator is also temperature dependent but is 
assumed to be constant. Table 5.8 shows this data. Figure 5.6 shows the plot of log k, 
against reciprocal temperature. This figure indicates the straight line relationship that 


exists between these quantities. 
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Table 5.6 
RUN 


WT6HA 
WT6HB 
WT6HC 
WT6HD 
WT6HE 


Avgs 


WT6MA 
WT6MB 
WT6MC 
WT6MD 
WT6ME 


Avgs 


WT6LA 
WT6LB 
WT6LC 
WT6LD 
WT6LE 
WT6LF 
WT6LG 


avgs 


Space Time Data at 654 K, 623 K and 592 K 


PRESSURE TEMP 


Atm 

1.4676 
1.4712 
1.4700 
1.4629 
1.4712 
1.4687 


1.4709 
1.4709 
1.4650 
1.4639 
1.4639 
1.4670 


1.4648 
1.4531 
1.4742 
1.4664 
1.4702 
1.4708 
1.4696 
1.4637 


K 
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624.15 
624.15 
624.15 
624.15 
624.15 
624.3 
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WATER 
percent 
6.0065 
S708 
me S99 
5.9485 
pre0Ss 
6.6557 


62729 
6.1704 
b:6092 
ep ehs aye! 
5.8065 
5.8908 


2), 7 O}e, 
5.4334 
Do9a/ 
Begod 
SG t6Z 
5.67S2 
4.9378 
eiayer ke, 


SULPHUR 
percent 
0.2882 
0.2948 
O72 
0.3007 
0.2971 
0.2956 


0.3668 
0.3788 
0.3800 
0.4016 
0.3780 
0.3810 


0.4968 
Otay, 
0.5039 
O:51'03 
0.5096 
Oso 
Ob 183s 
0.5080 


Wt/MOLES X 

g hr/mol Conversion 
S79 232741 
siete O1e592 
453 SUSE E RS. 
290 0.14256 
Pap os ON 2032 
Con Os 9677, 
5.49 0.15534 
4.57 013203 
EL Ele) 0.10042 
2.29 GO7618 
9.48 0.11416 
6:57 0.09410 
5.09 0.04214 
3.45 GODSZ7, 
2.04 0.04911 
6.30 0.06264 
5,70 0.06596 


Table 5. 


Data Set 


WTE6L 
WT6M 
WT6H 


Tabie 5. 


ih 


A 


0.127 
0.244 
O2229 


8 


Rate W 
mol/g hr 


0.019040 
0.03286 
0.056562 


Temperature Initial Rate Data 


Pressure = 1.4670 * 105 Pa 


Temperature = 592 K, 624 K and 653 K. 


B A*B EE oo: 
mol/g hr Pa 
*100 
1 0.1498 1.9040 75s 
1 0.1346 3:2856 566 
6 0.2437 5.5954 440 
Activation Energy Data 
Water Sulphur k In(k) 
pp Pa pp Pa VOT) 
8204 VAS 7.0154 =4a0170 
8756 566 13.088 -13.546 
B79 440 24.848 12905 
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PP water sum of sqs 
Pa 
8204 ESi17#107 
8756 2.454* 10-5 
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Figure 5.6 
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Taking the natural logarithm and calculating the slope of this line gives the quantity —E/R, 
the ratio of the activation energy to the gas constant. The intercept of this line is equal 
to the natural log of the Arrhenius parameter A. 


Thus; 


sein = ~/994 Kk 
AGe RO DOO es 2.) oft Ie oe 20 eee ean (5.4) 


5.5 Thermodynamic Equilibria 

The first 51 experimental runs were performed at low space velocities and high 
catalyst loadings. Under such conditions, the reaction goes virtually to completion and 
the equilibrium state is approached. These runs were basically of a preliminary nature, 
being used to investigate the operating characteristics of the experimental system. They 
are therefore subject to greater error than the later experimental runs, and there is 
considerable scatter. In particular, the measurements of temperature and compositions 
contain more uncertainty and these will affect the calculated equilibrium to a greater 
extent than, for example, the uncertainty in the actual flow rate, since the flow rate has 
little effect on the conversion in this regime. 

Figure 5.7 illustrates the data used to estimate the experimental equilibria) The 
experimental equilibria are estimated by a "free—hand” fit of the data, since the data 
displays a high degree of scatter. However, the trend of the data is quite clear. The 
conversion of water has been plotted in comparison with the thermodynamically 
calculated equilibria at the same temperature pressure and feed composition. Table 5.9 
illustrates the various equilibria that were estimated from Figure 5.7 (Exptl Xw) and those 


obtained from the thermodynamic calculations (Equil Xw). 
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Table 5.9 


Comparison of Equilibria 


Sulphur(8)% 
0.1654 
0.2239 
0.2758 
0.4000 
Go21s 
O6972 


Exptl Xw 
0.180 
0.190 
0.205 
0.210 
0250 
0.280 


Equil Xw 
0.128 
0.132 
0.138 
0.147 
0.156 
0.161 
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5.6 Discussion 


5.6.1 Discussion of Errors 

The study of the kinetics of a reversible reaction can be made from either side of 
the equilibrium. When the equilibrium of the reaction lies well to one side at the 
temperature of interest, the kinetic study is complicated by the necessity of obtaining 
precise data on the minority products which may be present only in small very amounts. 
In the case of the Claus reaction, the equilibrium eonversion of H,S is around 90% at 600 
K so that the maximum mole fraction of H,S is only about 0.1. In an initial rate study, the 
data must be obtained as far from equilibrium as possible and so in the case of this 
reaction, the reaction products H,S and SO,, are less than 1 % of the reactor product 
stream. As a result, the major uncertainties associated with the collection of the 
experimental data are contained in the product analysis. The average percentage 
uncertainty in the measurement of the G.C. peak areas are for H,S, SO, and H,O 
respectively, 2.46, 3.38 and 2.65 %. The uncertainties in the other experimental 
measurements are small, being less than 1 %, except in the case of the Differential 
Pressure recorder which included an uncertainty in reading the instrument of about 4 % at 


low scale deflections, viz, 
Scale reading 5 + 0.2 —~> 4 % error 


Since the volumetric flow rate is proportional to the square root of the scale deflection, 
the actual error is 2 %. 
Similarly : 


Pf, Scale reading 9.0 + 0.1 —~> 1.1 % error 


Patm, ” ‘ 700 mmHg + 0.1 —> 0.014 % error 
Pree 46.0 + 0.2 —> 0.4 % error 
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The average uncertainty in the sulphur and water feeds are 1.8 and 1.7 % respectively. 

The most significant uncertainty in the data was caused by the operation of the 
water and sulphur feeders. The operation of the water feeder alone produced a steady 
flow, as evidenced by the fact that the repeatability of the G.C. peak area test (Section 
4.2.5), which determined that the fluctuations due to the water feeder were 2.87%, was 
similar to the accuracy determined by the GC. calibration, (2.65%). The sulphur feeder 
calibration data (Appendix 1), show that the sulphur flow uncertainty was not affected by 
the choice of flow rate. As a result, it is not possible to explain the unevenness in the 
water product by considering the water and sulphur feeders individually. The source of 
the uncertainty in the product analysis must derive from either the use of the two 
feeders at once, or the effect that the vapourisation of sulphur had on the flow. Since it 
was not known whether the sulphur feed actually affected the water feed or not, the 
water feed was determined in the analysis by back calculating from the products, using 
stoichiometry and the G.C. calibrations. This calculation will not eliminate the error 
involved in the analysis since only the H,S component was used to calculate the H,O in the 
feed , and the overall mass balance could still be checked using the SO, component. 

There are two methods that may be used to estimate the uncertainty associated 
with the data ; 

13 Calculation of the H,S/SO, mole ratio. 
za Calculation of the oxygen g-atom balance over each run. 

The H,S/SO, mole ratio throughout the course of this study was within 5% of 2, 
the stoichiometric ratio. This ratio determines that the Claus reaction stoichiometry is 
being followed. The calculation of the product sulphur and the feed water using this 
stoichiometry is therefore valid. 

The difference in the g-atom balance on oxygen is around 3%. This indicates that 
the mass balance pees the reactor is correct since there was no source of oxygen 
except in the feed water and the oxygen balance has been calculated from the product 
mole flow rates assuming Claus reaction stoichiometry. The g-atom balance on 
hydrogen produces slightly larger differences than the oxygen balance, about 3.5%. This 
figure illustrates the error associated with calculating the mole flow rates from the H,S 


and SO, GC. peaks 
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As would be expected, because the water feed rate is increasing in the WI data 
sets, the percentage uncertainty decreases with increasing water feed, since the relative 
size of the H,S component to the water is much less. 

Alternative sources of error include questioning the validity of using the Gibbs 
free energy minimisation technique to estimate the equilibrium distribution of the sulphur 
species, especially as these results have been proven to be inexact. There is as yet no 
other method of directly measuring the sulphur species distribution and as a 
consequence, no method of measuring the uncertainty of this method. It is assumed to 
be small in comparison with the errors associated with the product analysis. 

It is evident that there is only a small amount of error associated with the mass 
balance across the reactor. However, an additional source of error is introduced by the 
uncertainty in the feed flow rates of water and sulphur. The average variation in the 
water feed rate for the SI data sets was 9.4%. The variation in the sulphur flow rates 
was 3.6%. The corresponding figures for the WI data sets were 4.0 and 2.9% 
respectively. The much smaller error in the water feed rates for the WI data sets is due 
to the much higher water feed rates employed. The average uncertainty in the S! data is 
therefore about 10% and in the WI data set, 5%. How these uncertainties affect the initial 
rates is not simply determined since each individual datum has a variable influence on the 
slope and shape of the fitted curve according to its position. To estimate the uncertainty 
in the initial rate data it is simpler to look at the data itself. 


The variance of the slope of the line is given by, (40), 


Sbi={Syixyey s(ci=o)2) 0 Si etary e me eet eee (5.5) 


and the confidence band for the slope is given by, 


epeiag Sj) 


where t is obtained from the "t” distribution. In the case of the log-log 
plot of the SI data, Figure 5.2, with 10 data points, there are 8 degrees of freedom. At 


the 90% confidence level, the slope of the line is, 
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0.397 + 0.047 (12% 


The errors associated with the WI data sets are less than those with the SI data 
sets. The average fit variance of Figure 5.5 is 2.6%, the maximum variance being in the 
point closest the origin. Such a low variance is a result of the fact that 3 parameters 
were used to fit the data. This is evidenced by the fact that only a small improvement in 
the value of the sum of squares is achieved with increasing power of the denominator 
term. The uncertainty in the finite rate data is 7%, so that the uncertainty in the derived 
kinetic parameters must be at least 7%. As a consequence, the values of the parameters 
will be rounded off to two significant figures. 

The WT data are subject to the same magnitude of error as the WIO6 data set. 
The rate constants are therefore subject to a compounded error of 15% (85). The 
activation energy is calculated from only three points, the linear alignment of which is 
plainly fortuitous, as can be seen from the foregoing analysis. The error in the activation 
energy may be calculated using the maximum and minimum values, obtained from the 


maximum error in the k, values of 15 %, thus, 


ER 7 9947 1] 200K 


5.6.2 Sulphur Order 
The order of the reaction with respect to sulphur was determined to be 0.397. 
This value was obtained by performing a least squares analysis on all the sulphur initial 


rate runs. Since the partial pressure of water was 3721 Pa, stoichiometric restrictions, 
2 H,O +3/8 S, = 2H,S + SO, 

might be thought to limit the range of the validity of the data for sulphur order derivation 

to a maximum sulphur (8) partial pressure of (3721 * 3/16) Pa, i.e. 698. In taking the 


logarithms of the rate equation, the following relationship is obtained : 


log(—Rs) = log (k) + a log (Ps) + log [f(Pw)] 
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and so the use of all the initial rate data is justified in obtaining the order of the reaction 
with respect to sulphur, since the water partial pressure term does not affect the slope 
of the line. 

The order with respect to sulphur may be rounded off to 0.4. The orders of the 
forward reaction have been found to be equal to one half the stoichiometric coefficient 
of the species as given by equation 2.7. Although the forward rate equations are 
essentially empirical , it is reasonable to expect that such a relationship would be found 
for the reverse reaction orders, since there are thermodynamic restrictions which 
govern the form of the kinetic equations. These are described by Denbigh (24). The 
order with respect to sulphur that would be expected following this relationship, is 
0.1875. The experimental order is almost twice this. Indeed, it would require an average 
sulphur species of about 4 to make the order 0.4 by this relationship. Since there is a 
distribution of sulphur species, it is unreasonable to presume that only S, is reacting with 
water. This shown by the Comparison of Rates section below (5.6.5). If it is assumed 
that there is decomposition of the sulphur molecules, then the actual Claus reaction may 
be taking place with sulphur whose average number of atoms is around 4, even although 
the species S, is but a small fraction of the sulphur vapour. There are thermodynamic 
limitations which will determine the nature of the sulphur species on the surface but as 
has been said previously, the exact distribution of sulphur species in the vapour is 
unknown, far less on a solid surface. By this reasoning, the order of the reaction with 
respect to sulphur would be temperature dependent, increasing with increasing 
temperature. This effect could be used to determine whether the above argument is 


valid. 


5.6.3 Water Order 

Previous work (30,53,56) had established that the rate of the forward Claus" 
reaction was retarded by the effect of water adsorbing on the catalyst surface. It was 
to be éxpected, then, that this effect would be observed in this study, and that, as a 
result, it would not be possible to determine the order of the reaction with respect to 
water simply by holding the sulphur partial pressure constant and varying the water partial 


pressure. Therefore, the expected form of equation is fitted to the initial rate equation 
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and the kinetic parameters are estimated using a least squares analysis. 

Table 5.5 contains the results of this technique. The "best fit’ of the data is given 
in the case of n=3, which is marginally more exact than the case of n=2. These two 
functions are apparently significantly better than the next best fit of n=1. However, the 
application of the statistical method of discrimination known as the "F” test determined 
that it was not possible to discriminate between the fits given by n = % and higher. The 
lack of fit of each model is characterised by the sums of squares of the deviations, 
hence, the ratios of the sums of squares of the residuals to each model to that of n=3 


(SSi/SS3) are, 


SS0/SS3 SS 21SS3 Sou/Sss SS2/SS3 
4.50 2.00 33 OS 


At the 90% confidence level and with 7 degrees of freedom, the test statistic is 2.78. At 
this confidence level, therefore, it is only possible to state the water is inhibiting the 
reaction, and that it is not possible to determine the order of the inhibition. 

The derivation of these parameters is a matter of some subjective judgment, 
since numerical techniques cannot guarantee that the global minimum least squares 
estimate of the parameters will be found. A number of initial guesses were employed 
and each minimum found, was compared both to the previous minima and the initial guess 
which achieved the lowest minimum was noted. 

lt was necessary to scale the data in order for the method to converge. It would 
converge rapidly and reliably for these functions when the maximum value was below 10. 
For this reason, the rates were multiplied by 10 and the partial pressure divided by 
10000. The allowance that must be made for the scaling values is illustrated in 
Appendix 2 in the Calculation of rate constants section. 

Uncertainty in the data for the water runs is less than that contained in the sulphur 
run data. This was probably due to the low flow rate of sulphur employed and this would 
lessen errors due to sulphur vapourisation and equilibration uncertainties. In addition, the 


higher water flows employed allowed a larger margin of absolute error. The greater 
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accuracy permits the simultaneous evaluation of three kinetic parameters with some 
degree of confidence. 

Allowing for the 7 % error that attends the fit of the model to the water initial 
rate data and noting that the tendency of the model will be to reduce the order of the 
reaction with respect to water, the order of the reaction with respect to water may be 
rounded off to 1. This agrees with the one half times the stoichiometric coefficient 
relationship that is suggested by studies on the forward reaction. This lends a measure 
of credibility to the suggestion that the sulphur order may be related to a lower order 
sulphur molecule than S,. The log-log plot, Figure 5.4, although showing a curve 
indicates that the order with respect to water is greater than 0.7, since this is the limiting 
slope of the experimental data) This supports the non-linear regression result that the 
order is 1. This is evidence that the parameters derived by the method are not spurious. 

The denominator parameter K,, has the value 2.472 * 10 -* where the units are 
(Pa). The corresponding values from Liu and McGregor’s work on the forward reaction 
are 0.00600 and 0.00433 respectively, units (mm Hg). In comparable units these are 
egivalent to 4.5 * 10-5 and 3.18 * 10-°. The adsorption terms are of a similar order of 
magnitude .The term is temperature dependent since it corresponds to an adsorption 
term, and so it is expected that with increasing temperature it will decrease. This work 
was performed at 592 K whereas the previous work was carried out at about 530 K, so 
this trend is reflected. This is in spite of the fact that three different catalysts, Porocel, a 
commercial bauxite, and 2 gamma-aluminas, Alon and Kaiser S~201 (this work), are 


employed. 


5.6.4 Temperature Variation 

Figure 5.6 illustrates the variation of the rate constant with temperature. The 
activation energy derived from this is 66.5 kJ. Only three temperatures were available 
and so only three data points are used to derive the activation energy: this value, 
therefore, has a limited range of validity. The lower temperature was limited by the dew 
point of the sulphur vapour because it was desirable to operate significantly above this in 
order to ensure that condensation of sulphur would not occur. The maximum 


temperature was limited by the high temperature cut-out protection of the reactor. This 
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could have been adjusted but would have required a significant rebuilding of the 
controller — an undesirable step at the late stage in the study. 

Calculation of the rate constant at the lowest temperature can be done in two 
ways: either directly from the measured rate of reaction, or else from the a” parameter 
derived from the water data reduction. 

k(rate) = 6.4606 * 10” 
k(a) = 7.0154 * 1077 


The difference, which is 8 %, measures the difference between the rate function, which 
includes the reaction order with respect to sulphur, and the experimental function. This 
is a measure of the consistency . 

The degree of uncertainty that should be ascribed to the activation energy is a 
function of the accuracy with which the rate constant is measured. This is estimated to 
be about 15 %, which leads to approximately 26 % error in the activation energy. 

The activation energy may be compared to the activation energies for the 
forward reaction. These have a value of 30.1 kJ in the case of McGregor’s studies and 
30.75 for Liu’s work. Since the Claus reaction is exothermic there ought to be a higher 
activation energy for the reverse reaction. Allowing for the degree of error that is 
associated with the activation energy derived in this study, these are of comparable 
magnitude. The reverse reaction activation energy ought to differ from the forward 
activation energy by the heat of reaction which is in this case is about 100 kJ/mol, so 
there is plainly a large discrepancy involved, which is not explained by the fact that a 
different catalyst is being employed. The error associated with the activation energy 
calculation is large however, the maximum value of the activation energy is only 80 kJ 
which is 50 kJ less than the value that is implied from the forward reaction activation 
energies and the heat of reaction. Obviously, it would be advantageous to repeat the 


initial rate runs using more temperatures in an effort to re-assess this data. 
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5.6.5 Comparison of Rates 

It would be expected that the calculated rates with respect to either component 
ought to be in stoichiometric ratio to the other since only the Claus reaction is taking 
place. However, this is not the case since the rate with respect to sulphur (8) neglects 
the components of the overall sulphur reaction rate that is due to the other sulphur 
species. The rates of all the sulphur species were calculated and the overall 
stoichiometric ratio of water and sulphur rates was calculated. On average, there was a 
difference of 12.8% in the rates. This difference was found to be correlated with the 
water feed rate. 

A positive deviation shown in Table 5.10 reflects a higher water reaction rate than 
that expected from simple stoichiometry. The higher the water feed rate, the more the 
water reaction rate exceeds the stoichiometric rate. This relationship is probably due to 
uncertainties in the product analysis, particulaly in the water component. It is probably 
exacerbated by the calculation of the water feed rate from the hydrogen sulphide 
component, since the error is probably associated with the tailing of the water peak at 
high mole fractions, and this error would have been minimised if the water feed rate had 
been measured directly. This deviation represents the limit of accuracy of the 
experimental data, ie. 13%. Since the error is not random, it will have an effect on the 
calculation of the kinetic parameters. Since the rate is exaggerated at high mole 
fractions, the inhibiting effect of water will be reduced. There will also be a tendency to 


increase the apparent reaction order with respect to water. 
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Table 5.10 Comparison of Sulphur and Water Rates of Reaction 
Data Set % Deviation Data Set % Deviation 
S!i04 alias WIO 1 19:0 
S!I05 SAT) WI02 =™3.0 
S106 -14.5 WI04 =e5 
Sl08 =O WI06 +10.0 
Siig Hee WI08 we eds) 
Sis 0:0 W110 15:5 
S!20 = 1 20 W112 +18.0 
Sizo -o.0 W115 H1e.5 
Si30 ay Wee 

S135 ralvies 


Appendix 2 contains an example calculation of the rate of the forward reaction. 
This is compared to the rate of the reverse reaction at identical conditions. The rate of 


the forward reaction is 9.75 times that of the reverse reaction. 


5.6.6 Comparison of Equilibria 

As shown by Table 5.9, the thermodynamic equilibria are consistently below the 
experimentally determined equilibria) This is the same relationship that was noted for the 
approach to the equilibrium for the forward reaction. It is not possible therefore, to 
explain the discrepancy between the equilibria in terms of the sulphur species data being 
too low or too high. These results indicate that a more complicated mechanism may be 
at work which causes the discrepancy, however, as these data are not very precise and 


are greatly scattered, they can only be used as an illustration that this relationship may be 


found. 
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6. SUMMARY - CONCLUSIONS AND RECOMMENDATIONS 


Conclusions — Kinetic results 


de The empirical rate equation for the Claus reaction indicated by this study is ; 


a b 
—Rw = A exp(-E/R T) (Pw) (Ps,) /[1+K, PwP 


where, 
A = 0.5004 
-E/R = -7994 
a= 0.89 
b = 0.39 
K, = 2.742 * 10° 


Temperature in Kelvin, pressure in Pa, rate in mol/hr g cat. 


2 The order of reaction with respect to water is in agreement with the order 
required by thermodynamics. This is not the case with the sulphur component. 

S Inhibition by water is confirmed by statistical evidence. It was not possible to 
determine the order of the inhibition. 

4. The relative rates of reaction are within 13% of the stoichiometric ratio. There 
was a correlation between the water feed rate and the error that would depress 
the inhibition and increase the water order slightly. 

By The parameter K,, corresponding to an adsorption equilibrium constant, is of the 
same order of magnitude as those values measured by Liu and McGregor for the 
case of the forward reaction. This indicates that this work is consistent with the 
previous work on the forward Claus reaction. 

6. The activation energy of the reverse reaction is found to be of comparable order 
with that for the forward reaction. The activation energy that is calculated 


indicates that the Reverse Claus reaction is endothermic but the difference in the 
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forward and reverse activation energies is less than that predicted by the heat of 
reaction. 
The experimental equilibria are significantly higher than the thermodynamic 


equilibria. 


Conclusions — Equipment Performance 


The use of two liquid syringe feeders for the reactant feeds had the advantage 
that the feed rates could be adjusted rapidly and precisely. However, the 
combination of the two led to large fluctuations in the water feed, probably due 
to the sulphur feeder. 

The Berty reactor operated smoothly and without problems throughout the 
course of this study, except that a small amount of sulphur back-diffused past 
the impeller shaft seal and condensed on the water cooled shaft. This was a 
negligible amount of sulphur and was readily removed. 

The performance of the sulphur condenser was very sensitive to temperature 
variation. This problem was successively eliminated by adding secondary 
condensers. The condensing system as described in this work performed well 
with only minor plugging occuring in the GC. line. This was related to a 
temperature rise in the main condenser. 

The performance of the gas chromatograph and electronic integrator was in 
general very good. The Chromosorb 108 column provided excellent separation 
of the reactor product components. In particular the water response had less 
"tailing” than had been experienced with other columns such the combination of 
Poropak S and T used earlier. 

The GC. calibration apparatus was awkward to use. !n particular, the piston that 
expelled the gas mixture was extremely difficult to move, even with the aid of a 


jack. The system allows fairly rapid calibration of the GC. but is clumsy. 
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Recommendations 


iff 


The principle source of error in this work is the uncertainty in the feed 
composition. It is likely that both liquid feeders contributed to this, and in 
addition, the vapourisation of sulphur may be partly (or wholly) responsible. 
Further studies should have as a primary objective the establishment of steady 
feed flows, as evidenced by the product analysis, which ought to be constant, and 
consistent with a mass balance employing the calibrations. It is therefore 
suggested that the feasibility of a feed system involving the use of water and 
sulphur saturators be investigated. 

The gas chromatograph is a satisfactory method of obtaining the mole fractions 
of the products on a sulphur free basis. However, with the reverse reaction, 
sulphur is a significant fraction of the gas mixture and it would be advantageous 
to be able to measure the sulphur component in order to complete the mass 
balance. This could be achieved by condensing the sulphur component of a 
sample product stream over a measured period of time. This would permit the 
calculation of the mass flow rate of sulphur 

With regard to the uncertainty in the feed composition, it would be advantageous 
to measure independently the feed composition rather than relying on the feeder 
calibrations. This would also enable the smoothness of the feed stream to be 
monitored. It is therefore recommended that provision be made to measure the 
mass fraction of sulphur in the feed stream. A similar method to that suggested 
for the sulphur product measurement might be employed, including G.C. analysis 
to obtain the water component. 

It was found in the early stages of this work that much of the tubing and vessels 
used in building the early versions of the experimental equipment was catalytically 
active. It is strongly recommended that great care is taken to ensure that only 
316 stainless steel is used in the fabrication of equipment for Claus reaction 
studies and that all tubing and vessels in contact with the reactants be kept 


scrupulously clean. 
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NOMENCLATURE 


Pre-exponential factor 

Initial rate parameter 

Surface Aluminium atom 

order with respect to water 

exterior surface area * shape factor, m? / m3(bed) 
order with respect to sulphur 

initial rate parameter 

rate equation parameter 

specific heat, J/Kg K 

Catalyst particle diameter, m 

diffusion coefficient, m*/s 

activation energy 

equilibrium 

mole flow rate of a, mol/s 

fugacity 

heat transfer coefficient. W/m? K 

heat of reaction, J/kg mol 

rate constant 

rate constant of forward reaction 

mass transfer coefficient, kg mol/m? s 

rate constant of reverse reaction 

equilibrium constant 

equilibrium constant of adsorption term, forward 
reaction 

equilibrium constant of adsorption term, reverse 
reaction 


catalyst bed depth, m 
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he 


Patm 


Ph,s 


Psn 


molecular weight 

mass flow rate, kg/s 

absolute reaction rate, kg mol/s 
mass flux, kg mol/s 

partial pressure, Pa 


unit of pressure, Pascal (N/m?) 


atmospheric pressure, mmHg 


critical pressure, atm 
feed pressure 

partial pressure of H,S, Pa 
Prandtl number 


reactor Pressure 


partial pressure of sulphur, where n is the average 


size of the sulphur molecule 

partial pressure of SO,, Pa 

partial pressure of total sulphur, Pa 
partial pressure of S,, Pa 

partial preassure of H,O, Pa 

heat evolved, J 

reaction rate, kg mol/s kg cat 


rate, g mol/hr g catalyst 


rate with respect to hydrogen sulphide, 


rate with respect to sulphur (8) 

rate with respect to sulphur dioxide 
rate with respect to water 

gas constant 

Reynolds number 

impeller speed 

standard deviation of the slope 
standard deviation of y over x 


temperature, K 
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critical temperature, K 
feed temperature, K 
reactor temperature, K 
chemical potential 
weight of catalyst, g 
conversion 
conversion of a 

mean x 

fugacity coefficient 
viscosity 

density 

shape factor 


stoichiometric coefficient 
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APPENDIX 14 Calibrations 


Calibration of Equipment 


iN Sulphur Syringe 


Figure A~1.1, Sulphur Feeder Gear Arrangement 


Six speeds are possible with the above gear train. These are arranged in 
the following progression : 
Gear A engaged, B disengaged, ——» high range 
Gear A disengaged,B engaged, ——> low range 


Gear S engaged with 1, 2 or 3, ——* 3 speeds in each range 


wed eee are * 


Ta ae 


| 3 A 5 
SOHGUCHEL RET UPARE 
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Calibration was achieved by using the micrometer scale located on the piston rod. 
Using this scale, it was possible to measure the rate of travel of the syringe piston. This 
speed could then be correlated to the weight of sulphur expelled from the syringe over 
the same period of time. The ratio of sulphur weight to distance should be the same for 


each gear. The Table A-1.1 illustrates the calibration data obtained. 


Table A-1.1, Sulphur Feeder Calibration Data 


Gear Mass flow Piston speed Ratio Average flow 
g/hr cm/min g/cm mol S, /hr 
a7 S3027 4.486 13.86 1.4937 
H1 5.3429 77376 14.23 2.3952 
L2 8.8395 $.6619 14.14 30201 
H2 10.570 1.4817 14.10 4.6413 
ES 13.588 1.9420 14.28 SRE 
H3 21.008 2.9729 14.15 9.2763 


The deviation from the mean ratio does not appear to be related to a particular 
gear or range of gears, and so a “corrected” sulphur flow was calculated using the mean 
ratio. 

The calibration data were obtained at an oven temperature of 135 °C. The length 
of time for each sampling was not less than 30 minutes. A single comparison run was 


performed at 173 °C, that is, above the liquid sulphur transition point (75). 


H2 1.0507 g/hr 0.14847 cm/min 14.13 g/cm 


From this result, it can be seen that the sulphur oven temperature has little effect 
on the piston operation over this temperature range. It was necessary to ascertain this 
as the variation of the density of sulphur was not known. In addition, it indicated that the 
performance of the syringe was independent of the sulphur viscosity, since the viscosity 


transition temperature is about 159 °C. 
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Water Syringe Feeder 

The calibration of this piece of equipment was carried out over the 
bottom two ranges of the capacity of the instrument, designated 1/1000 and 
1/100. These ranges covered the entire range of flows to be used in the kinetic 
experiments. The advantage of using a large capacity syringe in conjunction with 
the pump, eg. infrequent refilling, is limited by the ability of the pump motor to 
generate sufficient power to overcome the added friction. As a result of this 
limitation, a 20 mi. syringe was used throughout the period of this study. This 
size was sufficient to accommodate at least one complete run, even at high flow 
rates. 

The calibration was carried out by filling the syringe with distilled water 
and running the pump over various times of up to six hours in the case of very 
low flow rates. The discharge was collected in a narrow necked flask. Care was 
taken to reduce evaporation by using a cotton wool stopper. The loss due to 
evaporation was estimated by applying Fick's law, assuming saturated vapour 
above the water surface. The error was estimated to be less than 1%, which was 
negligible in comparison to the error associated with the pump itself, especially at 
low flow rates, as shown by Tables A-1.2 and A-1.3. 

Eight measurements were made over the entire 1/1000 range, from 100 
to 900 percent maximum flow. Six measurements were made over the 1/100 
range, from 50 to 500 percent maximum flow. The two sets of data were fitted 
by a linear least squares regression. The results are shown in the following 


Tables, A=1-2andiA—1:3 


Reactor Pressure Transducer 

The reactor pressure transducer was calibrated against a mercury 
manometer open to atmosphere. It was calibrated with the reactor system at the 
operating temperature in order that variations due to temperature would be 
minimised. The atmospheric pressure was measured by a barometer. The 


following Table A- 1.5 contains the calibration data obtained. 
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Differential Pressure Cell 

The differential pressure cell was calibrated using a precalibrated dry test 
meter. The system pressure at which the calibration was carried out was set to 
8.5 psig, which was approximately the average operating pressure. The 


calibration data is shown in Table A-1.4 


Gas Chromatograph 

The response of the gas chromatograph to the hydrogen sulphide and the 
sulphur dioxide components of the product streams was calibrated using sample 
gas mixtures of known composition. The apparatus used for this is illustrated by 
Figure A-1.2 and was originally constructed by McGregor(51 ). 

The sample mixing chamber was a five litre capacity lucite cylinder and 
was equipped with a moveable piston and mixing fan. The lucite cylinder had been 
calibrated by McGregor and this was checked and found accurate, thus the 
volume of the cylinder was known as a function of piston position. 

Following the calibration procedure as described by Cho (15), the entire 
calibration apparatus was purged with N, for some two hours before use. The 
cylinder was allowed to fill with pure nitrogen via line A, and was allowed to 
come to room temperature and pressure. The gas burette and line B were purged 
with the gas of current interest and a volume of gas was trapped in the gas 
burette. The height of the mercury reservoir was adjusted to bring the gas in the 
burette to atmospheric pressure, as was indicated by the water manometer. 

The trapped gas was compressed to about 10 psig by increasing the 
height of the mercury reservoir, thus forcing the gas into the lucite cylinder via 
line B. The residual volume left in the burette was measured and the volume of 


gas that had been forced into the cylinder was known by difference. 
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Figure A~1.2 G.C. Calibration Apparatus 
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The cylinder mixture was allowed to mix thoroughly using the fan for 
about 20 minutes. The gas was then forced through the GC. sampling valve by 
raising the piston. Several samplings were taken of each mixture and the average 
peak areas calculated. 

Since only small amounts of these two gases were expected in the 
product stream, all calibration mixtures were kept below 1%. The following 
Tables A-1.6 and 1.7 summarise the calibration data taken. 

Calibration of the water peak area was somewhat different, since the 
above apparatus was useless for high water concentrations. The water syringe 
feeder that has been described above was used in conjunction with the nitrogen 
D/P cell to establish the GC. water response factor. Using the syringe feeder and 
D/P calibrations, mixtures of known H,O0/N, composition were obtained. These 
were mixed using the reactor impeller and fed directly to the GC. An attempt 
was made to calibrate the G.C. water reponse using gravimetric analysis using 
"Drierite” and molecular seive, however, the results were subject to large 
variations. The calibration data obtained are given below in Table A-1.8 

The response factors for H,S and SO, are constant. The response factor 
for H,O is shown to be a function of the mol percent of water, since the 
intercept of the least squares straight line is much larger than zero. (0.09237). At 
small concentrations of water there would be a large error due to the assumption 
of a constant value for the response factor. The calibration has been performed 
over a greater range of water partial pressures than that to be encountered 
during the experimental runs, from less than 1% to greater than 14.5. In order to 
minimise the effect of the variation in the response factor, the water peak 
calibration data were fitted using a quadratic equation. This data is contained in 
Table A-1.9. The mole fraction of water must now be obtained by an iterative 
technique. A response factor is assumed and the mole fraction calculated. The 
response factor is recalculated and the iteration is performed until the mol 


response factor converges. This procedure is described in more detail in 


Appendix 2. 
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Table A-1.2 Water Syringe Feeder Calibration, 1/1000 range 
X = Percentage flow / 100 


Y = Mass Flow, g/min * 100 


The Coefficients of the Polynomial 
A, = —0.10665 
A, = 0.45620 


Calibration Data 


X measured Y observed Y calculated 
1.000 0.3624 0.34955 
2.0000 0.7827 1 0.80575 
4.0000 1.72394 L7S5 
5.0000 2.17974 2.1744 
6.0000 26173 2.55099 
6.0000 S93 3.08678 
8.0000 3.524 3.54295 
9.0000 oysiely © 399915 


Variance = 0.000335 
Standard Deviation = 0.018297 
Maximum Pct Error = 3.548072 


The calibration equation is therefore, 


PCT error 
255 
2.94 
0.34 
0.25 
0.48 
1.04 
0.55 
0.05 


mass flow (g/min) = 4.5620 * 10-5 X(percent flow) — 0.10665 10° 
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Assay) 


Table A-1.3 


Water Syringe Feeder Calibration, 1/100 range 
X = Percentage flow / 100 


Y = Mass Flow, g/min * 10 


The Coefficients of the Polynomial 


A, = —0.09326 


A, = 0.46653 


Calibration Data 


X measured Y observed Y calculated PCT error 
0.5000 0.14293 0.14000 2405065 
1.0000 0.38800 0337326 S/9595 
2.0000 0.82693 83379 1.55478 
3.0000 P29730 1.3063 1 1.16246 
4.0000 1.76930 1.77284 0.19987 
5.0000 2.25310 2235956 0.60977 


Variance = 0.000164 
Standard Deviation = 0.012788 
Maximum Pct Error = 3.798551 


The calibration equation is therefore, 


mass flow (g/min) = 0.46653 * 10° X (percent flow) — 0.009326 
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Table A-1.4 Calibration of Differential Pressure Cell 


X = Chart Percent 


Y = Flow Rate squared, (I/min)? 


The Coefficients of the Polynomial! 


A, = —0.05350 
A, = 0.08998 
Calibration Data 
X measured Y observed Y calculated PCT error 
oe 5U00 0.22966 23895 4.04475 
7.2500 0.59247 0.59889 1.08349 
7.7800 0.63072 0.64388 2.08628 

13.750 1.19806 1.18379 (teh iets) 
18.500 1.64552 160120 2.08484 
2-500 OTSA LeSoiy 1.72629 
22.000 1.89946 192616 1.40536 
29750 2.59748 2.62354 100325 


Variance = 0.000595 
Standard Deviation = 0.024393 
Maximum Pct Error = 4.044746 


The calibration equation is therefore, 


Q (I2/min) = 0.08998 X (chart percent)- 0.05350 32 ........ (A- 1.3) 
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Table A-1.5 Pressure Transducer Calibration 
X = Chart Percent 


Y= Pressure Cm Hg 


The Coefficients of the Polynomial 
A, = 0.01383 
A, = 0.89557 


Calibration Data 


X measured Y observed Y calculated PCT error 
32500 3.00000 2.92444 2.51869 
10.000 9.00000 8.96954 0.33841 
32.000 28.4000 28:672)1 0.95808 
40.000 36.0000 35 :3567 0.45374 
607750 54.4000 54.4198 0.03632 
65.000 58.0000 Sa2209 O:38955 
70.500 63.4000 fey}. eh! 0.39184 


Variance = 0.036751 
Standard Deviation = 0.191707 
Maximum Pct Error = 2.518685 


The calibration equation is therefore, 


Pr (Cm Hg) = 0.89557 X (chart percent)+ 0.01383 ....... (A- 1.4) 


it +a 
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Table A-1.6 Calibration of Gas Chromatograph, Hydrogen sulphide 


X = (H,S/N,)area 
Y = (H,S/N,)mol 


The Coefficients of the Polynomial 


X measured 


0.06540 
0.08110 
0.08753 
0.09820 
0.15440 
0.16160 
0.16510 
0.26600 
0.29430 
0.41690 


A, = —0.00198 


A, = 1.84419 


Calibration Data 


Y observed 
0.12840 
0.14860 
0.16290 
0.17500 
0.26990 
0.30210 
0.30430 
0.47700 
Vaso 
0.77720 


Y calculated 


O71 1S6s 
0.14759 
0.15944 
O73 12 
0.28277 
0.29604 
0.30250 
0.48858 
0.54077 
0.76686 


Variance = 0.000065 


Standard Deviation = 0.008073 
Maximum Pct Error = 7.606923 


The calibration equation is therefore, 


Response Factor H,S = 1.84419 


PCT error 
7.60692 
0.68208 
ZZ TS 
2.35547 
4.76669 
2.00482 
0.59245 
2.42695 
0.72026 
e2zoee 
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Table A-1.7 Calibration of Gas Chromatograph, Sulphur Dioxide 


X = (SO,/N,)area 
Y = (SO,/N,)mol 


The Coefficients of the Polynomial 


A, = 0.00508 
A, = 1.20770 
Calibration Data 

X measured Y observed Y calculated PCT error 
0.09162 0.11260 OAS 7 2.77797 
0.12760 0.16820 O.T59 16 5.36201 
0.14700 0.17250 0.18261 5.8616 
0.20310 0.25660 0.25036 2.43080 
0.25050 C3202 0.30761 4.38058 
0.25700 0.30920 0.31546 2.02384 
0.27630 0.32810 033877 S25093 
0.37100 0.44870 0.45314 0.98856 


Variance = 0.000079 
Standard Deviation = 0.008881 
Maximum Pct Error = 5.861156 


The calibration equation is therefore, 


Reponse Factor SO, = 1.20770 


Table A-1.8 Calibration of Gas Chromatograph, Water 


X = (H,O/N,)area 
Y = (H,0/N,)mol 


The Coefficients of the Polynomial 


A, = 0.09237 
A, = 1.04680 
Calibration Data 

X measured Y observed Y calculated PCT error 
0.7648 0.86370 0.89296 S277.0 
1.2280 1.02530 1.37248 3.43216 
2.1471 2.21600 2.34005 5.29994 
2.8680 2.93040 3.09460 5.30658 
S.ce2d 3.46080 3.58200 3.38499 
6.3273 7500950 7.23924 3.37341 
9.5966 10.9890 Om aoa Sues df, 
13.830 13.8340 14.5696 4.36255 


Variance = 0.035997 
Standard Deviation = 0.189729 
Maximum Pct Error = 8.19377 


The calibration equation is therefore, 


Response FactorH,O= 1.0468 ..... 
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Table A-1.9 


The coefficients of the polynomial, 


X measured 


0.75900 
1.20800 
1.80200 
2.78800 
3.62600 
6.39100 
9.43700 


12.15000 


X = mol % H,O 


Water Response Factor Variation 


Y = Response Factor, (H,O/N,)area / (H;O/N,)mol 


AO = 0.87473 
Al = 0.01918 
A2 = -0.00091 


Calibration Data 


Y observed 


0.88550 
@.922011 
0.82810 
0.87870 
108745 
SISTOSANE) 
0.94822 
099250 


Y calculated 


0.88876 
O:89657, 
0.90638 
0.92110 
0.93226 
C5907 
0.97430 
OS7278 


Variance = 0.005181 


Standard Deviation = 0.071976 


Maximum Pct error = 14.2471 


The calibration equation is therefore, 


PCT error 


0.36826 
2.84383 
9.44654 
4.82561 


14.24710 


BS957¢6 
2/0076 
ASS) 


Response Factor H,O = 0.87473 + 0.01918 (mol%) — 0.00031 (mol%)? 


V2 


: (A=1.8) 


a b=: 


APPENDIX 2 Sample Calculations and Analysis 


This Appendix begins with a definition of the names used to identify each of the 
variables necessary in the data-reduction. The experimental raw data from a single run is 
used to illustrate the application of the calibration data and the determination of the feed 
and reactor conditions. The establishment of the actual feed and reactor compositions is 
complicated by the fact that sulphur exists as a distribution of several different species. 
The derivation of the equations used to estimate the actual feed and reactor 
compositions will be given below and the iterative technique used to obtain the solution 
to these equations will be described in detail. The various kinetic parameters are then 
calculated. 

There then follows the calculation used to obtain the volume of the reactor and 
the estimation of the validity of the assumption of ideal (gradientless) Berty reactor 
operation as developed by Kuchcinski and Squires. 

Examples of the calculation of the rate constants are then given, followed by the 
calculation of the forward rate of reaction. Finally, the comparison of the two feasible 


sulphur species distributions is given. 


Sample Calculations 

The following calculation uses the experimental data obtained from the run 
WIO1A. The calculation follows the scheme of the data analysis program EXPTCAL, 
which is listed in Appendix 3. The variables are subject to calibration, whose derivations 


can be found in Appendix 1. 


In general, the variables can be identified by this key : 


A = area F = feed G = mass flow 

HO = water HS =H,S N=N, 

P = pressure R = reactor RT = R * Temperature 
S = sulphur SO SO, 


V = volumetric flow X = conversion Z = molar flow rate 
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The raw data for run WI01A were: 


Variable (name) Expt! value 
Weight of catalyst (WT) 0.6079 g 
D/P cell chart % (DP) 1.64 % 
Atmospheric Pressure (ATMP) 705.0 mm Hg 
Feed Pressure (FP) 8.35 Ibs/in’g 
Reactor Pressure (RP) 46.5 % 

Feed Temperature (TF) 34.0 °C 
Reactor Temperature (RT) 313.026 
Feed water area % (AWF) 2905 


Sulphur feed rate (ZS 1F) 


0.0014937 mol/min 


H,S area % (AHSP) 0.10306 
SO, area % (ASOP) 0.0883 
H,O area % (AWP) Ii22Ze 
Volume of reactor (VR) 664 cm: 


The first step in the data analysis was to apply the calibration equations to the raw 


data and obtain the actual experimental conditions. 


Convert temperatures to Kelvin ; 
Feed Temperature, TF = TF + 273.15 = 307.15K. 
Reactor temperature, TR = TR +273.15 = 591.15 K 


Convert pressures to Newtons/metre’ 
Atmospheric pressure, ATMP =ATMP * 101325/760 = 93922.3 N/m? 
Feed pressure, FP = (FP/14.17) * 101325 + ATMP = 149328.7 N/m? 


Apply calibration equation A- 1.4 to obtain the reactor pressure 
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Reactor Pressure, RP = (0.89957*RP-0.32023)*101325/76.0 + ATMP 
= 149264.1 N/m? 


Volumetric flow of nitrogen, apply calibration equation A-1.3 


VN = (0.14167 * DP - 0.096994) * 0.001 = 3.678924 * 10-4 m3/min 


Calculation of N, flow from the ideal gas law, 


molar flow N,, ZN = FP * VN/(8.314TF) = 2.182272 * 10-2 mol/min 


To calculate the mol percentages of H,S, SO, and H,O in the product stream ona 
sulphur free basis, the response factors (RF) from the calibration equations A-1.7 to 
A-1.9 are applied : 

RFHS = 1.92045 
RESO = 132023 
RFHO = 0.87473 + 0.01918 WP — 0.00091 WP? 


where WP = mol% water in G.C. product stream. 
Then, mol% H,S (HSP) =100 * RFHS * AHSP / DEN 
where, 
DEN = (RFHS * AHSP + RFSO * ASDP + RFHO * AWP +1.0 * 
(100-AHSP-ASDP-AWP) 
Similarly, mol% SO, (SDP) = 100 * RFSD * ASDP / DEN 
and, 


mol% H,O0 (AWP) = 100 * RFHO * AWP / DEN 


Since the response factor for water is a function of the mol% of water, the 
following iterative procedure was performed. 
il Assume the value of RFHO, the water response factor 


oa Calculate the mol% of H,O, H,S and SO). 
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3: From equation A-1.9, calculate the corresponding response factor. 

4. Compare this response factor with the previous one, if they differ by more than 
5 * 10°, repeat steps 2 to 4 until RFHO converges. 

Sf Stop. 


The mole percent of water in the feed (WF), is calculated from AWF by the same 


procedure, ie., 


mol% H,O feed = 100 * RFHO * AWF/ (RFHO * AWF * 1.0 * (100 - AWF) 


Once again iteration is required to determine RFHO. 


Following this procedure, the results that are obtained are : 


For the feed, 
RFHO =0.89576 
WEF = 1.15754 
For the product, 
RFHO = 0.893073 
WP = 1.00200 


HSPi— 0.197 155 
Sire 0;16572 
The mol flow of water in the feed is therefore, 
ZWE = WF * ZN / (100 -— WF) = 2.13182 * 10-4 mol/min 
The total mass flow rate of the feed ts, 


GFT = ZN * 28 + ZWF * 18 + ZSIF * 32 = 0.561620 g/min 


This completes the preliminary calculations on the raw experimental data. The 
calculation of the flow rates of the product stream components is left until after the 


distribution of sulphur species is calculated. 
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Calculation of Feed Composition 

Because sulphur vapour is an equilibrium distribution of several species, it is 
necessary to use an iterative technique to establish the actual mole fractions and mole 
numbers. The method used in this work was developed by Razzaghi (74 ). 

Only the sulphur species S,, S,, S; and S, have been included, since these species 
are most likely to be present and since there is no definitive data on the equilibrium 
composition of sulphur vapour. The data available on the even numbered species are the 
most reliable (74). 

Assuming the sulphur vapour to be at equilibrium, the Claus process consists of 
the following set of reactions : 


2s + SO; = 3/8'S;, 2,0 K1 


S,= 4S, K2 
S,= 25, K4 
S,= 4/3 S K6 


The Ki are the equilibrium constants of each reaction. These reactions form a consistent 
set of non-redundant reactions and are used to calculate the thermodynamic equilibrium 
of the Claus reaction. 

From thermodynamics, K, = exp(—A G°/RT), hence the the equilibrium constants of 
these reactions can be calculated from the Gibbs free energy data on the species in the 
reaction system. For any system at equilibrium, the Gibbs free energy is a minimum, and 
since the equilibria of the above reactions can also be calculated in terms of the partial 
pressures of the species, there is a system of equations which can be solved, and the 
equilibrium mole fractions of the species obtained. 

The mass balances on the total mixture and on the sulphur component form two 
equations in several unknowns. However, applying the constraints imposed by 
stoichiometry and thermodynamics, the number of unknowns can be reduced to only two. 
This is a pair of non-linear equations in two unknowns, therefore they may be solved by 
a suitable iterative technique. The derivation of these equations and the method of their 


solution employed in this work, is given below. 
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Assume unit mole of inlet gas, composition Yw, Yn), Ys, 
Let X be the conversion of H,O at any time. 
From amass balance, the following equations are derived : 
Yiw = Yw - X Yw 
Yn = Title 
Y'h,s = X Yw 


Y'so, = 2 X Yw 


Y's, = a; 
NS = 2 
¥'s, = ay 
Ys, = a; 
where (!) denotes products, and the a's are the unknown amounts of the sulphur species 


present in the product. Total moles (Nt), in the product, is therefore, 

he Gee a Alu ze ae oP Ga aie car: a a= age Page < Pn ape yee Apes (A-2. 1) 
Performing a mass balance on the sulphur species, 

Sin = 8 Ys, 

Sout X YW % YW + 2a, +4:a, + 6 a, 6a; 


BevS n= hey Wie & YW + 2a) Peay tO a, 125 Go eae eee (A-2.2) 


For ideal gases, Keg. = Kp =) Phi/) Pi 


where Pi = the partial pressure of species i. 


==> K2 = a,t » Nt —> a, = )K2 el si 

“a, ¥ PE Ei acre mes eee (A-2.3) 
Se A= a ve NE ——> 2d, - |) Ko eae 

“a, *P aes Beret ee Od ay ae (A-2.4) 
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substituting (A-2.3) (A-2.4) and (A—2.5) into (A-2.2), 


8 Ys, = 3/2 X Yw + 2 (K2 a, P%) 5 + 4 (K4 a, P)S + 6 (KE * a, PO?) 75 + Ba, 


Nt? Nt Nt®-33 
sie St ee eee ee (A-2.6) 
This equation is a function of two variables only, a, and Nt. Let, 
A = (K2 P3)0.25° 
B= (K4 P)0.5 
C = (K6 P2-33)9-75 
then substituting in A-2.5 and rearranging, the following function is obtained, 
F,(®D) = (3/2 X Yw Nt®75) + (2 A a,%5) + (4 B a,%5 Nt25) + (6 C a,°-75 Nt®5) 
pal Salas NS. INTC) ot) Cee eee eee ae {A=2:7) 
substituting similarly in (A-2. 1), 
Fitd)i= Ni — NEY w © Yn, +4200 Yw) aA 2, BapNt'} — Ca Nv 
ria, NGS 10). A 5 eee eee (A-2.8) 


==> two non-linear equations in two unknowns. 

The Newton-Raphson technique is used to obtain the solution. This method 
requires the Jacobian matrix, and hence the partial derivatives of the equations must be 
evaluated. 

These are found by straightforwardly differentiating equations A-2.6 and A-2.7 


above, viz. 


Jised ae ea te? Beare eNitia ea. Gage Ni + 8 Nt°.75 


ore 25-A wi gn 


' ot 
GD owe ae 
a9 6 3) + MOA hm She ttn oe x SNE IT 
iN SA Ooratteth ie 
s —_ 
. sale Mahia ambi 
Lana = 9k Cer 9 at ora will = 1 
(@A-Ai ons + 


y i) 
borisety ait inspsin ited 
ed taum Ae ame atic” 


a", 


TSrA ons OE-A erent whit 


120 


Of,/0 Nt = 1.125 X Yw Nt®25 + 0 + B a,%5 Nt-%75 + 3 .C a,%75 Nt-%5 


= O Nt? (a= 1's) 


ee rs ea tA At (A Zan) 
Oly Ody OO ae Ane BINT Ge CON tee so NICO: (A-2.1 1) 
Of,/A Nt = 1.75 Nt%75 — 0.75 Nt-25(Yw + Yn,+ % X Yw) - 0 - 
O72Z5) Branint iC ay Nits Oe OrdaiNcme 
Sy ener Mie yee (Aq2a1 2) 


To solve equations A-2.6 and A-2.7 for Nt and a,, the constants in these 
equations must be evaluated and so the equilibrium constants K1, K2, K4, K6, must first 
be calculated. These are obtained from the Gibbs free energy of the species. The Gibbs 
free energy of each species is calculated using the standard free energies and specific 
heat data as given by Rau(7 1), in the case of the sulphur species, and JANAF (36) for the 


others. The Gibbs free energy of each species is calculated from the relationship : 


Gi=-a (1,0 -l6g 1)" 1/2 bile el Gpc sie tr 1 
1/20 e T' +f/T - g - h/T? 


where T = temperature in Kelvin, and the constants a to h are the empirical constants of 
the specific heat polynomials as given by Rau and in JANAF. The derivation of the above 


relationship may be found in the CRC handbook (14) (page D-6 1). 


ole ime 
Gh,s = 3.9163 * (1 — log(593)) + 3.5138 * 10-*593/2 
a2 19 Sie 0S 5937/6 
+2 7454 4710593312 =4 85834593720 
-3609./6/593 - 2.366 - 0 
Gh,s = -29.6403 
Similarly, 
Gso, = -91.00975 
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Gh,o = -72.51437 
Gs,= -2.039625 


Gs, = -9.48614 
Gs, =42538013 
Gs7 = =S5:306 12 


The changes in the free energy are calculated for each reaction. For example 
ASh SS, 
AG =4" Gs, -— Gs; 
Thus, G= —2.039625 * 4 + 35.30612 = 27.14762 ( = DFRT2) 
Similarly, 
DFRT4 = 16.338 
DFRT6 = 1.46596 
DFRT81 = 7.97693 
The equilibrium constants are then calculated from the following equation, 
K2 = exp(-DFRT2/RT) 
nus A hens 62 1579154 a).074 
K4 = 8.0593723 * 10 * 
K6 = 0.2308563 
Keite=3.4329075 4107" 


The feed composition can now be calculated using multiple iterations. 

The mole fractions of hydrogen sulphide and sulphur sige the feed are zero. 
The mole fractions of nitrogen and water are estimated assuming the sulphur to be 
entirely S,. The total mole fraction is assumed to be 1.0. For the case of the feed, the 
conversion is 0.0. The assumption, as a first estimate, that the feed consists entirely of 
S; Bila the calculation of the feed molar flow rate, 


total feed moles, ZFT= ZN + ZWF + ZS1F/8.0 = 0.01860356 mol/min 
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The initial mole fractions can now be caiculated, 
YW = ZWF/ZFT = 0.01145918 
YS8 = ZS1F/8.0/ZFT = 0.01993638 
YN = ZN/ZFT = 0.9785044 


The subroutine NLSYST is called, having first determined the various constant 
inputs required, A, B, and C. The other constants D and E are required but both of these 
are functions of the water and nitrogen mole fractions and are subject to change as the 


correct mole fractions of S, is obtained. 


Operation of NLSYST 

The two functions and their partial derivatives are evaluated using the first guess 
estimates. An external subroutine LEQT 1F is called from the IMSL program library to 
solve for the roots of the functions. The subroutine obtains the solution of the set of 


linear equations, 


Matrix A contains the Jacobian matrix and hence this scheme follows the Newton 
Raphson iterative technique where, 


Kean TO fe OG) 


The system of equations is solved repeatedly until f, and f, = 0, to within 10°. 

The functions are re-evaluated and when the value of each is within a certain 
tolerance, the routine returns to the main program with the new estimates of Nt and Ys;. 

The mole fractions of every species are now calculated using equations (2.2) (2.3) 
and (2.4). : 

Mteor= Oey Sita Neo 
The partial pressures and mass flow rates of each component are calculated : 
PPi = mfi * Press(atm) * 10 1325(N/m?/atm) 


GFi = mfi *«MWi * ZFT, the total mole flow rate 
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The actual sulphur feed mass flow rate is calculated : 
SGFA = ZS1F * 32.0 g/min 
and the sulphur feed mass flow rate is calculated, 
SGF = SGFi 

The constants D and E contain the initial values of YN YW and YS8. and so the 
Newton Raphson technique must be re-applied using the new estimates of the mole 
fractions. Two tests are therefore performed to test for overall convergence. 
Ie The total mole fraction is evaluated. Since unit mole of feed has been assumed, 

the total mole fraction will converge to 1.0. The actual mole feed is adjusted 


according to the component mol fraction, viz. 
ZFT, = ZFT, * TFME (total mol fraction feed) 


The mole fractions are now recalculated using this new estimate of the total mole 
feed rate. 
Pa The mass balance on sulphur is evaluated in order to obtain a new estimate of 


Yoo, vize 
(YS8), = SGFA/SGF * (YS8), = actual mass flow S/ calculated mass flow 


The constants D and E can now be recalculated and the routine NLSYST is called 


repeatedly until both TFMF and SGFA/SGF = 1.0. 


Calculation of actual flows 
The feed composition and the relative amounts of the various sulphur species are 
obtained in this manner. 
The volumetric feed flow is calculated from the ideal gas law, 
VET = ZFT * RT/RP 
Before the rate of reaction can be determined, the distribution of sulphur species for the 
reactor conditions must be calculated. This is done in a similar manner to the feed 


composition, only in this case the conversion does not equal O. 
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ZN = as before 
ZWP = WP * ZN/(100 — WP) 
Therefore, conversion = (ZWF — ZWP)/ZWEF 
The initial guess to be fed to NLSYST is calculated. The total moles is again 
assumed to be 1.0. The mole fractions of nitrogen and water are calculated again 
assuming all the sulphur is in the S, form. 
Total product moles = ZN + ZWP+ZS1F/8.0 *(1.0 - 3/16 X) 
==> YN = ZN/ZP I) 
VWs ZEAE IN oo = 25 1h, 60" | 1093/1 Oey eed 
This initial estimate of the variables is passed to the iteration loop and the solution is 
obtained as in the case of the feed composition. 
The volumetric flow is calculated from the ideal gas law : 
VET = Zea Rae 
The atom balances on each of the elements are calculated, for Hydrogen for 
example, 
AHF = mf H,O feed * ZFT * 2.0 
AHP = mf H,O product * mf H,S product * ZPT* 2.0 
Evaluation of kinetic variables 


The various kinetic variables can now be evaluated, 


Residence time, =0.001 * Volume of the reactor/(volumetric flow prod. *60) s 


Space velocity -> sv = 1/t 
Conversion of water, as before, = (ZWF-ZWP)/ZWF 
Conversion of sulphur 8, = (MF S,*ZFT — MF S,p * ZPT) /(MF S,f * ZFT) 
Water rate of reaction, = (ZWF — ZWP)/WT mol/s g cat 
Sulphur rate of reaction = (MF S, * ZFT — MF S,p * ZPT)/ WT mol/s g cat. 
X - axis of initial rate plots = WFS = WT/(MF S, * ZFT * 60) 

WEW = WT/(ZWE * 60.0) g mol/hr 


The remainder of the program deals with printing the output. 


124 


voy orn ae tant 


oe 5 - a 7 7-rs —_—\ ae 
16g) vst we <stue kw in enatcett scone bf ad et beers 
> 7 a = : oo a va a 
erect. orth ni oh ae _ es 
a ' - ; 7 7 
sare -ohegeiararews - PL" 220" trains iy Nate an 
. - 
rai = AY aa eo oe 
i] 7 7 
*% sien vs ea plan ’ 
re Pe 2 sterey er ? @ g e166 a4 ago hg J oy %% armas | atl nn 


a 
& aw! vt tones ot es Rar 
) @ep Bh) ees SGA LIES wor a ai 


we 


fe eee - 


7 
cestterint wh inte cies qe werehn ns Ge 1S dane ip eexnled MRA? 


tra i> AA 


“7 


CT) mec Sw Huptke fn tn = FA 
evideney pirentt Yo nolieuis a 

Sets qury off wrens Mer YR ~Joree sonia ott a 

: 

ig meta’ © 198.0" ont woneiela@” | 7 


il ews ersotew onan 


TATA STWR eretedes “er to moveravnad: - 


e 


(CG es 4c}! o Tang e376 


(1.2 whe . Qe We TAte8 OH < S weine to NOS IRS 


ec 6 chee THING ~ THES = an? he ane a 
ANE 6 EN - TREw (2 Bee ees to et 

_ - . : a 
_ - 


aog (tor TW 
he ms ps a” vo ‘tala wn eis te 


125 


Berty Reactor Ideal Operation 


The set of equations derived by Kuchcinski and Squires were developed using a 
reactor similar to the one used in this work. They caution that the constants of these 
equations will be different for other Berty reactor configurations. However, since they 


used an Autoclave reactor it is probable that these equations are applicable to the reactor 


used in this study. 
For Reynold’s number < 50, the following equations apply, 

Z/3 Ye Ye 1/4 
Sle=s bir per Pr (T/M P) (1500/rpm) (4/(L/Dp) )/ 


[CO (a Oi 825 re ern ere anes ae (6.1) 


AX = 47000 r p pv (T/P)-16/ 


[la W )%5 (Dim)®-§7 (M/ )°46 (rpm/1500)5 (4 Dp/L)®.5] 


Slightly different equations are applied for Re > 50. For this reaction, 


M = 29 a= 6/D = 1890 
B—F1.5:48) 08 Pa Viet 

T= 600K M = 3% 10-5 Kg/ms 
rpm = 1200 L/Dp = 1 


rp =667* 10-5 kg mol/s 


<— 


=) 
(ay) = 1890 m? / m(bed) Dim = 1.3 * 10-4 
H = 10% J/kmol Pr = 0.68 
Cp = 1000 
Substituting in the equations ——-> 
Re = 16.8 
X = 4.87 * 10% 
a= 2. U0 102 16 
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Calculation of Reactor Volume 


The volume of the reactor was determined by experiment in the following 
manner. A large vessel of known volume was connected to the the sealed reactor via a 


valve. A pressure gauge was inserted to measure the pressure as shown in Figure A-2. 1. 


Figure A-2.1 Measurement of Reactor Volume 


The large vessel was pressurised using the building compressed air system to 
about 70 psig and connected to the reactor as shown. The pressures in the vessel and 
the reactor were recorded, the reactor pressure being assumed to be the same as 
atmospheric. The vaive was then opened slightly, to allow the gas to flow slowly into the 
reactor, and the pressure was allowed to equilibriate. This final pressure was then 
recorded. The following assumptions allow this process to be used to calculate the 
volume of the reactor. The air is assumed to behave as an ideal gas, and the process is 


isothermal. The volume of the reactor can be obtained viz. 


oF ema: peezerqua> 

eye, 

oF in 

stream att! ow emcee eis eh WARE 
gnibiiua 
pare 
| ae ncveragel ent 
er ni 


5 
sacar ect 


Ea 


Conditions for vessel 1, P,,V, and n, 
Conditions for the reactor, P, V, andn,. Let (1) denote final conditions. 
Applying the ideal gas law, 
P, V,/n, = P, V,/n, = RT 

now, 

nis ne inyaned 9 Via Ve VV, 

ee ea aa ined) aes Nein 
and, 

Ay = (Pa ey Va) 

ae PV NIM tee ee Vale, Val = Piva 
which on rearranging yields, 


V, = Vi(P? — P,P, — PY) 


The volume of vessel 1 was 500 in?, atmospheric pressure was 700 mmbg and 


temperature was 25 °C. 
Calculation of Rate Constants 


Equation 5.2 is the model which has been fitted to the initial rate data, 
b 


n 
—Rw 10 = a (Pw/10%) /(1 +c (Pw/10%) ) 


With n=2, the parameters a, b, and c have the values, 0.3358, 0.8911 and 0.2740 
respectively. | 
Thus, 

—Rw = (Pw/104)989 0.3358/ (10 (1+ 2.74 *10° * Pw)?) 


—Rw =9.1555 * 10-6 Pw / (1 + 2.74 * 10-5 Pw)? 


then, 
k= 9 1555 10-/(Ps,)*>’ 


me Sait seg ea > 
. - ~~ AMA - ae 


bre BHM GOT ver rae mu 


>" 6S saw srumeqmns 
- a 


atiisatm eee tiv tae ont i $2 rita 
( “pekywea =} 5 = OT alt 


/ORTEO yae)\ 1880 RACED amily em geet a kee <2 


Since Ps, = 753 N/m’, and the order with respect to sulphur is 0.39, 

k= 9 1555.%910-*/(753)? = 6.4606 #1077 
This will equal the least squares value of the rate constant since it is derived from the 
least squares fit to the model. 


The rate constant may be calculated directly from the each initial rate viz., 


k = —Rw (1 + K Pw)?/ (Pw) Ps,"9) 


1.904 * 10°71 + 2.74 * 10-* 8204)?/ (8204°9 7539-39) 


7O1S4.*. 107 


Comparison of Forward and Reverse Reaction Rates 


From Table 5.8, the initial rate of the reverse reaction at a partial pressure of 


8204 N/m? water and 753 N/m? sulphur, and a temperature of 592 K, is 


0.01904 mol/g hr. 


Using equation 2.2, the rate of the forward reaction is given as, 


-Rh,s = 1.121 exp(—7440/RT) (61.5) (5.6)%5/ ( 1 + 0.00423 61.5) 


SO Tasks, 


Thus the forward reaction is 9.75 times the rate of the reverse reaction. 
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Comparison of Sulphur Species Distributions 


At a temperature of 5S2 K, and a pressure of 1.46 atmospheres and a mole 
percentage of sulphur of 0.388, the distributions of the sulphur species that are 


predicted from the thermodynamic data are as follows, 


Species mol % (8) mol % (even) 
S; 0.0 

3) S03 5.74 

S; 0.44 

Sy 0.26 0.27 

Sr Eh Ele 

S< 36.44 60.06 

S; PEAS. 

Sy 35.4 Sas 


The distributions are very similar although that which includes only the even 
species has an average molecule size which is slightly smaller, 6.4 as opposed to 6.7 for 


the complete distribution. 
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APPENDIX 3 Computer Programs 


This Appendix contains listings of the computer programs used in the analysis of 
the experimental data| There are three main programs, EXPTCAL IRANAL and ADSREG, 
each of which uses a subroutine. The subroutine NLSYST is employed by EXPTCAL, and 
the subroutine NLREG is used by both IRANAL and ADSREG. 

The raw data analysis and kinetic data calculation were performed by the 
FORTRAN computer program EXPTCAL (for EXPerimenTal CALculations). The calculations 
performed by this program have been illustrated in Appendix 2. The basic outline of the 
program is shown by the flow diagram below. 

EXPTCAL employs a subroutine NLSYST (for Non-Linear SYSTem). This routine 
solves the non-linear pair of equations in Nt, the total mole fraction and Ys,, the sulphur 
mole fraction using the Newton-Raphson technique. A flow chart outlining the main 
sections of this subroutine is given. NLSYST calls a subroutine LEQT 1F from the IMSL 
(International Mathematical and Statistical Library). This routine solves a linear system of 
equations of the form, A X = B. The values of the kinetic variables of interest are 
calculated by EXPTCAL. A full listing of the output from this program can be found in 
Appendix 4. 

The main program IRANAL (for Initial Rate ANALysis), sets up the non-linear curve 
fitting subroutine NLREG (Non-Linear REGression). IRANAL calculates the form of the 
function being fitted and its derivatives. The important features of this program and its 
subroutine are elaborated on below. 

The main program ADSREG, (ADSorption REGression), fits the water rate function 
to the data using the subroutine NLREG. Several versions of this program were used 


depending on the power of the denominator term in the model. 
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Flow Chart of EXPTCAL 


Begin 


Read in expt! data 


Apply calibrations 


Calculate Sulphur reaction equilibria 


Calculate constant input to NLSYST 


Make initial guess of feed composition 


Calculate variable input to NLSYST 


Call NLSYST 


Re-estimate YS8, 


sulphur mole fraction 


Evaluate S mole fractions 


Evaluate Partial pressures and mass flows 


y 


Is S-feed = S mass flow calc? 


Evaluate all mole fractions 


Is the total mole fraction = 1.0? 


Calculate conversion 


evaluate first guess to NLSYST for products 


Calculate variable input to NLSYST 


Call NLSYST 


as above 


Calculate Atom balances 


Kinetic Calculations 


Print output 
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Flow diagram of NLSYST y 


Call from main 


Input from main program 


Evaluate functions 


Compare functions 


Tolerance met? RETURN 


Calculate derivatives 


Scale derivatives 


Call LEQT 1F 


Estimate X(1) and X(2) 


Tolerance met? 


RETURN 


soe 
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Experimental Calculations 
FER EE OK 25 AK OK KOK 2K OK 2K 2K ok 2K OK 2K 2K OK 2K OK 2K 2K 2K 2K OK ok OK KOK KOK KKK KKK KKK 


MAIN PROGRAM FOR EXPERIMENTAL ANALYSIS 
FEED ASSUMED TO CONSIST OF N2, H20 AND S8 
CEBU SpREACTUON@ Ss 2enaogt S02 <=> 27420 +. 376,585 31 


Ge S22 <=> 58 28 
2, ohe<=27 758 48 
A /3e Sb <=> O65 68 


INPUR@RG)THIS PROGRAM SHOULD BE IN THIS 
RORMAr: = 
NREPSs NOSORSCOPLTES OF OUTPUT REQUIRED 
IREP,; TERMINATION; CODE, NEGATIVE IF NO 
MORE DATA 
RUNDAT(I), STRING WHOSE FIRST FIVE MEMBERS 
ARE THE RUN NUMBER FOLLOWED BY THE 
DATE IN THE FORM 02/04/66. 
Wi wise@ineewelGhIOF CATALYST USED 
RONG oe THES LMPELUER SPEED 
DBeSmuGeeD/ POCELL SCHART READING =, 
ATMP IS THE ATMOSPHERIC PRESSURE IN MM.HG 
Gee ome EEDe PRESSURE SINGPSIG 
REC OMB EBREACTOR PRESSURE INSCHART, 
ieee Op tee ELD eEMEPERATURE MING CELSIUS 
Rev S WHET REACTOR TEMPERATURE IN GEUSTIUS 
iWECmUREeAREAS SSORSWATER SENS THESE EED 
Est bE eioeine GM=ATOM FLOW. OF SSULPHURS REED 
ABSPSeASOPy AWP ARE THE AREA® 4S. 0F 
Ze, o02, AND WATER RESPECRIVERY 
T NST GESPRODUCES 
THE DATA REQUIRED FOR THE CALCULATION OF 
TREY SULPHUR® SPECIES EQUILIBRIA AREMREAD 
FROM UNIT 7 WHICH IS ATTACHED TO THE 
GEN1, A LISTING OF WHICH IS GIVEN 
FOLLOWING THIS ONE 
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IMPLICIT REAL*8(A-H,O-Z,$) 

DIMENSION X(2),FRAC(7),F(2,1),1SPN(7),SNAM(7),FRE(7,8) 
#,FME(8) ,EMF(8),FPP(8),EPP(8),GF(8),GP(8) ,CHEM(8) ,WMOL (8) 
#,RUNDAT( 16) ,FRT(7) _ 

PAT EMCUEMGIN2 = day tH 2 Oinde, GH2 Spat SOCRMieS an oA eae 


&SOlRae Some e/ 
BATA WMOL/28.01,18.02,34.08,64.06,64.13, 128.26, 
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CONTINUE 


EXPERIMENTAL DATA INPUT 


A0=0.87473 
A1=0.01918 
A2=-0.00091 
AQ2=1.15950 
A12=0.00227 
A22=0.00106 
RFHS=1.92045 
REoD=Imo202¢ 
RFHS2=1.84419 
RESD2Z=1820770 
DAS amo 9, 
ITA=0 

ICT=5 

ICuU=4 

WREPEL R73) 
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PORN ATHMSeeSEONT=12000B0LUD. 9° FIXEDR LANDSCAPES Wve t = 


(Pato ak hWeDAtA, [ABLES (= sa 60k oRAWeDEIEA 
TABEES ///// 45x Wt. Cat:— Units are, Grammes 
//,45X,’RPM - impeller speed, 
//,45X,'D/P - differential pressure cell 

reading, Chart Percent’ , 

//,45X,’ATMP - atmospheric pressure, mmbhg’ , 

TPR a heeke - feed pressure, psig’ , 

PEO XeARP - reactor pressure transducer reading, 


‘Chart Percent’ , 


pia aA AE - feed nitrogen temperature, 
py Oey nae - reactor temperature, 
//,45X,’AWF - Area Percent water 
Wy PAD eA7 Sa © =emo lat tow 'S1 
//,45X,' AHSP - Area Percent H2S in product, 


sulphur free basis’ , 


//,45X,’ASDP - Area Percent SO2 in product, 


sulphur free basis’ , 


Revs per Minute’ , 


feed, Mol/Hr’ , 


/ 
, 
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//,45X,’ AWP  - Area Percent water product, 


‘free basis’ ,/) 
READ(5,12) NREP 
FORMAT ie 
READ (5 
IF (IREF 


REP 


pal 
a 0) GOTO 28 
R 


| 
6. 

Oe rime ana ed BR 
6M 2a v2 2h wey) 


PELs i+ i 
DPehCT NE 6 GOIO.3 
IC=6 

IPAGE=1 

ICT=0 

GOTGeS 

IC=0 

IPAGE=0 


WRODEES sh) IPAGE, IC, (RUNDAT( 1) b= 
ANPP O RP CRS AWE 25 lms Adore yaoDP, AW? 


FORMAT CI 7 Pht, 


TAK OMRON Seok eeUA TE AaX  WheCat.: | 
AED PEA Near i ME ee oka PP VAX, REPU OX, 
POR On (re Oman niet Uwe) wo, bO Hob O us ie areas 
PASM AN ENO Amel 1 EEA ANSP 
Pa oNe noe ae Ene REO! 5 YEO o> eR Oe 7) 


ENDSGE SINPUT 
FEED VARIABLE MANIPULATION 


VR=664.0 

LR 2hS so oer 
LRe2 hor ott 
RT=8.314*TR 

ATMP= ATMP*101325/760 
FP=FP*6895+ATMP 


RP=(0.89557*RP-0. 1383) *101325/76.0+ATMP 
VN=DSQRT(0.08998*DP-0.05350)*0.001 


ZN=FP*VN/8.314/TF 
DEGEREPYEOQ>2)G010 86 
RFHO=1.0 


WF =AWF*RFHO/ (RFHO*AWF+(100.0-AWF) )*100.0 


RFHOP=WF*(WF*A2+A1)+A0 


IF (DABS(RFHOP-RFHO).LT.0.00005) GOTO 43 


RFHO=RFHOP 
GOTO 44 


sulphur’ 


WPM DPS 


Bia 
5X,’ AWP! 


b] 
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CONTINUE 

GOTO 87 

RFHO=1.0 

WF =AWF*RFHO/ (RFHO*AWF+(100.0-AWF))*100.0 
RFHOP=WF* (WF*A22+A12)+A02 

IF (DABS(RFHOP-RFHO).LT.0.00005) GOTO 33 
RFHO=RFHOP 

GOTO 34 

CONTINUE 

ZWF=WF*ZN/( 100. -WF) 
GFT=ZN*28.01+ZWF*18.02+ZS1F*32.064 


CALCULATION OF FEED COMPOSITION 

Zi @lSsiGRPMGEE (PEOW*Cr =SUPPRHUR RIE FALE Si 
USEDSASSANSESUIMATE 10 FIND 

THE EQUILIBRIUM COMPOSITION 


PRESS=RP/ 1013825. 
T=TR 


CALCULATIONS INCLUDING THE CALIBRATIONS 
AREAE Geo ePORNOLE 4° S 


HeeUReP ee Once GOLO 8B 

RFHO=1.0 
DEN=RFHO*AWP+RFHS*AHSP+RFSD*ASDP+ 
1.0*(100.-AWP-AHSP-ASDP) 
WP=AWP*RFHO/DEN* 100.0 
SDP=ASDP*RFSD/DEN* 100.0 
HSP=AHSP*RFHS/DEN* 100.0 

RFHOP=WP* (WP*A2+A1)+A0 

REUDABS REHOP-REHO) 17 .0.00005) GO10 46 
RFHO=RFHOP 

GOTO 45 

CONTINUE 

GOTO 89 

RFHO=1.0 
DEN=RFHO*AWP+RFHS 1*AHSP+RF SD 1*ASDP+ 
as er00 AWPSAHSP-ASDP) 
WP=AWP*RFHO/DEN* 100.0 
SDP=ASDP*RFSD1/DEN*100.0 
HSP=AHSP*RFHS 1 /DEN* 100.0 

RFHOP=WP* (WP*A22+A12)+A02 

IF (DABS(RFHOP-RFHO).LT.0.00005) GOTO 89 
RFHO=RFHOP 

GOTOG3S 
CONTINUE 
ZWP=WP*ZN/(100-WP ) 
ZHSP=HSP*ZN/ ( 100-HSP) 
ZSDP=SDP*ZN/(100-SDP) 
ZWF =ZWP+ZHSP 
ZWF=WF*ZN/ (100. -WF ) 
GFT=ZN*28.01+ZWF*18.02+ZS1F*32.064 
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C 
: CALCULATION OF REACTION EQUILIBRIUM CONSTANTS 
DOg20 ais 7 
FRT(U)=FRE(U,1)*(1.-DLOG(T))-FRE(uJ,2)* 
#7/2.-FRE(J,3)*T*T/6. 
#-FRE(U,4)*(T#*3)/12.-FRE(d, 5)*(T**4)/20. 
See une fie rREt Ue ERE 61/171 
20 CONTINUE 
Dik Pel 2eeG Rie iE R121 5*FR1 (4) 27 0tER To) 
DFRT41=2.0*FRT(1)+FRT(2)-0. 75*ERT(5) <2 *FRT (3) 
DRAG d= 25 O4ERI) FRI (2)-0.5*FRT(6)-2 0+ ee 
DPRDGi=2 SOFER TTI tERTG2) SONS 754ERM( 7h-2 ERT (3) 
RK2T=DEAPG=DFRI21) 
RK41=DEXP(-DFRT41) 
RK61=DEXP(-DFRT61) 
RK81=DEXP(-DFRT81) 
DFRT2=4.0*FRT(4)-FRT(7) 
DERE4=2¥0*FR1(5)-FRI(7)} 
DFRT6=4.0/3.0*FRT(6)-FRT(7) 
RK2=DEXP(-DFRT2) 
RKA=DEXP ( -DFRT4 | 
RK6=DEXP(-DFRT6) 
C 
c CALCULATION OF FEED COMPOSITION 
C 
C 
C ASSUMPTION THAT FEED IS S8, (FIRST APPROXIMATION) 
C 


16 


YHS=0.0 

YSO=0.0 

YW=ZWF/ZFT 

YSS=2515/8.07Z5 1 

YN=ZN/ZFT 

W=1.0 

PRESS**3*RK2)**0.25 
RESS*RK4) **0.5 

RESS*t4 1073.0) *RKG) +0..75 


IW=1,50 
LRuiWesQgo0 heGOLO R32 

D=YN+YW+0.5*CONV*YW 
=1.5*CONV*YW-8.0*YS8 

CRIRINGS STuey ARB CD. Ey pW) 

FMF(8)=X(1) 

FMF (5)=(RK2*(PRESS/X(2 


(2)) 
FMF (6) =(RK4*PRESS/X(2)*X(1))**0.5 
FMF (7) =(RK6*(PRESS/X(2) )**(1.0/3.0)*X 
TFMF=YN+YWtFMF (5)+FMF(6)+FMF(7)+FMF(8 


FMF(1)=YN 


**3%*X(1))**0.25 
(1 


eso) es 
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ou iui 
(e5e> T< 
(DS) 


eco an) ioZor 
MOL (1)*ZFT 


a 
= * 


CONTINUE 

SGF=0.0 

DO 18 Pao 8 

SGF =SGF+GF (1) 

SGFA=ZS1F*32.064 
PReUVabStsGRAE SGE) Et 20.21D-08) GOO 19 
DS=SGFA/SGF 

X(1)=XIP*DS 

YS8=X(1) 

XIP=X(1) 

GOTR@e21 

Zl = Zee 

YN=ZN/ZFT 

YW=ZWF/ZFT 

YS8=FMF (8) 

X(1)=XIP 

X(2)=1.0 

PE OARS TEMES mOU00} et OeiDi- CS). GOTUs22 
CONTINUE 

VFT=ZFT*RT/RP 


CALCULATION OF REACTOR CONVERSION 


CONV=(ZWF-ZWP ) /ZWF 

X(2)=1.00 

RIPEAU4 } 

XW=CONV 
ZPT=ZN+ZWP+(1.0-0.1875*XW)*ZS1F/8.0 
YN=ZN/ZPT 

YW=ZWP/ZN 

EMP GSy=ZHSP/ZPT 

EMF (4)=ZSDP/ZPT 

Woo= (ie OMe B75 xW) ZS 7/8. 07 ZPT 
DOR25 IW=1,50 

PeeWee Qeo0) eGOTOe 32 
D=YN+YW+0.5*CONV*YW 
E=1.5*CONV*YW-8.0*YS8 

CAINMMRNESY SIUX¥AYBSCUD,E, PW} 


EMF(8)=X(1) 

EMF (5) =(RK2*(PRESS/X(2) )**3*X(1))**0.25 

EMF (6)=(RK4*PRESS/X(2)*X(1))**0.5 

EME (7) =(RK6*(PRESS/X(2))**(1.0/3.0)*X(1))**0.75 
EMF(1)=YN 

EMF (2)=YW 
TEMF=YN+YW+EME (3) +EMF(4)+EMF(5)+EMF(6)+EMF (7) +EMF (8) 
GPT=0.0 
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28) 


24 


IGG, 


DBD) 


Dore I=1,8 

EPP(I - EMF(1)*PRESS*101325, 
GP(1)=EMF(1I)*WMOL(1)*ZPT 
GPT=GPT+GP(1) 

CONTINUE 


SGP=(EMF(3)+EMF(4)+2.0*EMF(5)+4.0*EMF 
&8.O*EMF (8) )*32.064*ZPT 

IF (DABS(SGFA-SGP).LT.0.1D-08) GOTO 24 

DS=SGFA/SGP 

X(1)=XIP*DS 

YS8=X(1) 

XIP=X(1) 

GOLOP25 

ZPT=ZPT*TEME 

YN=ZN/ZPT 

YW=ZWP/ZPT 

EMF (3)=ZHSP/ZPT 

EMF (4)=ZSDP/ZPT 

YS8=EMF (8) 

X(1)=XIP 

X(2)=1.0 

PSROABSViEMES1.0C0) ar. ORD 20a GOTUy2 

CONTINUE 


ATOM BALANCES 


AHF=FMF(2)*ZFT*2.0 

AOF=FMF (2)*ZFT 

ASF=ZS1F 

ATF=AHF+AQF+ASF 
AHP=(EMF(3)*2.0+EMF(2)*2.0)*ZPT 
AOP=(EMF (4)*2.0+EMF (2) )*ZPT 


(6) +6 702EN Ey) 


6 


ASP=(EMF(3)+EMF(4)+2.0*EMF(5)+4.0*EMF(6)+ 


& 6.0*EMF(7)+8.0*EMF (8) )*ZPT 


ATP=AHP+AQP+ASP 

VF T=ZFT*RT/RP/60.0 
VPT=ZPT*RT/RP/60.0 
Ghk=Gri7 S040 
GPT=GPT/60.0 
AHD=(AHP-AHF ) /AHF*100.0 
AOD=(AOP-AOF) /AOP* 100.0 
ASD=VASPASP )/ASP¥100. 0 


KINETIC CALCULATIONS 
TAU=0.001*VR/VPT/1000. 


SV=1.0/TAU 
XW= (ZWF-ZWP) /ZWF 


XS8=(FMF(8)*ZFT-EMF(8)*ZPT)/(FMF(8)*ZFT ) 
XS6=(FMF(7)*ZFT-EMF(7)*ZPT)/ (FMF (7) *ZFT) 
XS4= (FMF (6) *ZFT-EMF(6)*ZPT)/(FMF(6)*ZFT) 
XS2=(FMF(5)*ZFT-EMF(5)*ZPT)/(FMF(5)*ZFT) 
RATEW=(ZWF-ZWP)/WT/60.0 
RATES8=(FMF(8)*ZFT-EMF(8)*ZPT)/WT/60.0 
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RATES6=(FMF(7)*ZFT-EMF(7)*ZPT)/W1T/70.0 
RATES4=(FMF(6)*ZFT-EMF(6)*ZPT)/WT/60.0 
RATES2=(FMF(5)*ZFT-EMF(5)*ZPT)/WT/50.0 
RATEWS=RATES8*16.0/3.0- 

*+RATES6*4.0+RATES4*8.0/3.0+RATES2*4.0/3.0 
RATED=(RATEW-RATEWS)/RATEW * 100.0 
WFS8=WT/(FMF(8)*ZFT*60.0) 
WFS6=WT/(FMF(7)*ZFT*60.0) 
WFS4=WT/ (FMF (6)*ZFT*60.0) 
WFS2=WT/(FMF(5)*ZFT*60.0) 
WFW=WT/(60.0*ZWE ) 


C 

C OUTPUT 

C 

C 

C BUBEMGALCUBATED) VARIABLES 

C 
DOPS tea 1 NREP 
Whe (6, 90) (RUNDAT G1), D=2%6) RURENDAT GIs b= 715) 

90 PORMATVIRIG/ 74/7/77, 19x 7 RUN ] COA de Oh ee Mier eke One 
WRITE(6,91) WT 

Ca BORMAI M1 Ok; @COMARYST. KA S201 4X e024 e Gs) 
Wile vO, 92) PRESS) 1, REM 


ve BORMA MW) 7 19X. PRESSURE CATM)) 99 Fd255) 10x. 
(UGNMEERATUREMIK) Soe ri2n2 0k, PANS SPEED == eRe Oe heen eae 
WRIMEV6y2/ 0) RATEW; RATES, 
tRAMESG, RATESS RATES2, RATEWS, RATED 
27 BORAT S42R, WATER’ =10A,2 Se 12K0~ Sb ei) c 4a 
GIZA, soem) 
G/F AIK RAKE (MOL/S GAGA ve Ol 2k, El2 eo eioe 
GORAQEM (WS) @=) 9, EA2°5, 


&’ MOL/S.GM CAT! , 6 Sin '% DIFFERENCE = OF Ore 
WRITE(6,94) TAU, SV, PAT Ta eM boa Gl oi welenia (Mel wre dl 
94 FORMAT(/,19X,’ RESIDENCE TIME = cee 
DO rt) iA POPACES VELOCITY <> EI 2u4ay eS) 
iON OMOFESEEGW RPEEDPE6’ (E12, 454%, PRODUCTE Be ie74: 
&' MOO ee ae a, SMASS SOE COWS, 
ZOREED YS) FEL 2R4SX | PRODUCT SENQn a) CuSenes, 
SloAy AVOLe eh LOW FEED = PE 2e4 5X "PRODUCT / aan 
APOMEAS* A) 
WRITE(6,97) 


O7PmBFORNAT C46, 9B-FEEDY, 40%, “REACTOR a//, 19x, SPECIES 1% 
&’ MOLE FRACTION’ , 
OO) ie NRG AIS PRESSY, OX, MASS FLOW 
C7 Ae NOUMERERAG) TON@, 1X, PARTIAL PRESS 72%, 
AENASSEE HOW, /, 46%." (PASCAL) 7A, 1G/S15 724k | PASCALS 
DOXe CS) sh} 


DOSS TS 178 
Gawi=GEVE) 260.0 
GPii7=GR( Be) 760). 0 
WRITE(6,96) CHEM(I),FMF(1),FPP(1),GF(1),EMF(1I),EPP(1) 
1, GPT) 
( 


She CONTINUE 
96 FORMAT (4 
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99 
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AU 


80 
81 


85 


ee Ete s ON Eile wmore A e125.) 
WRITE(6,98) XW,XS8,XS6,XS4,XS2 
LORMAN (OCR WAVER ES IMXOASR? ORM SG on DKS Are 

Sek ole, 

&/,19X,’ CONVERSION pees Ey Gh) 
WRITE(6,99) WFW,WFS8,WFS6,WFS4,WFS2 
PORMAT LOX WimCAl/MOE FEED 8) 5(9,0b ao. Say) 

WRITE(6,37)AHF,AHP,AHD,AOF,AOP,AOD,ASF,ASP,ASD 


FORMAT(19X,'’ ATOM BALANCES PEED 0X7 ee RODUC aeeoe 


Sa eUUEr ERENCE ay, 
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et ORUGEN Eo yon 2X E lon oy or Onan) meter nO Ge Ne ae 
One et tO OA Oe 19h) SULPHUR abi Ome Ae ca ono 


Coe mote) 
WRITE(6,29 } RK21,RK41,RK61,RK81 


BORMAU 7 pa OX, SEQUIEIBRIUMPCONSTANES 7 7/, 24%. RKOT= | 


SDS mOmowy RAAT Se 
HO Uoe Oh AON aR KOmi= tO. 6, 544° RKSI=" 013.6) 
WRETE(G2 C0) RK2 RK4> RKG 
RORW AM eo) eeeRK2= 043.69 oks RAG= | Dts beaks 
PRKo= oo; Dlr Okey ¥ 
[CUSICUT1 
Dae CUSNESon GOTO 1 
IC=6 
IPAGE=1 
ICU=0 
GOTO 9 
IC=0 
IPAGE=0 
TEVRUNDAM Ri) EO2S) GOI0 82 
FA=WFW 
XA=XW 
GOTO 83 
FA=WFS8 
XA=KS8 
WRITE ( 
FORMAT 


LPAGE? 1C 

iA oe eRUN SS Ax, DATE TON eRe ome 
WT GAte | 3X5! Saver 

(RUNDAT (1), i. Hes PRESS,T,WT,SV,FA,XA, 
ae CEM Ei) l= 1, ee 

SAO cs sh ae 

Ox WATER® , We 
gt AG Oo! 


1 
8 
OF 


1 
1 
; AS ee oe Ore oon Omron 


eke 0 
8)*100.0 
TTASITAY 
CEC LIAS Gio GOTGes5 
WRITE(10,80) (RUNDAT(I) ,1=1,4) 
FORMAT(1H1,20X,‘DATA SET ‘ ,8A1,//) 
WRITE(10,81) 
FORMAT(17X, "RUN! ,2X,’PRESSURE(ATM)’ ,2X, 
&’ TEMPERATURE (K an we 'WATER(%) FEED’ ,2X, 


&’ SULPHUR (%) aa G.HR/MOL FEED’ , 4X,’ CONVERSION’ 


WRITE(10,84) TRUNDAT (1). 1=1,5) ,PRESS, 1 ,FMF (2) , 


ROA WEA CONV’ ) 
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&FMF (8), 
&FA,XA 
PORMAT AVE 2OX SSA TPE 1196,6X,F7.2,6x, 
ce OGOr ANE 10 CeO Et 386), F13.6) 
CONTINUE 
PEIIRERG ROU EGOT OR) 3 
WRETECOY 359) 


FORMAT | ae 20 mere (7, 0X,’ FEED CONCENTRATION NGT 
& CONVERGED’ ,//,20('*' )) 

SHOP 

END 


SUGROURINESNUSY STUN, ABYC, ORE. Fs Wy 
THIS SUBROUTINE SOLVES A SYSTEM OF TWO NONLINEAR 
EQS 4B MBNEWLONS 4S IMEIBODe 


IMPLICIT REAL*8(A- 
DIMENSION FSAVE(2 
WKAREA( 12) , XSAVE ( 
FORMAT (5D14.5) 
M=1 

N=2 

TA=2 

IDGT=0 

MAXIT=50 

polar Dao 
XTOL=.1D-08 

DO 48 I=1,MAXIT 


nh dO ber A nd 
Nanaia aon 
2 


XSAVE (u) =X(u) 


ZT4R=X(2)**0.25 

et ent aioe O*A*A4R/(ZT4R**3)+4.0*B*A4R/(ZT4R*ZT4R) + 
&6. O*C*AAR**3/ZT4R+8.0*X(1) 

F(2,1)=A*A4R/ (ZT4R**3)+B*A4R*A4R/ (ZT4R*ZT4R) 
&+C*A4R**3/ZT4R +X(1)4+D-X(2) 

PTEST=0 

Dior 28) weal! 

PCO ARS EGU eGie et OL) aEPES Foie iced 
FSAVE(JU)=F (ud, 1) 

CU ac EM Ul) 
CONTINUE 
IF(ITEST.NE.0O) GO TO 30 

REN 

F(1,1)=0.5*A/(A4R*ZTA4R) **3+2.0*B/ (A4R*ZTA4R) **2+ 

aa ST aera EEe 

DF(1,2)=(-1.0/X(2) )*(1.5*A*A4R/ZT4R¥*3+ 
OPE. TR aT aa 5*C*A4R**3/ZT4R) 

DF(2, 1)=0.25*A/(ZT4R*A4R) ¥*3+0.5* 
&B/(ZTAR*AAR ) **2+0.75*C/(ZT4R*A4R)+1.0 
DF(2,2)=(-1.0/X(2))*(0.75*A*A4R/ZTA4R¥*3+ 
lee. SHB*AAR*AGR/ZTAR**2+0, 95*Cx* (A4R**3)/ZT4R)-1.0 
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, IDGT ,WKAREA, IER) 


Wh sae [Wie WE fen) [ie a [Bw] (ual 
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GOT 33 

GOTO 47 
yD GIeV Ola) 

GO TO 49 


ISU plese Aste) (86 


CONTINUE 
48 CONTINUE 

RETURN 

RETURN 


Celie ORME T PRINT. LHe RESUET 
END 


50 
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Operation of IRANAL and NLREG 


Having obtained conversion versus mole flow rate data it remains to fit the data 
with a suitable function to obtain the initial rate from the slope at the origin. This slope 
may be obtained in several ways. A smooth curve can be drawn through the plot and the 
slope estimated by eye. A spline curve can be fitted to the data. A non-linear function 
can be fitted. This latter approach was used since it was consistent and also made use of 
all the data) The eye-fit would not be consistent, and the spline curve places a great deal 
of emphasis on the point closest the origin. All these methods were applied to the data, 
and the nonlinear curve fitting technique proved the most reliable. However, each 
method gave similar results, indicating that the use of each technique was valid. 

The curve, A = tanh(B x), was used to fit the experimental data. This curve has the 
advantage that at the origin the slope is equal to A*B, the two parameters which are 
estimated by the least squares technique. 

y = A tanhiB x) 
dy/dx = A sech?(B x) B 


Hence at y = x = 0, dy/dx = AB, since, sech?(0) = 1.0 


The least squares estimates of the parameters was obtained using the modified 
damped least squares (MDLS) non-linear curve fitting technique developed by Meyer and 
Roth (58). This method is illustrated by Wolfe (87). 

The program IRANAL listed below, reads the experimental data and calls the 
subroutine NLREG, which returns with the values of the parameters A and B, as estimated 
using the MDLS method. A flow chart of this method can be found in the references 
cited above. This program requires two subroutines MODEL and MDER. MODEL 


calculates the value of the curve being fitted to the data. 
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MDER evaluates the derivatives with respect to the parameters, viz. 


y =A tanh(B x) 
dy/dA = tanh(B x) 
dy/dB = A sech?{B x) x 


The MDLS algorithm is based on the Gauss~Newton procedure. It is reported by 
its authors to be a reliable and efficient method of solving a general non-linear least 
squares problem. In this work, the method proved to be fast and reliable. The fitting on 
the Tanh curves had no difficulties. The fitting of the water function was slightly more 
difficult since the variables required scaling in order that a solution could be obtained. 
This scaling of the data has been described in Appendix 2. 

The program ADSREG is included for completeness, although it is basically the 
same as IRANAL. The difference lies in the form of the function that is of interest and of 


course its derivatives. 
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EVALUATION OF THE FUNCTION 
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EVALUATION OF THE DERIVATIVES WRT THE PARAMETERS 


SUBROUTINE MDER(N1 
DOUBEE ERECISION A 
D(1)=(DEXP(2.0*X(2 
#(DEXP(2.0*X(2)*T | 
DUMI=2. eee 
FDEXPALME see SE. 2 in 
Di eAlerts 1) 
RETURN 
END 


148 


149 


aon feel clon cet eee ced cell ede aad alll andl nell nll ool 
-—_—=— SS Se eS Se 


ADSREG 


Initial Rate analysis 


MAINLINE PROGRAM FOR NON-LINEAR CURVE FITTING 


Y = C#W**N/ (1+K*W) **2.0 
PORS INIT TAL RATESDAT A ANALYSIS 


ID NEGATIVE IMPLIES NO MORE DATA 
FOR ANALYSIS 
LERNER EE SUNPUT INES ARE@ THE 
SReC UBECATTONSS@OF Site DAT ANSE 
NeeONUMBERSOR DATA POINTS 
Me= THE NUMBER OF PARAMETERS BEING ESTIMATED 
Peo unEeNUNBERS OP SKEPEATSDATASPOINES 
ALAM AND GAMMA ARE VARIABLES GOVERNING THE 
KE OReAr PROACH sl Osi Tea SOUT ON Paes En 
ALAM TO 0 AND GAMMA TO 10**(20) AS FIRST 
GUESS. 


X Gls) (ARERSRHESVINIT TAL (GUESSES (OF SIHESPARAMETERS 
Y(I) ARE THE DEPENDENT VARIABLE VALUES 

Tg ARE LAE INDEPENDENT AVARTABUE VARUES 
Wi ARE }UEAWELGH LING FACTORS "FOR EACH 


DATA POINT, NORMALLY EQUAL 10 1 


DIMENSION DAT(5) 
DOUBLE PRECISION X(3),Y(1 
DATA FSTAT,TSTAT,EXPVAR / 
READ(4,80) ID 

FORMAT (12) 

LEG Ovian sO PeSTOR 

READ (APO) DAT) l=ie5) 
FORMAT(5A1) 

READ(4,40) TEMP,ET,ATM, EAT 
FORMAT(F6. 1,56. ce DF7.4) 
SXP_ESXP,WP, EWP 
Ded) 

,M,L,ALAM, GAMMA 
Be? eas 2 | 
( X 

) 


55 om 
i 


N 

ap 
1) eel ar Sh 

2 


' & 


a a eye) 


ns 5 

ow 6 4 == oe ee ee “ =. = She TD 

~ vs 3 

‘Ti% Svar SOR 4 Ae KAREORD 4 an AN 4 

- ae we ‘et soones #7 he 

Vibe 933 is Epa, 27AR SAT TAL 7 3 

ed = aie atk =e Se ows Fa FAS > = she ne Rae as - 

(iq Xo Ch APLAR TeETADIM 
ajeyianA 704 
Ja en 8d FLERE Bernt TA3n 
R\as4¢ 247% 3 ASIOL TAD: 4193992 


iselow-STAG 30 S3EaNLKN = 
337 sie Toe eon PRT SAMAEAS SD BaeRUM “ SH ~M > 

At | Tan TARRAR 40 RSAMUM dT = J. 
» SVIEALRAY Bre AMMLD at Ae BALA 


\* ‘10? 2% rn oa "peveds 30 3749 
_ ifje0ed| OF een ah6. © OF: Ate 7 
RIE or 
1637 $A544: SUT * 10358 Antti HT am iT. 
: ay 21ae (Pay 74a PT CET 347 aaa (Liv: 
4799 SPOTIA4 Piao Siew 3HT ana Lie 
i 6? Jauge +) smut -THID9 ATAD © 


omt¢ Ge ¢e~oee Uee 4428 oS pee emedaseres oo cena eee owes 
a) Tag MOTeMgME 


19389 38 
FAiS4. Tad 


we ee &t 


Gute rt 
Were titagy 


Ad 
Tmt ; 

. 418) 

Ww ib - 


his t 


“13.08 1 
+“ ? 
E ab ue 


150 


(EOD ae O ee ome las Wl) 
SUs CE ORMAPUS Rt 2809 
WRastECOy / Ope DAs =10,5) 1 EMP Rede Mae 
&SAP,ESXP, WP, EWP 
70 FORMAT(//," DATA Se lea Oh ey 7, PEMBER AU Ric meee Oma 
* / / 
GEDSSty e *PRESSURE Slain ic Se Bad) Wc Pee ied aah, 
8’ SULPHUR or 
Ch 2. ee Le eee De be (eb WATER bal tal eh cael Shen idea 
Ger omy | 
WRITE(5,90) 
90 FORMAT(//,10X,'Y = Ca#We*N/((1+K*W)**2.0)’ ,//) 
CALL NLREGIN rN, BRK Gh G W ALAM. GAMMA, FSTAT,TSTAT, EXPVAR) 
Ra loe Gl O)@eGG Oat 
Sy iol 
END 
C 
C EVALUATION OF THE FUNCTION 
C 
SUBROUTINES MODER( NISMS NING Uy Xo ts 1s VARUED 
DOUBLEPRECISION X(N1), T(M, L) 
VALUE=(X(1)*T(1 Tyee DVN / C1 04x (3) +i epee) 2 0) 
RETURN 
END 
C 
GD EVALUATION®OF THE DERIVATIVES WRT THESPARAMETERS 
C 
SUBROUTINE MDER(N1,M1,N,M,L,X,1,1,D) 
DOUBLE PRECISION X(N1),D(N1),T(M,L) 
Dei p=ht 1s ee ee ae QEXNG) ST +4250) 
DUMI=X (1 ) T(1, 1) *#X(2) 
D(2)= DUMI*DLOG(T(1,1))/((1.0+X(3)*T(1,1))**2.0) 
DiesepUM! = (S2M0) F1ll, 1) / Vi eOe Xie) 1A leg 320) 
RETURN 


END 


a '- 

a aa | 
“' 2 re 
‘Vie, 19, 

\mayans, Tater. Tales 


om gH) WO WOTTAD. 


Pred Pe ‘Ts ql ber 
)|. gee Li ht boy AE tay aa oe = 


saarsana’ gat Tie “age wr 19 ore ; 


sade 7 
ae ee ie 33 

i/o, Beals cath Revit se aor _ 
‘tite + 


ie 


ieeelt ite mee ee 


G 
G * * 
C * NLREG * 
C * * 
G ** non-linear regression + 
C 
ee pe a at ee ee ee, ae ee ee ee 
C 
S SUBPROGRAM FOR NONLINEAR LEAST-SQUARES REGRESSION 
: MetRURP OR IMEYER AND ROTH PRESENTED! IN@SECT UNM (a /7e0F 
C WOLFE,M.A.,’ NUMERICAL METHODS FOR UNCONSTRAINED 
: OPTIMIZATION’ , VAN NOSTRAND REINHOLD, (1978). 
C R.B.NEWELL YUL alee 
C 
G Saat ecae Me Wess sas Most Maat Tan) ees tes dees tacts sae pi cs rece! Laas ies eae ce ae, San 10a fais eal tea ae aes fa eat, = acs al ms Lew ee pu eros Mae ee ee ee Se Dee 
C 
SUBROUTINE NLREG (N,M,L,X,Y,7,W,ALAM, GAMMA, 
*FSTAT, TSTAT,EXPVAR) 
INDVEGERSGTES TMI TEMP (10) SUG EMP (10) 
DOUBLESRRECTSIONEXUN|s oy T(M,L) ,W(M) 
ae), ee be Fe Se GU OWPR (OS 1OlexsvEUiO) 
*,Q(10,10),P(10), TEMPN( 10), TEMPNN( 10, 10) 
*, ALAM, GAMMA, F STAT, TSTAT, EXPVAR 
SP lala: Wiel SENSE Sens) | 9 Blanes eles’ 
PO POM pol Pe GNAAI 1 Ue, DABS, OSORT 
Da mN oN AM eer 05 2557) 
DAIAUBETAMEPSOSEPSINERS2UERSCa) BRR D>0D 9. 50-065. 50> 0b 
CESD- 067 1, D-04 / 
LT=6 
WRITER (W890) LT 
SOQeeRORMAT (2 TERMINAL’ , 13, “= ENTER LUN | FOR@EURL OGULPUT 
LU=7 
DOT ISSIK=E1gN 
3 XSVE(IK)=X(IK) 
NF=Q 
ITER=0 


C 
CeeUN TAPRVAERGESS OR Fh AND’ S 
C 


T=] 
8000 CONTINUE 
aula =X (1 
T=I+1 
Cols Gire NGO 10. 8003 
GO TQ 8000 


8001 CONTINUE 
SRG EANOSMUNGON A ONS M yb gkethtey geley Watt al sles lenin ta 


START OF A MAIN ITERATION 


QiG 


wales sRbIs 
40g) 5.t Se wee es 


| See Baa te aah 


ret aly 


is 
,5b0°Gé aoe. so-aih ‘ g2e7..S2%5. 


ph 
pt 06.1) 371fe)) 
; *yeyve 4u9 9? wu gata - §\ép ay "4 ik : : 
ty 


“,f=ai biog - 
dint List At bVEI eh ; 
Oa : 
O=83T Db» 


2OMs § WO 2a05A¥ JAITIML! 


Lm, 08, 0078 0 


J. -* _ 


IGE) 


8003 
8004 


8005 


302 


904 


906 


8006 


8008 


8009 


8007 
8010 
8012 


8014 


Toe 


ITER=ITER+1 


it 
Gr 
FORMAT 
seo EEL 
T=] 


ee 
WT; 


GT.0) 
NE.LU) 
( / 


WRITE (LU,900) ITER 
WRIDES (ET O0C)}s ITER 
GRR RK MAIN ITERATION’ ,15) 


CONTINUE 
X(I)=XTIL(1) 


T=I+1 


Leer ac) oN) GUs1U. 38003 


GO REG 


8002 


CONTINUE 


T=] 


CONTINUE 
BUY/SELL) 


Dede 
GOS10 


QL GiEM) GOSLET 8005 


8004 


CONTINUE 


SE 
CECE 
FORMAT 
iM 
ie Haley 
FORMAT 


SE UP 


as 
FORMAT 


rau. 
oNESEU) 


CU 
SNe ew tele 


FED: 


Giacje WRITE 
WRITE 
(oD.15 >) 


GT.0) WRITE 


eu ea eee alae 
bbl 2029 =A ASL} 


~ et 


(LU,904) § 

(on 04S 
(eeSUMSOP=SQUARES@ 1S! D2005) 
THE LINEAR EQUATIONS 


GT.O) WRITE (LU,906) ALAM 
(PADAMOIS ab d5.. 4) 


GAL ERAXS UNG MIS NGM, LAT Lot AST EMPN) 


T=1 


CONTINUE 


GUE) =0. 


K= 1 


DO 


CONTINUE | 
G(1)=G(1)+A(K, 1) *F(K)*W(K) 


K=K+1 


LE heGheM GOs a? es009 


GUatG 


8008 


CONTINUE 


TEMPN( I 


T=I+1 


= aG ely) 


HEV IeGieN) GOSTG 8007 


GO TO 


8006 


CONTINUE 


T=1 


CONTINUE 


J=1 


CONTINUE 


Q(I,u)= 


K= 1 


0.D0 


CONTINUE 


(at worrangtt | 


= 


tI uu nt 
aaa 97 eine 7D. ait 


re’ 7. c* of. 13 grey i? et 
iy! aa 7 a4) a7 oH igs. Im 


YP erg 
¢ 


a 
' = 
ha Pa : 3 ' a o.1TD 
a, (8a@ iy ay ud. 3} 
: a5 ; +. aaaiee« S15> i wis.) TR 


oitglibd ReSel) SUT Ww Tae 
: 7 


yucla is ba) ahi: (9. Taw) Se 
WG. (27.wei6 °) TAMAOR 
meet bE. JUTE LPR KA Jake 


sare. bai 
Og. one 


ts 7 

Suvi Tl 2 00 

wiwe( else], ~Ayeaei 11O=) : 
‘ / Hee) 


ene oF Doe. TO. * 
C046 nee Qf 00). 


2 AOD €00E 
iLip-= a) _ it 


$08 ar Sy ot ae 


8015 


8013 


8011 


DiGor@ 


8016 
C 
C 
C 


8019 


8020 
C 
C AN 
C 
8021 


8023 
C 


910 


8024 
8025 


Su22 


Toe 


py Nan IN ee SEE 0S, 
+ 

Ke GieN GG snGieco 15 

GO Terns 014 

CONTINUE 

IF (1.EQ.uU) Q(1,U)=Q(1,U)*(1.D0+ALAM) 
Jat] 

To 2GueN GO 1080 i 

GORI OWES 012 

CONTINUE 

T=1+1 

Pea eeobeN GO: TO) 8011 

GOAIG 8070 

CONTINUE 


TRY e1TOSSCEVEAEINBARS EQUATIONS 


CALL SOLV (N1,M1,N,Q,P,TEMPN, TEMPNN, ITEMP,JTEMP, 
&EPSO,ISING) 

GRELSINGTEQ. 0) GO TO) 8076 

GOmheae 807 

CONTINUE 


CAEGUBARE \G)1.(P) 


GTP=0.D0 

31 

CONTINUE 
GTP=GTP+G(1)*P(1) 
T=I+1 

LAGlnG ban) GOaTO) 8020 
GOR tOmec O19 

CONT INUE 


DECHEG Si Ems (Ze Ora (P)) 


T=] 

GTEST=0 

CONTINUE 

IF (DABS(P(1)).GT.GAMMA*DABS(X(1I))) GO TO 8023 
GO TO 8024 

CONTINUE 


GTEST=1 

WRITE (LT,910) 

FORMAT (/' ##### EXCESSIVELY LARGE GRADIENT #####' /) 
GOm 0 Res025 

CONTINUE 

CONTINUE 

T=]+1 

Pru GIaNeGRyatesie GI .0)G0) 10° 8022 
GORION S023 

CONTINUE 

GOR = 8018 


14, private ‘Bat 
Erde ot cinta 
“mnt 
- 4eg3 Ot de 48 
= “QUAL TTAG. 


ewnt hangs QAgvis BvsG2 at at 


imaseseg.d mt 
| r 


cW3TY .aM3ST1 Ae aMa7T. 4.0 4, 1a, t) Td 43 
37) .9M371 Uma WNT .5, tanta t. Oot 
era of o8'0. ea 
UK! 


emiurio STA 1UOJRD 


Od, osaTo 

tel 

2UMI TOD: € 
i7jaeift ThA tite) 
rele 

nta9g “OF o 4h ep bat 7 


4 9 = 


iQ 
(a9 3% sbr2 37 49H. OH a 


> aa 


(\-wenee TAPTGARD Seeds ave 

.cn : oy : 

_ sd, ee 
ee 


. ~ 


a 


C 
C 
C 


8017 CONTINUE 
8018 CONTINUE 


FIRST MAIN CHOICE ON SEARCH TYPE 


PONT SUNG EQS IT ORT GIP. GE. 0, D0,0R%G1ES) .Giv0 G0. 1G 


GOLMaN e027 


8026 CONTINUE 
C 


C 
C 


C 


Ce oECONDSMATNMCHOICE ON SEARCH OT YPE 
C 


DIO 


C 
C 
C 


(OMB) 


C 
C 
C 


DETERMINE GMAX 


GMAX=0.D0 
T=] 

8029 CONTINUE 
TT=DABS(G(1) ) 
IF(TT.GT.GMAX) GO TO 8031 
GOn)OsG8032 

8031 CONTINUE 


GMAX=TT 

LMAX=I1 

GUm Ome S033 
8032 CONTINUE 
8033 CONTINUE 


fin 
EatCinGilenNiGoO 10) 8030 
GONTONUSO29 


8030 CONTINUE 


PEUGMAXBLESEPS1)GO TO 8034 
COMUOmmnoU 3. 
8034 CONTINUE 
PERFORM A GRID SEARCH 


IGRID=1 


8026 


CHEMBGRIDTINI SMA NI MeLeX EXT TE Wer StL NE Lae 


STOP IF NO SOS IMPROVEMENT 


ISTOP=1 
Trestle GEsS) SG0ri0790 


OR GO ON TO ANOTHER ITERATION 


GOSyOeFs6036 
8035 CONTINUE 


PERFORM A LINE SEARCH 


IGRID=0 
T=1 


ato a.m 

£0 tT Opi, 10.1709 

; esos OT 
ay 


: : i | 
scv1 Wokate WO ZOTOHD HIAM 
a (1207, 24, XAMe 

a | 
wonase CISD A Mme 


aie 


(ys, WW, ART2, abe 5 at 1 ER ESM. remo 1199" 


8037 CONTINUE 
Bi S0. BO 
on +4 
IF(I.GT.N)GO TO 8038 
GUM] ORC OS / 
8038 CONTINUE 
P(LMAX) =-G(LMAX) 
CALE ETNERANGENASNGMMVE XOPNSE X Ta El OMST Ley a1. 
&W,NF,LU) 
ALAM=4.D0*ALAM 
8036 CONTINUE 
GUaTLUMFBOZE 
8027 CONTINUE 


C 
CeyrPERFORM*AVBISECTION SEARCH 
C 


IGRID=0 
CABESE LSC Ia NIN GONOM, 0X, PyS, X10 Po ar teste ae 
ye W, BEPAY GTP, NEY LU) 
ALAM=ALAM/4.D0 
8028 CONTINUE 


C 
Une ewes GR TDS SEARCH WAS DONE 
C 


ReelGRUDeE Oe) Goel. 8069 
GO TO 8040 
8039 CONTINUE 
C 


Geo Lene CONVERGENCE STE ST ENG 
C 
GO TO 8041 
8040 CONTINUE 


CHECK CONVERGENCE 


GiGye@ 


TP=0.D0 
T=1 
8042 CONTINUE 
TPaTee On ELL) -ACL) (ATI eh UL 


C= 
Cee GIN GOO. 8043 
GO TO 8042 


8043 CONTINUE 
PRCOSORT (Pee EEREPS2) GOr1Ores044 
GO TO 8045 
8044 CONTINUE 
C 


ICONV=1 
GO TO 8046 
8045 CONTINUE 
TERDABS (STIU-S)LESERSS*S) "GORI OENS047, 
GO TO 8048 
8047 CONTINUE 


155 


/1, 4, 4072, S084, oF 


a v 
dete Ate, StTH, E96 1. Mat ts i JIA 


. 


‘st 2a HonaTa 
£006 a coy 3,038 oa 
ety 
evisTRat BOMBDAAVIRGD. SHE 


Tear al 


C 
C 
C 


ICONV=2 

GO TO 8049 
8048 CONTINUE 

ICONV=0 
8049 CONTINUE 
8046 CONTINUE 


TRY A GRID SEARCH IF CONVERGENCE INDICATED 


IF(ICONV.GT.0)GO TO 8050 
GOMIG 1 8054 


8050 CONTINUE 
C 


2) (GE) QDi@ie> 


Sie @ 


QS 


8053 


ENG han, Gere 8054 
GOm Gre SUs3 
8054 CONTINUE 
GAEUAGR EGA CNG SMI WNeN LA XT ILS Y, DyWy hits ie NEw LU 
ISTOP=2 


STOP IF NO SOS IMPROVEMENT 
Dato wil eGerns | aGO STO +90 
OR GO ON TO ANOTHER ITERATION 


GUS LO Sago SZ 
8051 CONTINUE 
8052 CONTINUE 
8041 CONTINUE 


DO ANOTHER ITERATION IF WITHIN LIMIT 


ES LOR =s 
Dien eT ER. 200 GG TO eto 


STOP PROGRAM WITH APPROPRIATE PRINTOUT 


90 eee lenOa a WRostE: (eu OS Oa heR 
Ree reaaNeve 0 eeWwRH PEt VT) SOO) ML IDER 
950 FORMAT (' ###########% OPTIMIZATION HALTED AFTER’, 
Cael eeeuI ERATIONS | 
SOSom APG ESTOPR.EOs 1) GO 10 8057 
GO 88058 
8057 CONTINUE 
DERG ULGIEO SWRI CRLEUPS52) 
Poe eU WRITE CT 952) 
952 FORMAT (’ NO SOS IMPROVEMENT AFTER ILL-CONDITIONED’ , 
GaeetGRIDSSEARCH’ | 
GOmn0 Med 10S 
S0>S lm IS|OPSE0.2) GO TO 8059 


: 
a 


’ 
_ ns pine 


ayraatent yr03 yas. 14 HOMAgR Salta 
“OT GR 0, 19; 
‘ nae Ts d i 2 Ten 
iTve 


elt : = 
anes, 

bi are “Te 

peo. OT ai af Ae 7 


fotege 
au, J172, 4079. E, EARP Ae aise 2agtet _ 

19 
‘EWAVORONT 202.04 a1 x 29 


56 GY Ge 12.30; ATE). at 


“tf 
inst 


vOLtIARSTI @aNTORA OT HO 0 
secs. OF oo 7 
qn THOS 208 


iUMLTHOD 508 
AVALT 7 evs 


Cidi WEP ie, VM. portanatt ek ad 


rte > 
yt oF Me (Qc S + 


rgoria’ arateagnene HT iw MARDOF 07 


dat 
Te re et 
aa * ee dy Pi in (ua. Site W 


373k ST JAK wnireasear eet bate Tp eo 
¥ og | ae 


a. 
7 


: 


< 
~ 


Di Oi GO @iic 


tes 


GO TO 8060 
8059 CONTINUE 
Ee Gl 70) swe ee GU, 954) = Cony 
PRD LT ONE LUV WRITE (LT 954) ICONV 
954 FORMAT (’ NO SOS IMPROVEMENT FROM A GRID SEARCH’ , 
Cees ShrrERe CONVERGENCE «CHECKS 413} 
GOS Os ate 
8060 CONTINUE 
Ue EUG Oe Wild Bae 956) 
Reeth teNEsbU)«WRiTeeLEE; 956} 
956 FORMAT (’ MAXIMUM ITERATIONS EXCEEDED WITHOUT ‘, 
CBUSEONVERGENGE®)) 
8105 CONTINUE 
NC=M*NF+M*N*ITER 
DEC EUNGT. 0) WRITE (00,958) NC 
Lae OL WeNES Ee) WRIMES@UET *958) ING 
CoC oBGORMAT  fERPFECTIVENESS MEASURE 9157" 91410) 
WRITE (LT,960) 
960 FORMAT (’ ENTER LUN FOR FULL PRINTOUT (O=NONE)’ ) 


LU=5 
PoP ULUAEOr OF ReIuURN 
T=] 
8106 CONTINUE 
AUPeGl) =x 1) 
T=]+] 
HEGreGTAN GO, 10) 8107 
GOMTORGE 106 


8107 CONTINUE 
CAC EANDGRUN TE MIN Wo beT DIY el Wa Rieti odes 
CREMESTATLS BCNGs Ma Ne MS Lala Wy Avil ey eh toilette toes cretar, 
Ge eT EMPNeTEMPNN GE TEMPSU TEMP; EPSOSESTAT AT StATs 
COME PVAR EUR SVE) 
RETURN 
END 
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DOUBLE PRECISION XTIL(N1),F 
roy VALUE 
S=0.D0 
T=1 
8108 CONTINUE 
CALL MODEL (N1,M1,N,M,L,XTIL,T,1, VALUE) 
Py AU aL) 
SESeF Ql) sbi) *wil) 
l=I+1 
IF(1.GT.M)GO TO 8109 
Goa 28108 


5 » 
alee 


TOM IW O2092943 | 


nat l+jteleI ei 
) aTisw iO, Tai 


aa ‘yi eit (U4). 2 Pad 
ioet, 23 2 aa pavi 23943 *hF 


weiteh) TUDTHLAS 5m? S24 ws AaTHS y i Orr 
yagTan. 10.03.02) ae 


ae 


Tora. OF opti. ths es 
Ore nt eae 


40, SETS SETA TAIT 2) 20d JER : 
1.4.0. &, 11Te irra JLEK MT 1 ey! erate. 382 
SAT2TTARET Dees, Seat, NTL MMaNIT “INST Tse 
1 


vex, Ul WAVGAS 
WAUT roe : 
Serene 


-— «a2 2 & °° i-<«e26¢ 62 64° ° * 


arAGEiaTAOr KY = 3 sThudae 
‘1.422 ana 


ow 66 oes O-@ 244288 eee eee” 


(m2, 2,8,7. 4) tGma> JI TUOREUS 
lan Lwin a aan O12 10989. 318U0e 


| ieavaaerbaaa E} 10M 48g, 
; os ' yy : 4h 4 


Y) 
a — he 
ae rE 
‘ 


A 


QDae@® 


C 


8109 CONTINUE 
NF=NF+ 1 


SUBROUTINE AX (NT MI,N,M,L,XTIL, 1 Ay PEMPeN) 
ene PREGESTON WMP) AT REUNG) AMS Nee eA LNGs 
=| 
8110 CONTINUE 
CAE ESMDERS (NI MIS NOME LER tel ste M PN 
J=1 
8112 CONTINUE 
A(1,J)=TEMPN(u) 
Jet] 
EReUsGieN GOs 10" sete 
GOPHO - Satie 
Guo CONT LNGE 
T=1+1 
ERLE Gt NM) Gomi OES 144 
GONTONES 110 
8111 CONTINUE 
RETURN 
END 
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SUBROUDINESSOLY (N1;M1,M 0, Pot tEMP, ERGOWT eROW: 
*EPS,ISING) 

DOUBLE PRECISION Q(N1,N1),P(N1),T(N1), TEMP(N1,N1) 
FEODSORMMDABS, EL, AL, AMULET, EPS 

DIMENSION IROW(N1),JROW(N1) 

ISING=1 


Shar OL VeANUeLROWSTOeZERO 


T=1 
8114 CONTINUE 
P(I)=T(1) 
IROW(1)=0 
T=1+1 
Peel. GLIM)GOSTOR Sis 
GO TO 8114 
Ciao GON INOUE 
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C COPY Q INTO TEMP 
C 


J=1 
8116 CONTINUE 
T=] 
8118 CONTINUE 
TEMP OTs) =O01, J) 


T=I1+1 
LRCL.GieM)GO (C8 8119 
GON) OSS ic 

8119 CONTINUE 
J=adt] 
Pa GieM) GO TOs 8147 
SU ees 6 


8117 CONTINUE 
C 
GC) FOR PEACH TO SIHESM COLUMNS aN TEMP 
C 


JJ=1 
8120 CONTINUE 


C 
Cee NOS ies oRGnoh EEEMENT SiN) THE SCOEUMN 
C 


Ea isis 
T=] 

8122 CONTINUE 
AL=DABS(TEMP(1, Ju 
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eRe WiC) COMO eANDS Al Gib GO 1Ogeei24 


GOaTOe 6125 
8124 CONTINUE 


EL=AL 
IL=I 
GORI GC eo26 
8125 CONTINUE 
8126 CONTINUE 
T=I+1 
Reel Gia GON TUS e312 3 
GUBIOme S122 
8123 CONTINUE 


CHECK ON MAGNITUDE FOR CONDITIONING 
Tee De eee S) RETURN 
SAVE ROW POSITION 


(P)KIKSY) PQ) ONG 


IROW(IL)=uud 
JROW (UU) =IL 


C 
C FORWARD ELIMINATION PHASE 
C 


T=1 
8127 CONTINUE 
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: FOR EACH ROW NOT YET ELIMINATED 
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Giss 


FECEROW( TI SEQsC)GG 70 8129 
COPTER IS {30 
CONTINUE 


CALCULATE A MULTIPLIER 


AMULT=TEMP(1I,JJ)/TEMP(IL, Ju) 
TEMP(1,JuU)=0.D0 


ELIMINATION ON TEMP 


J=Jd+1 

CONTINUE 

Yee aoe metals ee Mele Me Tey 
= U/Gr 

DErUnGlen) GOaTO 6136 

GOR Ol 8432 

CONTINUE 


C 
C ELIMINATION ON P 
C 


8128 


8121 
C 


Pelt =B LOS AMUET +P ie] 
GO TO 8131 

CONTINUE 

CONTINUE 

Gey 

IF(I.GT.M)GO TO 8128 
GUeI0 5127 

CONTINUE 

JJ=uUdt1 

IF(UJJU.GT.M) GO TO 8121 
GO TO 8120 

CONTINUE 


C BACKWARD ELIMINATION FOR EACH COLUMN 
C 


JJ=M 


8134 CONTINUE 


IL=JROW( Ju} 


C 
C FOR EACH ROW DUE TO BE ELIMINATED 
C 


T=] 


8136 CONTINUE 


IF(I.NE.IL.AND.TEMP(I,JJ).NE.0.D0) GO TO 
GOmLO meoacd 


8138 CONTINUE 

C 

Ce CARCURAI EAS MUL GEBEDOR 
C 


AMULT=TEMP(1,Ju)/TEMP(IL, uu) 


8138 
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a TEMP(1I,uUU)=0.D0 
C ELIMINATE ON P 
C 


P(1)=P(1)-AMULT*P(IL) 
GOSTLORtS 1.40 

8139 CONTINUE 

8140 CONTINUE 
T=I+1 
Bri beGte M)GOM107 88137 
GOmrOs 6136 

G13/ CONTINUE 
JJ=uu- 1 
batuuwaes 1) GO TOM 8135 
GUPILGE S134 

8135 CONTINUE 


C 
C NORMALIZATION PHASE 
C 


T=1 
8141 CONTINUE 
IL=ITROW(1) 
Bees) SPH) TEMP (15 11) 
T=]+1 
IF(I.GT.M)GO TO 8142 


8142 CONTINUE 
C 
C REARRANGE ROWS 
C 


J=1 

8143 CONTINUE 
IL=1ROW(d) 
TEMP(IL,1)=P(J) 
J=dt1 
IF(uU.GT.M)GO TO 8144 
GUSTO 8143 

8144 CONTINUE 
J=1 

8145 CONTINUE 
P(J)=TEMP(u, 1) 


J=udt1 
TE(U,G).M)GO 70 8146 
GOaHGe 6145 
8146 CONTINUE 
ISING=0 
RETURN 
END 
C 
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SUBROUTINE GRIDS NGM to NG Malet l teeny alee Weta lie 


SSL Nip Us) 


DOUBLE PRECISION GAMMA(4),X(N),XTIL(N1),Y(M),T(M,L), 


*FTIL(M1),W(M),SMIN,STIL,S 
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DAT AA GAMMAY DOM DOr ONDOS=201D0/ 

Te (LUNGT? OF SWRET ES (1 UPS00) 

FORMAT (' GRID SEARCH: PRINTOUT OF STIL Kees) 
SMIN=1.D+35 

as 
CONTINUE 
PGi) = 
eS leay 
Reet GiheN) Goael Geo dt Sal 
GOe1O==8450 

CONTINUE 

w =f 

CONT INUE 

bev 

CONTINUE 


X(1) 


CALCULATE A SEARCH POINT 


XTIL(1)=XTIL(1)*(1.D0+GAMMA( J) ) 
Tole Gime y A Olah 0 oa) 


EVALUATE THE SOS 


CRIPRRAND Ss (NT SM GN Melee X YT Wok, eo Lee N Ee 
CEPCPURGTSOIMWRITE (LUS O10) “STATIS Kae ND 
FORMAT "5D1575) 


C 
Ge CHECK Tih 
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Bins 


GES ed oN) GOmnO) Gs. 
GOBLOn ci) / 
CONT INUE 


SMIN=ES TIE 

JMIN=d 

IMIN=I 

GO 70a88156 

CONT INUE 

CONTINUE 

[=i 

TECIGGTAN)GO 10 8155 
GOelLOe, 6154 

CONT INUE 
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8153 CONTINUE 


SET UP THE MINIMUM POINT 


920 


XTIL(IMIN)=XTIL(IMIN) *(1.DO+GAMMA(JUMIN) } 
CACERPANDS MEN TAM TENN YL eX PL Yee he W epee S TE NES 
CreevbUS GIO eWRITE SLU, 920) STIL 

FORMAT (‘ MINIMUM STIL IN GRID IS’ ,D15.5) 

Die SEU GReO EWR ITE (CUS 910) (XT GT eaten? 
RETURN 

END 
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SE=0.D0 
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CAEGUUATE 2SF sia (Py) 
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8160 


PTP=0.D0 

ta 

CONTINUE 
PTP=PTP+P(1)*P(1) 

iP ay 

Teele GT oN) GO 10 6160 
GORTOMSE 159 

CONTINUE 
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Ge -CALCUEAT Eo 


C 
C 


H=tS=SEy/PTP/2.D0 
TP=.5D0/DSQRT(PTP) 
Te ERE ale ar) esTkP 


FIRST SEARCH 


C 
8161 CONTINUE 


8163 


SBAR=STIL 

H=2.D0*H 

{i=3] 

CONTINUE 

One A= lee PL) 
Cael aa 

IF(I.GT.N)GO TO 8164 
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GOr 1G -Snigs 


8164 CONTINUE 


900 


910 


100 


C 


CALL FANDS (N1,M1,N,M,L,XTIL,Y,T,W,FTIL, STIL,NF) 

IF (LU.GT.O) WRITE (LU,900) H,STIL 

FORMAT (' FIRST LINEAR SEARCH: H =’ ,D15.5, 
erSOS Es. 4045.5) 

Ge MEU GHEG) WRITE (10,910) (XTILI 1) tet 

FORMAT (5D15.5) 

IF(STIL.GE.SBAR)GO TO 8162 

DO 100 IK=1,N 

sel erantnele 

S=STIL 

GO TO 8161 


8162 CONTINUE 
C 


Gaeoe COND SEARCH 
C 
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8167 


8168 


8169 
8170 
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C 


ITR=1 

CONTINUE 

DT R=20R+4 

ALPHA=H*H*PTP/(STIL-S+2.D0*H*PTP) 

TE CAEPHAS GE. 0. 25D0*H.AND ALPHA. LE 10. 75D0+4H |) GOTO 3167 
GOBIGHT S168 

CONTINUE 


DEL=ALPHA 

IDEL=0 

GO TO 8169 

CONTINUE 

DEL=0.75D0*H 

IF (ALPHA.LT.DEL) DEL=ALPHA 

IF (H/4.D0.GT.DEL) DEL=H/4.D0 

IDEL=1 

CONTINUE 
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CONTINUE 

XTIL(1)=X(1)+DEL*P (1) 

T=I+1 

IF(I.GT.N)GO TO 8171 

GO TO 8170 

CONTINUE 

CALL FANDS (N1,M1,N,M,L,XTIL,Y,7,W,FTIL,STIL,NF) 

Te (8U0Ghld) WRITE (LU{920) DELIsTIO 

FORMAT (’ SECOND LINEAR SEARCH: DEL = ,D15.5, 
me S0Se=" 901505) 

eRe Gia) WRITE GLU o10uen x Uli) pete=ceN) 

CEST S) GO 10 18172 

GO TO 8173 

CONTINUE 

DO 110 IK=1,N 

X(1IK) =XTIL (IK) 

S=STIL 


ICONV=1 
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G73 CONTINUE 
ICONV=0 
TF WIDER. EQraneH=DEM 

8174 CONTINUE 
DEWECONVEEQRIZOR DIR? GItH5) 1GOlTO  AS166 
GORTONTE 165 

8166 CONTINUE 
LEC 1T RE PEs 
DOVS14 751 
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1,N 
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S=STIL 


8148 CONTINUE 
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RETURN 
END 


SUBROUTINESRI SEVEN IYM IGN MEL, Xo Pe Sh Atl ig eile oils 
C i lee bag G Le, NF, LU) 
DOUBLE PRECISION X(N), B(N1) .XTIL(N4) ,ETIL(M1) ,Y(M), 
*T(M,L),W(M), ALPHA, STIL,S,BETA, GTP 
Pee LOsG) 20) sWRIMER (LU, O00) SBETA, GIP 
OO OPP EORN ATS sep SECTION SEARCH WITHS BEI AS=7 9) Di oe 5s 
C meAND eG ie = 6.015. 57) 
ALPHA=2.D0 
8175 CONTINUE 
ALPHA=ALPHA/2.D0 
T=] 
8177 CONTINUE 
KG =X 1) FALPHA*P (1) 


T=]+1 
Pal Glan)! GO tO" 6178 
GO LEME iat 


8178 CONTINUE 
CAGE EANDS SUNT NIN: 18 Se Ys lea Water ieee eee Nee 
is (eehaleny a8) WRITE (L 910) ALPHA, Si au 
910 FORMAT (' ALPHA Tari 5 SOSmaes DiS eon 
Ree eu Ol W Relate (LU, 920) (XTIL (Ted) 
920 FORMAT (5D15.5) 
IF(STIL.LE. (S+2.DO*BETA*ALPHA*GTP)) GO TO 8176 
CORO Masao 
8176 CONTINUE 
RETURN 
END 
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SUBROUTINE SS carom. NIN TaN ON eels Wer el bein Lerso Tei, 
AO, ey DENN | EMPNNG TTEMPSUTEMPREPSO, 
Ron Alo Ge eV AR LU SAS Ven 

DOUBLES PRECISION SV(M) 11M, LE) oWUME XT EEGN Ty E eLeChiia 


*A(M1,N1), 

*CQ(N1,N1),P(N1),TEMPN(N1), TEMPNN(N1,N1),F(M1),XSVE(N1) 
Se OO won Ase EV hoe seo eel Oho O eevee 
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DIMENSION ITEMP(N1),JUTEMP(N1) 

WRITE (5,900) N 

eG RNA acai aes NONLINEAR REGRESSION’ 
PROGRAM 377, sMELHOD OR MEYER ANDS ROLES EROM: 
Petree eh Ch eWOCeER Nahe 97S INONSRG Cre 
mE LHODSeEORSUNCONST RATNED SORT INEZ ALTON, 
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SCORE EGe Cr PARAMETERS ARES, /) 
Wilt] Beets SiO) eee eile Lae Ni) 


FORMAT (6D12.4 

WRITE(5,915) 

FORMAI(//e-" JHE INITIAL PARAMETER VALUES WERE 7 | 
WRibeb EE Cul CAV Elle) Sal be N 

WRITER OC 20) sroiiL 

FORMAT (/’ WEIGHTED SUM-OF-SQUARES 1S‘ ,D12.4,////) 
DOSZO0 5 T= 1M 

2 dbl ah tS ipaldilllio elie) 

WRITE(5,930) 

PORMAIMesPREOLGTEDaN 3 AORUALERT  % 


‘WEIGHTING’ ,’X VALUE’ } 
DUeZoml Kaa 
WRITE (SY S40) PITS YOLKS WOK eel el Kort?) 
FORMAT(11D11.4) 
WRITE(5,950) 
FORMAT (’ 1° ) 
RETURN 
END 
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APPENDIX 4 Results 


This Appendix contains the full experimental results as given by the data analysis 


computer programs. The raw experimental data that was obtained has been summarised 


and only the feed and product mole fractions have been given. There are four sections, 


each of which is preceded by a description of that section. 


i 


The first section contains an example of the full output that was obtained from 
the program EXPTCAL and illustrates the kinetic and thermodynamic data that 
were calculated using this program. 

The next section contains a summary of all the experimental data that were 
obtained throughout the course of this study. 

The third section presents the output from the non-linear least squares program, 
NLSYST using the main program IRANAL. It contains the best fit estimate of the 
parameters of the Tanh curve, A and B, and hence the initial rates. 

The final section contains the output from the non-linear least squares program 
using ADSREG as the main program. ADSREG is a main program similar to 
IRANAL but containing a different function. This function gives the best fit 
estimates of the parameters which mode! the water dependence of the reverse 
Claus reaction. 


A listing of the computer programs that were used is given in Appendix 3. 


a 
7. 
| . = 
ve = et —— : a md er este 
\ ae ewends cae ttt Ge Ganeeny puxe waved? 3 were 2 
+ ent ide Ia noninbrasi slat 
agarcge hemi 1G Gealensse ge 3 “a wanna 69 
nator one (orl aapuo Vall oT Sqr oe encima reitaan 12 eft 
rena 20G0T be 261 OS AS “PD mwtbovg © 
ory we orev wrth n 
ie» latrines eaus ore 4 1@ OuQarnS 3 (ewe naitoes ae 
mure.citt PO apwics et nace vil ident 


b S) ssh) en ov ren) ave @F ty } a noobs @wit ent - 
aa 33 1900 a are tf JAA rege saaicn nena ! 
a 
sum whet ot aerate E tne Crue wn? ov ' owen oa 
" i 2A Gea + wa ti-nge ere mart tues oO Grete pacvae lat RAT 4 
3 n ' i> celtpore sen as ipa YA pried 
awe & a wieuly MOUS RN = 
segs) ofl) Ub Erin aeqatah tataw fF gar (eae elena of! ‘co eran : 
| notnee) aud — 
f miakaeé ry op Son Pe ow ares aaxrce at to gral A 7 


a 
a. 
- : 


OVE Yaa 


Output from EXPTCAL, full experimental results from data set, WI01A 


Pressure and Temperature conditions are those inside the reactor. 


RATE(WS) = total sulphur rate * stoichiometry. 


The volumetric feed is calculated at the reactor conditions. 


The feed partial pressures are those at the reactor conditions. 
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The following reactions have their equilibrium constants designated as follows, 
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2H,0 + 3/2 S, = SO, + 2H,S 

4S,=5S, 

2S5,= 5S, 
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Summary of All Experimental Data 


Explanation of terms, 


SIO4A, sulphur initial rate data set, 0.4% total sulphur. 


WI0O1A, water initial rate data set, 1% water in feed. 


WT6HA, Temperature data set, 6% water, High temperature. 


DO600, Diffusion test data set, impeller soeed 600 rpm. 
PRESS = pressure in atmospheres. 

TEMP = temperature in Keivin. 

WT CAT = weight of catalyst, g. 

S.V. = space velocity, 1/s. 

WT/FA = inverse mole flow rate, g hr/mol. 


CONV = conversion 


The feed and reactor mol fractions of each component are given. 
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Results of Fitting Experimental Data to the Tanh Function 


The identifying code for each data set was illustrated at the beginning of the last 
section. A summary of the experimental conditions at which the data was obtained is 
given. The average percentage deviation from the mean value is given after the '#' sign. 
The factors by which the experimental data have been scaled in order for the program to 
run are noted next. These will affect the actual value of the slope at the origin. The 
estimates of the parameters A and B follow. The initial guess of the parameters is 
recorded since it is necessary to determine that the solution is the global rather than the 
local minimisation of the sums of squares. It should be noted too, that the actual value of 
the weighted sum of squares is dependent on the scaling factors. 

The table that follows lists the experimental data and the least squares estimate of 
the dependent variable. The X values are the scaled inverse mol feed rates. The Y values 


are the scaled conversions 


bi ’ 
vo oulew igus ott PEP SAD CTL: vn ee 
avcait qoepag ant ne:eabnegyly Dachientha.' 
in aeanites an \auce Tree! ary a ke enaeRgNe AMF eI uO? a oat 7 
pniew Yael aura ee? cnn ae wt 


= 
> 


DATA SET S104 


TEMPERATURE 


PRESSURE 


SULPHUR 
WATER 


METHOD OF MEYER AND ROTH FROM THE TEXT OF WOLFE,M.A. 
NUMERICAL METHODS FOR UNCONSTRAINED OPTIMIZATION 
VAN NOSTRAND REINHOLD, 


SCALING FACTORS 
* O. 10000E+01 


CONVERSION 
PARTIAL PRESSURE x 


0. 10000E-01 


oS lezek AT# 50 
1.4767 ATM # 0.1300 
0.17840 % # 4.00000 
Ze O6O200 uae 6.70000 


NONLINEAR REGRESSION PROGRAM 


i973) 


ESTIMATES FOR THE 2 PARAMETERS ARE: 


0.1385E+00 0.5300E+00 


THE INITIAL PARAMETER VALUES WERE: 


0.5000E+00 0O.5000E+00 


WEIGHTED SUM-OF-SQUARES IS 0.1388E-03 


0 
0 


oOocoOoO 


PREDICTED Y 
0 


.4450E -01 
POO On r= Oi 
BO ts Je Oi 
. 9840E -01 
=) 174s) = 618) 
mi oOSE+00 


ACTUAL Y WEIGHTING 
0 


0. 1000E+01 
.4630E-01 0. 1000E+01 
eoSG 1 E=OtrC 1O00RtOI 
.8760E-01 0.1000E+01 
eYO0350E=01 70. 1O00E+01 
. 1090E+00 0. 1000E+01 
.1343E+00 0. 1000E+01 


X VALUE 
0 


.6284E+00 
.8091E+00 
2 OREO 
~16/75E+01 
Pe SASS AO) 
pooOEtO 
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DATA SET SI05 


TEMPERATURE = 592.3 K # 0.1 
= 1.4789 ATM # 0.1500 


PRESSURE 
SULPHUR 0.24325 % # 6.12000 
WATER 2.69740 % # 10.82000 


SCALING FACTORS 
CONVERSION x 0. 10000E+01 
PARTIAL PRESSURE * 0. 10000E-01 


NONLINEAR REGRESSION PROGRAM 
METHOD OF MEYER AND ROTH FROM THE TEXT OF WOLFE,M.A. 
NUMERICAL METHODS FOR UNCONSTRAINED OPTIMIZATION 
VAN NOSTRAND REINHOLD, (1978). 
ESTIMATES FOR THE 2 PARAMETERS ARE: 


0.1370E+00 0.6201E+00 


THE INITIAL PARAMETER VALUES WERE: 
O.5000E+00 0.5000E+00 
WEIGHTED SUM-OF-SQUARES IS 0.4806E-04 


PREDICTED Y ACTUAL Y WEIGHTING X VALUE 
0.0 : 0.1000E+01 0.0 
O26254E-0 1) 0, 6102E-01 0. 1000EF01 90. 7950 E+00 
0.4925E-01 0.5416E-01 0.1000E+01 0.6070E+00 
0.1169E+00 0.1129E+00 0.1000E+01 0.2044E+01 
OPiS hVEFOOSO21S42E+00 OC S1000EtO1 Orsi E+0i 


DATA SET SIO6 
TEMPERATURE = 591.4 K # 0.1 


PRESSURE 1.4813 ATM # 0.1300 
SULPHUR 0.29798 % # 2.21000 
WATER 2.54520 % # 13.87000 


SCALING FACTORS 
CONVERSION x 0. 10000E+01 
PARTIAL PRESSURE * 0.10000E-01 


NONLINEAR REGRESSION PROGRAM 


METHOD OF MEYER AND ROTH FROM THE TEXT OF WOLFE ,M.A. 
NUMERICAL METHODS FOR UNCONSTRAINED OPTIMIZATION 


VAN NOSTRAND REINHOLD, (1978). 


ESTIMATES FOR THE 2 PARAMETERS ARE: 


0.1247E+00 0.8152E+00 


THE INITIAL PARAMETER VALUES WERE: 


0.5000E+00 0.5000E+00 


WEIGHTED SUM-OF-SQUARES IS 0.4640E-04 


PREDICTED Y ACTUAL Y WEIGHTING 

0.0 0.0 0. 1000E+01 

0.3749E-01 0.4059E-01 0. 1000E+01 

O.5554E-01 0.5022E-01 0.1000E+01 

0.6783E-01 0.7074E-01 0.1000E+01 
0. 


OLSS73E-O0f 02 9353E-01 1000E+01 


X VALUE 
0.0 
O.3807E+00 
0.5876E+00 
0.7481E+00 
0.1198E+01 
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DATA SET SI08 


TEMPERATURE = 592.0 K # Q.2 

PRESSURE = 1.4755 ATM # 0.3600 
SULPHUR 0.40130 % # 3.20000 
WATER 2.15840 % # 10.40000 


METHOD OF MEYER AND ROTH FROM THE TEXT OF WOLFE,M.A. 
NUMERICAL METHODS FOR UNCONSTRAINED OPTIMIZATION 


SCALING FACTORS 
CONVERSION 
PARTIAL PRESSURE * 


* O. 10000E+01 
0.10000E-01 


NONLINEAR REGRESSION PROGRAM 


VAN NOSTRAND REINHOLD, 


(1978). 


ESTIMATES FOR THE 2 PARAMETERS ARE: 
OFStO7E Ci 


0.9947E+ 


00 


THE INITIAL PARAMETER VALUES WERE: 
0.5000E+00 0.5000E+00 
WEIGHTED SUM-OF-SQUARES IS 0.2329E-03 


Oooogoo00o0000 


PREDIGIED TY 
.O 


FOZ Ze OT 
. 4607E-01 
Fo590E-01 
so/GIE-OT 
.6794E-01 
LOvgoE-U1 
. 8688E -01 
.8940E-01 
SEI] 2 (07) 
.9106E-01 
. 9069E-01 


Sore ee oe ioe 


ACTUAL Y 
0 


-4017E-01 
.4946E-01 
FO OZE-O1 
.5447E-01 
60 6SE-01 
.8018E-01 
.8744E-01 
.8642E-01 
On OE-01 
.8870E-01 
. 9990E-01 


(Slee (el(ea(e(ele(e\(e)(eie) 


WEIGHTING 


. 1OOOE+01 
. 1000E+01 
. 1000E+01 
. 1000E+01 
. 1000E+01 
. 1000E+01 
. 1000E+01 
. 1000E+01 
. 1000E+01 
. 1000E+01 
. 1000E+01 
. 1000E+01 


X VALUE 
0 


. 3705E+00 
.5601E+00 
. 7189E+00 
. [496E+00 
. 9692E+00 
. 1478E+01 
. 1885E+01 
.2354E+01 
i2338Et01 
. 4695E+01 
eo LOOEFO 1 
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DATA SET S110 


TEMPERATURE = 592.0 K # 3.6 
PRESSURE = 1.4680 ATM # 0.1000 
SULPHUR 0.60076 % # 3.60000 
WATER 2.991208% BF 14. 10000 
SCALING FACTORS 
CONVERSION * 0. 10000E+01 
PARTIAL PRESSURE * 0.10000E-01 


METHOD OF MEYER AND ROTH FROM THE TEXT OF WOLFE,M.A. 
NUMERICAL METHODS FOR UNCONSTRAINED OPTIMIZATION 


NONLINEAR REGRESSION PROGRAM 


VAN NOSTRAND REINHOLD, 


(1978). 


ESTIMATES FOR THE 2 PARAMETERS ARE: 
0.7081E-01 


0. 1658E+ 


01 


THE INITIAL PARAMETER VALUES WERE: 
0.5000E+00 0.5000E+00 
WEIGHTED SUM-OF-SQUARES IS 0.1603E-03 


OOO99990000000 


PRED UGie Dey, 
20 


eSO23E-01 
.3682E-01 
.4523E-01 
-0030E-01 
poo20E 0 
.6412E-01 
. 6964E -01 
. 7043E-01 
SHO T2E-0 1 
SAO SE—U | 
. 7080E -01 
. 7081E-01 


ACTUAL Y WEIGHTING 
0 


0. 1000E+01 
USCE Ole Uo. FOUUEEOT 
.4138E-01 0.1000E+01 
S432 7E-01 Oa1000E+OT 
20607E-01 0. T000E+01 
.5086E-01 0.1000E+01 
,6182E-01 0. 1000E+01 
.6676E-01 0.1000E+01 
-O970E-01 0. TOO0E+0 
.7620E-01 0.1000E+01 
.6969E-01 0.1000E+01 
-fOS0E-01 0. fO00ESOi 
»/250E-01) 0] 1O00E+01 


OOODOD9000000000 


X VALUE 


0 
e275 0E +00 
. 3476E+00 
. 4560E+00 
sJoooEtO0 
. 7012E+00 
- 9058E+00 
. 1443E+01 
.1787E+01 
22o0E*0) 
.2246E+01 
. 2882E+01 
.4437E+01 


sartiut aHT 


oan f' 


‘ 2.0 9+30002.0 
egeaegar.o 2x asthe “90-12 OaTHOLW i 


Bt oe y eatorees 
25 ae 1o-aesae 


DATA SET S115 


TEMPERATURE 


=4 59135) kK 
PRESSURE = 


1.47283 


fal 
ATM # 0.6500 


On90153 % # 6.91000 
238 1009s 6.60000 


SCALING FACTORS 
CONVERSION x 0. 10000E+01 
PARTIAL PRESSURE * 0.10000E-01 


SULPHUR 
WATER 


NONLINEAR REGRESSION PROGRAM 


METHOD OF MEYER AND ROTH FROM THE TEXT OF WOLFE,M.A. 
NUMERICAL METHODS FOR UNCONSTRAINED OPTIMIZATION 


VAN NOSTRAND REINHOLD, 


(O75 le 


ESTIMATES FOR THE 2 PARAMETERS ARE: 


0.5499E-01 


0.2567E+ 


01 


THE INITIAL PARAMETER VALUES WERE: 


0.5000E+00 0.5000E+00 


WEIGHTED SUM-OF-SQUARES IS 0.2129E-03 


PREDICTED Y 


OODODO00000000 
Oo 


Sahetsia (or 
Sao =O 
4443E-01 
pOOSEa | 
ROO OEE 
.D491E-01 
. 5499E-01 
. 5499E-01 
. 5488E -01 
.4750E-01 
. 3842E-01 


ODOOD0O00000000 


ACTUAL Y WEIGHTING 
0 


0. 1000E+01 
.2474E-01 0.1000E+01 
Pouce Ode On OOUREO 
.4405E-01 0.1000E+01 
.4725E-01 0.1000E+01 
.5543E-01 0.1000E+01 
.5864E-01 0.1000E+01 
.5741E-01 0.1000E+01 
POIs tE-O0180.1000E+0) 
5/60E-0 i@0e 1 OO0E+O 
.3790E-01 0.1000E+01 
.4550E-01 0.1000E+01 


X VALUE 
0 


. 1649E+00 
. 2540E+00 
-4366E+00 
.6730E+00 
.8681E+00 
mi Soe O 
.2133E+01 
.2754E+01 
mI sSOE+01 
. 5094E+00 
. 3368E+00 


é M28 
WoiTAas 


shea 


ran 


see 
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DATASET S120 
TEMPERATURE = 591.5 K pO 


PRESSURE 1.4590 ATM # 0.4300 
SULPHUR 1.25030 % # 2.00000 
WATER 2.69050 % # 4.30000 


SCALING FACTORS 
CONVERSION x 0. 10000E+01 
PARTIAL PRESSURE * 0.10000E-01 


NONLINEAR REGRESSION PROGRAM 
METHOD OF MEYER AND ROTH FROM THE TEXT OF WOLFE,M.A. 
NUMERICAL METHODS FOR UNCONSTRAINED OPTIMIZATION 
VAN NOSTRAND REINHOLD, (1978). 
ESTIMATES FOR THE 2 PARAMETERS ARE: 


0.4586E-01 0.3516E+01 


THE INITIAL PARAMETER VALUES WERE: 
0.5000E+00 0.5000E+00 
WEIGHTED SUM-OF-SQUARES IS 0.9622E-06 


PREDICTED Y ACTUAL Y WEIGHTING’ X VALUE 

0.0 GEO 0.1000E+01 0.0 

5 ZO20E-01F Of 2271-01 OF TOCOE- OI Ca t5S5E+00 
0.3194E-01 0.3271E-01 0.1000E+01 0.2447E+00 

O235/709E-O01,0rco7Se-01 70 f000E+OIR Ors OGE+OD 


wOOne a 


ODE A . 
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PAL) 


DATA SET S125 
TEMPERATURE = 592.2 K ie AOE 


PRESSURE 1.4795 ATM # 0.3600 
SULPHUR 136200) 4. 4s 4.10000 
WATER 2.51840 % # 8.00000 


SCALING FACTORS 
CONVERSION x 0. 10000E+01 
PARTIAL PRESSURE * 0.10000E-01 


NONLINEAR REGRESSION PROGRAM 
METHOD OF MEYER AND ROTH FROM THE TEXT OF WOLFE,M.A. 
NUMERICAL METHODS FOR UNCONSTRAINED OPTIMIZATION 
VAN NOSTRAND REINHOLD, (1978). 
ESTIMATES FOR THE 2 PARAMETERS ARE: 


0.3802E-01 0.4459E+01 


THE INITIAL PARAMETER VALUES WERE: 
0.5000E+00 0.5000E+00 
WEIGHTED SUM-OF-SQUARES IS 0.3802E-17 


PREDICTED Y ACTUAL Y WEIGHTING X VALUE 
0.0 0.0 0.1000E+01 0.0 

0.2324E-01 0.2324E-01 0.1000E+01 0. 1594E+00 
O233604E-015073364E-01 OF 100UE+0 1 O2SIS5E+00 


MeRDOAS WGESE BADR Angel wie 
¢ That BT AGO TOM Oak Hava 30 GOHTSM 
nite deat ATE\RCIMD 2" wOONT SM IAD ‘MM 
Giat) .OgePITRA GAARTZON Ms 
72 ee 
an, gaatamise > aT N02 earAM TER -_ 
iaeee.6 r- accne.d 7 


970 Sasa asi TuAAS satTuas HT 
- poegagOe.n se~0008.0 - 
(<a00ge.c 27 eal 0-08 OFTHE - 


[seed | 


7 
sh 4A & sl sao y garore 
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cousin : sag § ioe. cod 
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| 


218 


DATASET Sisd 
TEMPERATURE = 592.2 K # 0.2 


PRESSURE le5o25 AIM 99 #920.0600 
SULPHUR Pto200 % # 0.90000 
WATER 2.21460 % # 13.90000 


SCALING FACTORS 
CONVERSION * 0. 10000E+01 
PARTIAL PRESSURE * 0.10000E-01 


NONLINEAR REGRESSION PROGRAM 
METHOD OF MEYER AND ROTH FROM THE TEXT OF WOLFE,M.A. 
NUMERICAL METHODS FOR UNCONSTRAINED OPTIMIZATION 
VAN NOSTRAND REINHOLD, (1978). 
ESTIMATES FOR THE 2 PARAMETERS ARE: 


Uys SO3E—Onm Oe 5oGDEr04 


THE INITIAL PARAMETER VALUES WERE: 
0.5000E+00 0.5000E+00 
WEIGHTED SUM-OF-SQUARES IS 0.1147E-05 


PREDICTED Y ACTUAL Y WEIGHTING X VALUE 
0.0 0.0 0.1000E+01 0.0 
On2280E-015072239—-01> 0. T000E+0 1h 03 1545E+00 
0.2866E-01 0.2952E-01 0.1000E+01 0.2374E+00 
OfsdiniE-O1 07. 3062E-01 07 1000E+O 10. 3060E+0G 


ie 


= + A, ) 
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DATA SET S135 
TEMPERATURE = 593.2 K # 0.0 


PRESSURE 1.4776 ATM # 0.1000 
SULPHUR 2.00550 % # 2.70000 
WATER 2. Old O vo 7.20000 


SCALING FACTORS 
CONVERSION x 0. 10000E+01 
PARTIAL PRESSURE * 0.10000E-01 


NONLINEAR REGRESSION PROGRAM 
METHOD OF MEYER AND ROTH FROM THE TEXT OF WOLFE,M.A. 
NUMERICAL METHODS FOR UNCONSTRAINED OPTIMIZATION 
VAN NOSTRAND REINHOLD, (1978). 
ESTIMATES FOR THE 2 PARAMETERS ARE: 


O72937E-Ol) 026353E+01 


THE INITIAL PARAMETER VALUES WERE: 
O.5000E+00 0.5000E+00 
WEIGHTED SUM-OF-SQUARES IS 0.8954E-07 


PREDICTED) Y “AGTUAE Y WEITGHTANG Xx] VALUE 
0.0 0.0 0.1000E+01 0.0 
O22192E=01 022182E-01 02 1000E*0 1) 021519E+00 
05 2659E-01950.2679E-01 07 1O00E+01 02 2351E+00 
OnzS18E-01 0. 2803E-01 0. 1000E+0 1) 07 3052E+00 


DATA SET WI01 


TEMPERATURE 


PRESSURE 


SULPHUR 
WATER 


METHOD OF MEYER AND ROTH FROM THE TEXT OF WOLFE,M.A. 
NUMERICAL METHODS FOR UNCONSTRAINED OPTIMIZATION 


OOZE Ze ik 
W422 


0.50240 % 
1.03250 % 


SCALING FACTORS 
CONVERSION 
PARTIAL PRESSURE x 


fe O2 
AT™M # 0.3100 
# 2.10000 
# 1.00000 


7 cE TOOOOE-Oy 


0. 


10000E+01 


NONLINEAR REGRESSION PROGRAM 


VAN NOSTRAND REINHOLD, 


(1978). 


ESTIMATES FOR THE 2 PARAMETERS ARE: 


0.1262E+00 0O.5544E-01 


THE INITIAL PARAMETER VALUES WERE: 


0.5000E+00 0.5000E+00 


WEIGHTED SUM-OF-SQUARES IS 0.2397E-03 


(SOS) (SSeS) 


PREDICTED: Y 
nO 


Ke} | 74010) 
.9439E-01 
. 1150E+00 
he oesO0 
. 1250E+00 
a rZi2e+OO 


ACTUAL Y WEIGHTING 
0 


0.1000E+01 
.8434E-01 0.1000E+01 
295 7-6E— Ob 0; LOO0ESOi 
PiitoerOO 0. TOOOEFOT 
.1208E+00 0. 1000E+01 
.1357E+00 0.1000E+01 
.1113E+00 0.1000E+01 


X VALUE 
Ao 


wis GE+02 
. 1746E+02 
-2169E+O2 
TOOSIE+O? 
.4769E+02 
erste) ts (ASS (02 


ea 


DATA SET WI02 
TEMPERATURE 


PRESSURE 


SULPHUR 
WATER 


METHOD OF MEYER AND ROTH FROM THE TEXT OF WOLFE,M.A. 
NUMERICAL METHODS FOR UNCONSTRAINED OPTIMIZATION 


D922 
1.4708 


0.48652 % 
(2355082 


SCALING FACTORS 
* 0. 10000E+01 
PARTIAL PRESSURE * 


CONVERSION 


pA 
AT™ # 0.3400 
# 8.50000 
# 7.80000 


GF 


10000E+00 


NONLINEAR REGRESSION PROGRAM 


VAN NOSTRAND REINHOLD, 


(1978). 


ESTIMATES FOR THE 2 PARAMETERS ARE: 


0.1218E+00 0.7327E+00 


THE INITIAL PARAMETER VALUES WERE: 


O0.5000E+00 0.5000E+00 


WEIGHTED SUM-OF-SQUARES IS 0.2098E-03 


OODOO0O0000 


PREDICTED Y 
0 


. 5680E-01 
oO 2n sO. 
.9493E-01 
-SMiE=O1 
. 1085E+00 
. 1088E+00 
_ a eelzes a6) 


ACTUAL Y WEIGHTING 
0 


0. 1000E+01 
soa m/E-O1 OS 1000E+01 
pOCOr Ole On 000s G 
.9618E-01 0.1000E+01 
.9950E-01 0. 1000E+01 
.1044E+00 0.1000E+01 
FOo(oe- O10. O00EsOr 
.1243E+00 0. 1000E+01 


OODOOOD0000 


X VALUE 


0 
.6900E+00 
.9844E+00 
.1425E+01 
. 1492E+01 
. 1947E+01 
. 1963E+01 
.2455E+01 


2 « 


DATA SET W104 


TEMPERATURE = 
PRESSURE = 
SULPHUR 

WATER 


SCALING FACTORS 
* 0. 10000E+01 


CONVERSION 


PARTIAL PRESSURE x 


0. 10000E+00 


592.4 K # ae 
1.4721 ATM # 0.3200 
0.49985 % # 1.70000 
3.68195 % # 6. 10000 


NONLINEAR REGRESSION PROGRAM 


METHOD OF MEYER AND ROTH FROM THE TEXT OF WOLFE,M.A. 
NUMERICAL METHODS FOR UNCONSTRAINED OPTIMIZATION 


VAN NOSTRAND 


REINHOLD, 


(1978). 


ESTIMATES FOR THE 2 PARAMETERS ARE: 


OSItS7E+00 


OS 1427E+07 


THE INITIAL PARAMETER VALUES WERE: 


0.5000E+00 


0.5000E+00 


WEITGHIED “SUM-OF-SQUARES #15507 1233E-024 


PREDICTED Y 

0.0 

0.4875E-01 
-6598E-01 


.9807E-01 
. 1097E+00 


0 
0 
0 0 
OFS 98E-01 70; 
0 0 
0 0 


ACTUAL Y WEIGHTING 
0 


: 0 

-OO1ZE-01 0 

POS JOrmU demu 
0 
0 
0 


8381E-01 


.9644E-01 
. 1104E+00 


. 1000E+01 
. 1000E+01 
1000E+01 
. 1000E+01 
. 1000E+01 
. 1OOOE+01 


X VALUE 
0.0 
0.3837E+00 
0.5460E+00 
0.7438E+00 
O.025E+01 
Gels 1E+01 
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DATA SET W106 


TEMPERATURE = 592.2 K # 0.2 

PRESSURE = 1.4637 ATM # 0.5000 
SULPHUR 0.50760 % # 1.50000 
WATER 5.53100 % # 3.80000 


SCALING FACTORS 
CONVERSION + 0. 10000E+01 
PARTIAL PRESSURE * 0. 10000E+01 


NONLINEAR REGRESSION PROGRAM 


METHOD OF MEYER AND ROTH FROM THE TEXT OF WOLFE,M.A. 
NUMERICAL METHODS FOR UNCONSTRAINED OPTIMIZATION 


VAN NOSTRAND REINHOLD, (1978). 


ESTIMATES FOR THE 2 PARAMETERS ARE: 


0.1272E+00 0.1496E+00 


THE INITIAL PARAMETER VALUES WERE: 


0.5000E+00 0.5000E+00 


WEIGHTED SUM-OF-SQUARES IS 0.1837E-04 


PREDIC IED Y AGIUAE.Y (WELGHEING 
0.0 GEG 0. 1000E+01 
0.4572E-01 0.4911E-01 0.1000E+01 
0.6042E-01 0.5878E-01 0. 1000E+01 
0.9588E-01 0.9410E-01 0.1000E+01 
0.1131E+00 0.1142E+00 0. 1000E+01 


X VALUE 
OG 
0.2516E+01 
0.3453E+01 
0.6563E+01 
0.9490E+01 


a20 


DATA SET W108 


TEMPERATURE 


PRESSURE 


SULPHUR 
WATER 


SCALING FACTORS 
CONVERSION 
PARTIAL PRESSURE * 


METHOD OF MEYER AND ROTH FROM THE TEXT OF WOLFE,M.A. 
NUMERICAL METHODS FOR UNCONSTRAINED OPTIMIZATION 
VAN NOSTRAND REINHOLD, 


* 0. 10000E+01 


0. 10000E+01 


S250 Ke # “OR? 
1.4666 ATM # 0.6300 
0.52170 % # 5.10000 
7.35400 % # 5.40000 


NONLINEAR REGRESSION PROGRAM 


(1978). 


ESTIMATES FOR THE 2 PARAMETERS ARE: 


O29271E-01 


0.2442E+00 


THE INITIAL PARAMETER VALUES WERE: 


0.5000E+00 0.5000E+00 


WEIGHTED SUM-OF-SQUARES IS 0.2882E-04 


PREDICTED Y 
0.0 


0.3944E-01 
OL ones =O 
0.6667E-01 
OS 7972E-01 
0.8865E-01 


QO0000 


ACTUAL Y WEIGHTING 
0 


0. 1000E+01 
.4358E-01 0. 1000E+01 
SOSGe-0 | OS 1000EFOt 
.6405E-01 0. 1000E+01 
.7834E-01 0.1000E+01 
.9029E-01 0. 1000E+01 


X VALUE 
0 


. 1860E+01 
e200 EF0) 
3 LOQJE+01 
. 5294E+01 
7 ORLSEtOT 
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DATA SET W112 


TEMPERATURE 


PRESSURE 


SU 
WA 


METHOD OF MEYER AND ROTH FROM THE TEXT OF WOLFE,M.A. 
NUMERICAL METHODS FOR UNCONSTRAINED OPTIMIZATION 


LPHUR 
TER 


SoZzn0 K 
1.4692 


O. 59258 0% 


11.82300 % 


SCALING FACTORS 
* 0. 10000E+01 
PARTIAL PRESSURE * 


CONVERSION 


# Ue 
AT™ # 0.5400 
# 0.40000 
# 3.80000 


GF 


10000E+01 


NONLINEAR REGRESSION PROGRAM 


VAN NOSTRAND REINHOLD, 


(1978). 


ESTIMATES FOR THE 2 PARAMETERS ARE: 


O27 It9S8E+00 (0. 205TE+00 


THE INITIAL PARAMETER VALUES WERE: 


0.5000E+00 0.5000E+00 


WEIGHTED SUM-OF-SQUARES IS 0Q.1657E-04 


PREDICTED Y 
0.0 


0 


le (a )\(eo (es) 


- SO00OSE-01 
eSO95E-01 
poe ooe=O1 
76352E-01 
76790E-01 
-O9S59E-01 


ACTUAL Y WEIGHTING 


0 

Po200E- Ut 
POO onESOt 
.5043E-01 
76362E=01 
s66o2E-U1 
.9014E-01 


0 
0 
0 
0. 
0 
0 
0 


. 1000E+01 
. 1000E+01 
. 1000E+01 
1000E+01 
. 1000E+01 
. 1000E+01 
. 1000E+01 


X VALUE 
0 


ra2OZEtO" 
wi OODOESO 
.2296E+01 
. 2880E+01 
SSE tO 
-4702E+01 
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DATA SET W115 


TEMPERATURE = 592.2 K # 
PRESSURE 2 


D2 
1.4677 ATM # 0.3300 


SULPHUR 0.50290 % # 3.00000 
WATER 14.31800 % # 3.30000 


SCALING FACTORS 
CONVERSION * 0. 10000E+01 
PARTIAL PRESSURE x 0.10000E+01 


NONLINEAR REGRESSION PROGRAM 


METHOD OF MEYER AND ROTH FROM THE TEXT OF WOLFE,M.A. 
NUMERICAL METHODS FOR UNCONSTRAINED OPTIMIZATION 


VAN NOSTRAND REINHOLD, (1978). 


ESTIMATES FOR THE 2 PARAMETERS ARE: 


0.1069E+00 0.2447E+00 


THE INITIAL PARAMETER VALUES WERE: 


0.5000E+00 0.5000E+00 


WEIGHTED SUM-OF-SQUARES IS 0.3402E-04 


PREDICTED Y ACTUAL Y WEIGHTING 
0.0 20 0. 1000E+01 
0.2446E-01 0.2700E-01 0.1000E+01 
0.4759E-01 0.4323E-01 0.1000E+01 
0.5858E-01 0.6148E-01 0.1000E+01 
OF (B4SE-0 1 sO oi E-O1 On 1000Es01 


X VALUE 
0.0 
0.9519E+00 
0. 1956E+01 
Oe2515E+01 
OF 36s 1E+01 


226 


+ ee layer 0 : cae on see 
am er air NO? | 
cgeVASE.0 OO+RR80r.0. 
200-2 MTR GALTON AT 

‘woeapcine. ¢ nena: 


DATA SET W110 


TEMPERATURE 
PRESSURE 


SULPHUR 
WATER 


METHOD OF MEYER AND ROTH FROM THE TEXT OF WOLFE,M.A. 
NUMERICAL METHODS FOR UNCONSTRAINED OPTIMIZATION 
VAN NOSTRAND REINHOLD, 


oo 17cak 
1.4695 


0.51620 % 
9.39510 % 


SCALING FACTORS 
= 0. 10000E+01 


CONVERSION 
PARTIAL PRESSURE x 


#, BOR 
ATM # 0.6300 
# 0.70000 
# 2.60000 


0. 10000E+01 


NONLINEAR REGRESSION PROGRAM 


(1978). 


ESTIMATES FOR THE 2 PARAMETERS ARE: 


0.1883E+00 0.1046E+00 


THE INITIAL PARAMETER VALUES WERE: 


0.5000E+00 0.5000E+00 


WEIGHTED SUM-OF-SQUARES IS 0Q.1086E-03 


OOMOO90O000 


PREDLETED VY 
a0 


. 2864E -01 
POUL ISES01 
Poved ESO 
m26/7E-0)1 
.7412E-01 
. 1092E+00 


ACTUAL Y WEIGHTING 


m0) 

5 SURES 10M 
-4456E-01 
SD LO/ELOT 
BOW Gree ui 
. 7024E-01 
Mu SoetO0 


0 
0 
0 
0. 
0 
0 
0 


. 1000E+01 
. 1000E+01 
. 1000E+01 
1000E+01 
. 1000E+01 
. 1000E+01 
. 1000E+01 


X VALUE 


0.0 


0 
0 


0 
0 
0 
0 


. 1465E+01 
SIS1OEEO1 
7SO02Ez01 
. 3890E+01 
pO VAS |= a0) 
.6330E+01 
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DATA SET WT6L 


TEMPERATURE = 592.2 K 

PRESSURE = 1.4637 AT 
SULPHUR 0.50800 % # 
WATER 5.53190 % # 


SCALING FACTORS 


# 


O22 
M # 0.5000 


1.18000 
3.09000 


CONVERSION * 0. 10000E+01 


PARTIAL PRESSURE * 0.1 


OOOOE+01 


NONLINEAR REGRESSION PROGRAM 


METHOD OF MEYER AND ROTH FROM THE TEXT OF WOLFE,M.A. 
NUMERICAL METHODS FOR UNCONSTRAINED OPTIMIZATION 


VAN NOSTRAND REINHOLD, ( 


1978). 


ESTIMATES FOR THE 2 PARAMETERS ARE: 


0.1271E+00 0.1498E+00 


THE INITIAL PARAMETER VALUES WERE: 


O.5000E+00 0.5000E+00 


WEIGHTED SUM-OF-SQUARES IS 0.1917E-04 


PREDICTED Y ACTUAL Y WEIGHTING 


Oe 31E+00) 01 42E+00) Oe 
OL9596E-0ilF O-S410E=01 O- 
OL6041E-O1F O53 07E-O01 0; 
O24566E-017 OL49TIE-Oir 0: 
0.0 0.0 OF 


1000E+01 
1000E+01 
1000E+01 
1000E+01 
1000E+01 


X VALUE 
0.9480E+01 
0.6570E+01 
0.3450E+01 
0.2510E+01 
0.0 
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DATA SET WT6M 


TEMPERATURE = 624.3 K 
PRESSURE = 1.4670 
SULPHUR 0.38104 % 
WATER 5.89080 % 


SCALING FACTORS 


# 1.0 
AT™ # 0.5000 
# 2.97000 
# 5.28000 


CONVERSION * 0. 10000E+01 
PARTIAL PRESSURE * 0.10000E+01 


NONLINEAR REGRESSION PROGRAM 


METHOD OF MEYER AND ROTH FROM THE TEXT OF WOLFE,M.A. 
NUMERICAL METHODS FOR UNCONSTRAINED OPTIMIZATION 


VAN NOSTRAND REINHOLD, 


(1978). 


ESTIMATES FOR THE 2 PARAMETERS ARE: 


0.2441E+00 0.1346E+00 


THE INITIAL PARAMETER VALUES WERE: 


Q0.5000E+00 0.5000E+00 


WEIGHTED SUM-OF-SQUARES IS 0.2454E-04 


Oe 
.1970E+00 0. 1968E+00 
P1534E+00 021553E+00 
nissiet005021320E200 
.1031E+00 0.1004E+00 
2 293E-07 0. 7616E-0) 


OOdN000 


0 
0 
0 
GF 
0 
0 


PREDICTED Y ACTUAL Y WEIGHTING 
0 0 


. 1000E+01 
= LOOCOESOi 
. 1000E+01 
1000E+01 
. 1O00E+01 
. 1000E+01 


X VALUE 
0 


-6o10E+01 
.S490E+01 
.4570E+01 
POSOUEZOI 
. 2290E+01 
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DATA SET WT6H 


TEMPERATURE = 653.2 K 
PRESSURE = 1.4687 
SULPHUR 0.29560 % 
WATER 5.88570 % 


SCALING FACTORS 


flO 
ATM # 0.5000 
# 1.40000 
# 1.49000 


CONVERSION + 0. 10000E+01 
PARTIAL PRESSURE * 0. 10000E+01 


NONLINEAR REGRESSION PROGRAM 


METHOD OF MEYER AND ROTH FROM THE TEXT OF WOLFE,M.A. 
NUMERICAL METHODS FOR UNCONSTRAINED OPTIMIZATION 


VAN NOSTRAND REINHOLD, 


H197S8))e 


ESTIMATES FOR THE 2 PARAMETERS ARE: 


0.2296E+00 0.2437E+00 


THE INITIAL PARAMETER VALUES WERE: 


0.5000E+CO 0O.5000E+00 


WELGOTEDeSUM-OR-SOUARES si oem Uro 1S5E-03 


PREDICTED Y ACTUAL Y WEIGHTING 


Cee Zone tOORO Zoe t00 
On 1S41E+00)051912E+00 
0.1415E+00 0.1426E+00 
Omi 5EtOORO 1 2Z03E200 
0 0 
0 0 


20 8) 
20a 1ETOOR0n I GoSEs00 


0 
0 
0. 
0 
0 
0 


- 1000E+01 
. 1000E+01 
1000E+01 
. 1\000E+01 
. 1000E+01 
. 1000E+01 


QBN) (ee; 


X VALUE 
.8710E+01 
.4530E+01 
72950E201 
meOZOETO1 


a0 
. 5S00E+01 
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Results of Estimating the Water Dependence of the Reaction 


This section follows the format of the previous one. In this instance, the X values 
are the partial pressures of water. The Y values are the rates of reaction with respect to 
water. The particular form of the rate equation that is being fitted is given after the data 
summary at the top of each page. The parameters, C, N and K are estimated in that order 


by the program. 


jot sana own eee 
; 7 dy: 
y (en ot: aes ores aT Ww so 
jo om ar 
. 2 bath? ort i Ht Rate ats ot 6 
= joer ) protemisad on? * asa esa tal i te saa 


- awyorg 


DATA SET WRATE 


TEMPERATURE 
PRESSURE 


SULPHUR 
WATER 


METHOD OF MEYER AND ROTH FROM THE TEXT OF WOLFE,M.A. 
NUMERICAL METHODS FOR UNCONSTRAINED OPTIMIZATION 
VAN NOSTRAND REINHOLD, 


Ne 


0: 
0. 


5OZeaiaKe. F® O8L70 
1.4684 ATM # 0.0 

50622 % # 2.90000 

0 * # 4.90000 


C#WEEN/ ((14+KeW) «0. 0) 


NONLINEAR REGRESSION PROGRAM 


(1978). 


ESTIMATES FORM THE 2 PARAMETERS “ARE: 


0.1930E+00 0.4885E+00 


THE INITIAL PARAMETER VALUES WERE: 


0.5000E+00 0.5000E+00 


WEIGHTED SUM-OF-SQUARES IS 0.1764E-02 


PREDIC LEDER 
ORL TSAE=0 1 


OOO0O000 


. 1051E+00 
. 1439E+00 
mi fSoE+O0 
. 2028E+00 
.2274E+00 
.2544E+00 
a2 LI 3ee00 


OO0O0O0000 


ACTUAL Y WEIGHTING 


?6996E—011 
.8924E-01 
. 1349E+00 
. 1903E+00 
. 2264E+00 
. 2420E+00 
.2457E+00 
. 26 16E+00 


0 
0 
0 
0 
0. 
0 
0 
0 


. 1000E+01 
. 1000E+01 
. 1000E+01 
. 1OO00E+01 
1000E+01 
. 1000E+01 
. 1000E+01 
. 1000E+01 


ODO9D00000 


X VALUE 


wLOSOETOU 
. 2880E+00 
. 5479E+00 
.8230E+00 
SL NOGEEOR 
. 1398E+01 
SEN 
P21 LEE 


Paes 


Shi 2net 2R4RR9 £ BMT FOF pore 


1936 ‘bday ATARI UATTIN a 7 


Anan > 954080 


ree 


Of iM, > G 
100 a. Tio 
Go* QUES. 9 


ihaaier? oO)! 


Negeee!..6 | 
4 : 
a 


‘ogres 
(Hew cre 


- = - 
(OC Certara 2 7 


a1 ¢OfeerecA BASH JOM 


4? BMDP T Ce Ca EY Se 30 OGHTSN 
BU hs lad fF EQOHTSM. JAI RG 


iv) .0 21 GERAOS+s0-me qsrHotaw 


4 


es 138 QUANTEOH “ay 
00+ JERR. D Go+s0Ee! .0 


_ gdsmpene.o EL, - 


= = 
a Ps sac ° duutoa | ¥ Gatoroas| 


A.D 19-31ETT..0 
ace 228.0 00+3tZ0r oO 
. eet 0 d<iBet ~{ 
aaron 


462006 
De 


(Qe349001 
2 Obes » - 


DATA SET WRATE 


TEMPERATURE 
PRESSURE 


SULPHUR 
WATER 


METHOD OF MEYER AND ROTH FROM THE TEXT OF WOLFE,M.A. 
NUMERICAL METHODS FOR UNCONSTRAINED OPTIMIZATION 
VAN NOSTRAND REINHOLD, 


v 


C*xW*N/ ( (1+K*W) **0.5) 


Deen # Orel O 
1.4684 ATM # 0.100 
0.50622 % # 2.90000 
0.0 if 4.90000 


NONLINEAR REGRESSION PROGRAM 


OIS7a)e 


ESTIMATES FOR THE 3 PARAMETERS ARE: 


0.2839E+00 0.7196E+00 0.1040E+01 


THE INITIAL PARAMETER VALUES WERE: 


0.5000E+00 O0.5000E+00 0.5000E+00 


WEIGHTED SUM-OF-SQUARES IS 0.1032E-02 


PREDICTED Y 


SeeV(SV(eie)(efe) 


.6847E-01 
. 1017E+00 
. 1470E+00 
. 1811E+00 
. 2082E+00 
. 2306E+00 
. 2534E+00 
ecole sOU 


eVel(e\ele)(@)(e@e 


ACTUAL Y WEIGHTING 


. 6996E-01 
.8924E-01 
. 1349E+00 
. 1903E+00 
. 2264E+00 
. 2420E+00 
. 2457E+00 
.2616E+00 


0 
0 
0 
O: 
0 
0 
0 
0 


. 1000E+01 
. 1000E+01 
. 1000E+01 
1000E+01 
. 1000E+01 
. 1000E+01 
. 1000E+01 
. 1000E+01 


Oo000000 


X VALUE 


. 1536E+00 
. 2880E+00 
.5479E+00 
-6200E+00 
7 LObE+O0] 
. 1398E+01 
a O9E+O1 
Aone Od 
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DATA SET WRATE 


TEMPERATURE 
PRESSURE 


SULPHUR 
WATER 


METHOD OF MEYER AND ROTH FROM THE TEXT OF WOLFE,M.A. 
NUMERICAL METHODS FOR UNCONSTRAINED OPTIMIZATION 
VAN NOSTRAND REINHOLD, 


Y 


SO2ea Ke 05170 
1.4684 ATM # 0.100 
0.50622 % # 2.90000 
GaG hb # 4.90000 


C*W**N/ ( (1+K*W) **1.0) 


NONLINEAR REGRESSION PROGRAM 


(1978) 


ESTIMATES FOR THE 3 PARAMETERS ARE: 


0.3557E+00 0.8790E+00 0.7461E+00 


THE INITIAL PARAMETER VALUES WERE: 


0.5000E+00 0.5000E+00 0.5000E+00 


WEIGHTED SUM-OF-SQUARES IS 0.6907E-03 


PREDICTED Y 


ODOO0000 


.6149E-01 
.9803E-01 
. 1488E+00 
. 1857E+00 
.2129E+00 
p2ooTEtOU 
72027E+00 
.2671E+00 


ACTUAL Y WEIGHTING 


. 6996E -01 
.8924E-01 
. 1349E+00 
. 1903E+00 
. 2264E+00 
. 2420E+00 
. 2457E+00 
. 26 16E+00 


O.1000E+01 
0. 1000E+01 
0. 1000E+01 
0. 1000E+01 
0. 
0 
0 
0 


1000E+01 


. 1000E+01 
. 1000E+01 
. 1000E+01 


X VALUE 


tosOE+O0 
. 2880E+00 
. 5479E+00 
.8230E+00 
. 1106E+01 
. 1398E+01 
. IVOSEFO 
2 181E+01 
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DATA SET WRATE 


TEMPERATURE = 592.1 K # LO 

PRESSURE = 1.4684 ATM # 0.100 

SULPHUR OF 50622) # 2.90000 

WATER OFC % # 4.90000 
Y = C*W**N/ ( (14+K*W) **2.0) 


NONLINEAR REGRESSION PROGRAM 
METHOD OF MEYER AND ROTH FROM THE TEXT OF WOLFE,M.A. 
NUMERICAL METHODS FOR UNCONSTRAINED OPTIMIZATION 
VAN NOSTRAND REINHOLD, (1978). 
ESTIMATES FOR THE 3 PARAMETERS ARE: 


O.3358E+00 0.8911E+00 0.2742E+00 


THE INITIAL PARAMETER VALUES WERE: 
0.5000E+00 0O.5000E+00 0.5000E+00 
WEIGHTED SUM-OF-SQUARES IS 0.5449E-03 


PREDICTED Y ACTUAL Y WEIGHTING X VALUE 

085824E=0190° 6G996E-01 Omf000E+Ods 0. 1536E200 
0.9513E-01 0.8924E-01 0.1000E+01 0.2880E+00 
0.1485E+00 0.1349E+00 0.1000E+01 0.5479E+00 
0.1879E+00 0.1903E+00 0.1000E+01 0.8230E+00 
0.2163E+00 0.2264E+00 0.1000E+01 0.1106E+01 
0.2365E+00 0.2420E+00 0.1000E+01 0.1398E+01 
0.2528E+00 0.2457E+00 0.1000E+01 0.1759E+01 
0.2625E+00 0.2616E+00 0.1000E+01 0.2131E+01 
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DATA SET WRATE 


TEMPERATURE 


PRESSURE 


SULPHUR 
WATER 


METHOD OF MEYER AND ROTH FROM THE TEXT OF WOLFE,M.A. 
NUMERICAL METHODS FOR UNCONSTRAINED OPTIMIZATION 
VAN NOSTRAND REINHOLD, 


uf 


Sek Fe (0). Ae) 
1.4684 ATM # 0.100 
OF 506223 out 2.90000 
0.0 % # 4.90000 


CxW**N/ ( (1+K*W) **3.0) 


NONLINEAR REGRESSION PROGRAM 


(1978). 


ESTIMATES FOR THE 3 PARAMETERS ARE: 


OF S2106+00N025751E+ 00m 0m 5/5E+00 


THE INITIAL PARAMETER VALUES WERE: 


0.5000E+00 0O.5000E+00 0.5000E+00 


WETGHTED SUM-OF-SQUARES IS 0.51/70E-03 


PREDICTED Y 


BD NVe(Hjeo(ye 


. 5808E -01 
. 9468E -01 
. 1481E+00 
AI SSlE+00 
HoOSEFOO 
Zo 2e+00 
. 2530E+00 
726 15E=00 


ACTUAL Y WEIGHTING 


TOOYOEZO 
.8924E-01 
. 134SE+00 
LOU SE LOG 
. 2264E+00 
. 2420E+00 
. 2457E+00 
*26016E+06 


0 
0) 
0 
0 
0. 
0 
0 
0 


. 1000E+01 
. 1000E+01 
. 1000E+01 
. 1000E+01 
1000E+01 
- 1OO00E+01 
. 1000E+01 
. 1O00E+01 


QV qaye 


X VALUE 


. 1536E+00 
. 2880E+00 
.5479E+00 
.8230E+00 
eTTOGE+OR 
.1398E+01 
.1759E+01 
P21G1E+01 
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